r - 用于计算 R 中功率的 2 个向量的 for 循环失败

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)
  }
}

标签： rarraysfor-loopvector

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.test 函数以及您要返回的对象中。

旧：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)

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