具有2个向量的for循环在R中计算功率失败

Question

我有2个包含数字的向量，我用来模拟我的研究能力，但在for loop部分一直得到这个错误

Error in pwr.2p2n.test(h, n1 = i, n2 = j, sig.level = 0.05) : 
  number of observations in the first group must be at least 2

如果你能给我一些建议，我将不胜感激

# #示例代码

grp1.n <- seq(30,150,5)  ##group 1, N
grp2.n <- seq(30,150,5)-15  ## group 2, N - 15
h=0.85 #specify large effect size
grp1.length <- length(grp1.n)
grp2.length <- length(grp2.n)
power.holder <- array(numeric(grp1.length*grp2.length), dim=c(grp1.length,grp2.length),dimnames=list(grp1.n,grp2.n))
for (i in 1:grp1.length){
  for (j in 1:grp2.length){
    result.pwr.2p2n.test <- pwr.2p2n.test(h, n1=i, n2=j, sig.level=0.05)
    power.holder[i,j] <- ceiling(result.pwr.2p2n.test$power)
    return(result.pwr.2p2n.test)
  }
}

Answer 1

我不完全确定这是否是你想要的，但我认为它是:

grp1.n <- seq(30,150,5)  ##group 1, N
grp2.n <- seq(30,150,5)-15  ## group 2, N - 15
h=0.85 #specify large effect size
grp1.length <- length(grp1.n)
grp2.length <- length(grp2.n)
power.holder <- array(numeric(grp1.length*grp2.length), dim=c(grp1.length,grp2.length),dimnames=list(grp1.n,grp2.n))
for (i in 1:grp1.length){
  for (j in 1:grp2.length){
    result.pwr.2p2n.test <- pwr.2p2n.test(h, n1=grp1.n[i], n2=grp2.n[j], sig.level=0.05)
    power.holder[i,j] <- ceiling(result.pwr.2p2n.test$power)
    return(power.holder)
  }
}

更改在pwr.2p2n。测试函数以及要返回的对象。

旧:pwr.2p2n.test(h, n1=i, n2=j, sig.level=0.05)
新:pwr.2p2n.test(h, n1=grp1.n[i], n2=grp2.n[j], sig.level=0.05)

注意，你的代码中还缺少一个}括号。