-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 に答える