r - R代码临床试验

Question

这是我运行临床试验的代码，以显示试验成功的可能性。我的问题是我需要通过引入第二组样本（n.2）来证明，需要多少样本才能产生高于 90% 阈值的值。请提供任何帮助，我知道我需要循环我拥有的代码，但在这样做时遇到了麻烦。

calc.quant = function( n, X.1, a, b, n.2, nsim, thr, p1=0.025, p2=0.975 )
{   
  a.star = a + n
  b.star = b + n - X.1
  theta = rbeta( nsim, a.star, b.star 
  X.2 = rbinom( nsim, n.2, theta )

  theta.p1p2 = matrix( 0, nrow=nsim, ncol=2 )
  for( j in 1:nsim ) {
    theta.p1p2[j,] = qbeta( c( p1, p2 ), a.star + X.2[j], b.star + n.2 - X.2[j] )
  }

  return( theta.p1p2 )
}

n = 117
X.1 = 110
a = 1
b = 1
n.2 = 50
nsim = 1000
thr = .90

res = calc.quant( n, X.1, a, b, n.2, nsim, thr )

sum( res[,1] > thr ) / nsim

Answer 1

[这不是一个完整的答案，只是为了澄清 OP 的目的。]

带有 for 循环的基本策略：

threshold <- somevalue
for(i in someseq){
    output <- somefunction(...)
    if(output > threshold)
        break
}
output

带有while循环的基本策略：

threshold <- somevalue
below.threshold <- TRUE
while(below.threshold){
    output <- somefunction(...)
    if(output > threshold)
        below.threshold <- FALSE
}