Question

这是我的代码：

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）的点如何更改脚本以使其以正确的方式工作？感谢

Answer 1

无法向您展示如何使用您的数据，因为它不可用;这显示了来自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),
                ")" ) )

标记不是标题上的所有点

1 个答案: