0
我使用來自qvalue
包(bioconductor)的qvalue
,並將其應用於來自t.test
的pvalues。然後我將結果寫到一個txt文件中,但是我有兩列,一列用於pvalues,另一列用於qvalues。爲什麼發生這種情況?qvalue結果混亂,是兩列?
PVal<-as.matrix(Pval)
str(PVal)
qobj <- qvalue(PVal)
qwrite(qobj, filename = "my-qvalue-results.txt")
@ LyzandeR謝謝您的回覆,我想上面的代碼,但IAM越來越錯誤?錯誤:估計的pi0 <= 0.檢查您是否有有效的p值或使用其他lambda方法。 – Azize
@ LyzandeR沒有上面的代碼作爲測試它產生的錯誤:錯誤:估計pi0 <= 0.檢查您有有效的p值或使用另一個lambda方法 – Azize
這很奇怪。 'set.seed(10)'產生一個等於零的值:P。現在再試一次。 'runif'產生0到1之間的隨機值。現在用'set.seed(9)'嘗試它,它會起作用。 – LyzandeR