10 求教高手指点,用R语言画的火山图报错

1、library(ggplot2)

2、biowolf<-read.table(file="uu.txt",header=TRUE,row.names = 1)

3、data$sig[(data$padj > 0.05|data$padj=="NA")|(data$logFC < 2)& data$logFC > -2] <- "no",data$sig[data$padj <= 0.05 & data$logFC >= 2] <- "up",data$sig[data$padj <= 0.05 & data$logFC <= -2] <- "down"

4、volcano<-ggplot(biowolf,aes(x= log2(foldchange), y=-1*log10(pvalues)))

报错截图如下
attachments-2018-06-9fpVZpp85b344b7a07560.png
 我想把fold change大于-2和小于2,以及padj > 0.05 的基因都作为不关注的基因。 
请先 登录 后评论

2 个回答

调研图

第三部改成这样:
data$sig[(data$padj > 0.05|data$padj=="NA")|(data$logFC < 2)& data$logFC > -2] <- "no"
data$sig[data$padj <= 0.05 & data$logFC >= 2] <- "up"
data$sig[data$padj <= 0.05 & data$logFC <= -2] <- "down"

请先 登录 后评论
shawn

R语言中 两句话不应该写在一行

如果一定要写在一行 用英文分号间隔 ;

如:

> a = 1; b = 2

> a

[1] 1

> b

[1] 2

>

所以你那么分成两行 要么’up‘后面用 分号 ;

请先 登录 后评论