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

[Это не полный ответ, а просто, чтобы получить разъяснения по поводу того, что ОП будет для.]

Базовая стратегия с для петли:

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

Базовая стратегия с временной петлей:

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