1

我有以下 pvalues 列表

pval.list <- list(list(a=c(0.05, 0.0001, 0.32, 0.45), b=c(0.1,0.12,0.01,0.06), c=c(0.1,0.12,0.01,0.06), d=c(0.01,0.02,0.03,0.04)),
         list(e=c(0.04, NA, 0.232, 0.245), f=c(0.05, 0.01, 0.22, 0.54), g=c(0.005, 0.1, 0.032, 0.045)),
         list(h=c(0.03, 0.01, NA, 0.4), i=c(0.5, 0.0001, 0.132, 0.045), j=c(0.005, 0.0001, 0.0032, 0.045), k=c(0.5, 0.1, 0.932, 0.545)),
         list(l=c(0.022, NA, 0.32, 0.45), m=c(0.0589, 0.0001, NA, 0.0045)),
         list(n=c(0.051, 0.01, 0.32, 0.45), o=c(0.05, 0.0001, 0.32, 0.45), p=c(0.05, 0.0001, 0.32, 0.45), q=c(0.05, 0.0001, NA, 0.45)),
         list(r=c(NA, 0.001, 0.32, 0.45), s=c(0.05, 0.0001, NA, 0.45), t=c(0.05, 0.0001, 0.32, 0.45)))

我正在尝试将一个函数(见下文)应用于此列表:

Fisher.test <- function(p) {
Xsq <- -2*sum(log(p))
p.val <- 1-pchisq(Xsq, df = 2*length(p))
return(p.val) 
}

按照@G.Grothendieck 的指示,我使用了这个命令:lapply(lapply(pval.list, Reduce, f = cbind), apply, 1, Fisher.test)但是当列表中缺少值时,该方法不起作用。我尝试将其合并na.omit为参数,lapply但它不能解决问题。

我应该在函数中添加na.omit作为参数吗?Fisher.test

谢谢

4

1 回答 1

1

您的 NA 是在sum. 如果你通过参数na.rm=TRUE,你的结果中不会有 NA:

Xsq <- -2*sum(log(p), na.rm=TRUE)

您可以为 设置一个na.rm参数Fisher.test,并从中传递 true 或 false lapply,但是当您总是想排除 NA 时,您应该只在函数中设置参数。

在 pvalues 上调用函数:

lapply(lapply(pval.list, Reduce, f = cbind), apply, 1, Fisher.test)

给出:

[[1]]
[1] 1.953968e-03 2.999509e-05 5.320831e-04 1.339104e-02

[[2]]
[1] 0.0007878665 0.0317662968 0.0457435481 0.1146067577

[[3]]
[1] 8.982382e-03 3.055250e-08 4.719676e-02 5.094879e-02

[[4]]
[1] 0.009911091 0.001021034 0.684618971 0.014584426

[[5]]
[1] 2.357950e-03 6.135981e-11 5.543601e-01 6.038488e-01

[[6]]
[1] 6.235072e-02 3.470933e-09 6.016501e-01 5.708789e-01

现在在计算总和之前删除所有 NA 值。这是做你想做的吗?

于 2013-05-24T13:02:04.420 回答