How convert q-value to p-value?
0
0
Entering edit mode
5.1 years ago
star ▴ 270

I have a vector as below , I like to convert q-value to p-value (estimate p-value base on q-value).

data<-c(0.35,0.52,0.12,0.46)


I used below code for converting but I am not sure about it:(I take this code from here)

qvalues<-as.numeric(as.character(data))

convert.qval.pval = function(qvalues) {
pi0 = max(qvalues)
m0 = length(qvalues) * pi0
return(qvalues * rank(qvalues) / m0)
}

qval_pval <- convert.qval.pval(qvalues)
qval_pval
[1] 0.33653846 1.00000000 0.05769231 0.66346154


Also, I like use "q-value" R package and but for this package I have to have a simulated or data-resampled (e.g., bootstrap, permutation) "null statistics", but I do not know how can I make it.

Code for converting using Q_value:

stat <- data
stat0 <- "creat using bootstar" ?!!!
p.testspecific <- empPvals(stat=state, stat0=stat0, pool=FALSE)

R p-value q-value • 3.3k views
0
Entering edit mode

What is the question? you are not happy with the conversion or you want to use q-value package with your own data?