【发布时间】:2013-09-27 17:24:35
【问题描述】:
这是我的代码:
A<-read.delim("/Users/macbookair11/Downloads/file.txt",header = T)
names(A)
attach(A)
P<--log10(fdr)
plot(lfc,P,type="p",
xlab="log2fc",
ylab="-log10fdr",
pch=16,
xlim=c(-5,5)
)
abline(h =(-log10(0.01)), untf = FALSE, col="red")
abline(v =(log2(2)), untf = FALSE, col="red")
abline(v =-(log2(2)), untf = FALSE, col="red")
text(lfc,P, labels=Transcript, cex=0.6, pos=4, col="red")
我需要标记的不是所有点,而是那些具有 y>2 和 (x1) 的点。
如何更改我的脚本以使其以正确的方式工作?
【问题讨论】:
-
Thread stackoverflow.com/questions/7611169/… 提供了多种解决方案。也许其中一个或几个也适合您的情况?
-
对不起,我有点误读了这个问题,试试 text(lfc[(lfc>1 | lfc2], P[(lfc>1 | lfc2], 标签=Transcript[(lfc>1 | lfc2], cex=0.6, pos=4, col="red") 。没有实际数据,这是未经测试的。