-1

这是我的代码:

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 和(x<-1 或 x>1)的点。

如何更改我的脚本以使其以正确的方式工作?

4

1 回答 1

0

无法向您展示如何处理您的数据,因为它不可用;这显示了 t 检验中“异常值”的一些条件标记。

x <- rnorm(100)
y<- rchisq(100,1)
plot(x,y, xlim=c(-3.5,3.5))
abline(v= c(-1.96, 1.96), col='red')
abline(h= 3.84, col='red')
text(x[ abs(x) > 2|y > 3.8], y[abs(x) > 2|y > 3.8]+.15, 
        paste0("(",round(x[ abs(x) > 2|y > 3.8],1), " , ", 
                   round(y[ abs(x) > 2|y > 3.8],1),
                ")" ) )
于 2013-09-27T21:53:37.973 回答