0

我有一个向量如下，我喜欢将 q 值转换为 p 值（根据 q 值估计 p 值）。

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

我使用下面的代码进行转换，但我不确定：（我从这里获取此代码“ https://stats.stackexchange.com/questions/51070/how-can-i-convert-aq-value-distribution- to-ap-value-distribution ")

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

此外，我喜欢使用“q-value”R 包，但对于这个包，我必须有一个模拟或数据重采样（例如，引导程序、排列）“空统计”，但我不知道我该怎么做。

使用 Q_value 进行转换的代码：

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

0 回答 0