1. 程式人生 > >qvalue: Check that you have valid p-values or use a different range of lambda

qvalue: Check that you have valid p-values or use a different range of lambda

have lambda ... support 重現 style range cond port

ERROR: The estimated pi0 <= 0. Check that you have valid p-values or use a different range of lambda.

重現錯誤的代碼:

ps <- runif(3e5)
library(qvalue)
ps <- ps[ps < 0.75]
qs <- qvalue(ps)

Error in pi0est(p, ...) :
ERROR: The estimated pi0 <= 0. Check that you have valid p-values or use a different range of lambda.


seq(0.05, 0.75, 0.05)
[1] 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 0.65 0.70 0.75
qs <- qvalue(ps, lambda = seq(0.05, 0.75, 0.05))

From:

https://support.bioconductor.org/p/74637/

qvalue: Check that you have valid p-values or use a different range of lambda