the usage of qvalue

[yueli]

hi,

I just want to import my data in the package qvalue under R environment.

The file name is test_soma.txt.

The data format is just as follow:


0.0642720365
0.0642720365
0.0642720365
0.0642810982
0.0642810982
0.0642810982
0.0643117914
0.0643117914

Please see the attached file.

How can I import my data in qvalue package?

Thanks in advance for great help!


Yue

[blancha]

library(qvalue)

# If the file is still zipped
data <- read.table(gzfile("testis_soma.txt.gz")) 

# If the file is unzipped, use read.table directly instead.
#data <- read.table("testis_soma.txt") 

qobj <- qvalue(p=data$V1)

qvalues <- qobj$qvalues
pi0 <- qobj$pi0
lfdr <- qobj$lfdr

summary(qobj)
hist(qobj)
plot(qobj)

[yueli]

Thanks, blancha!

Your answer really help me. I appreciate!

Yue