こんにちは、次の for ループをベクトル化することで (つまりsapply()
、 などを使用して) 計算を高速化しようとしています。
しかし、非常に複雑な構造のため、ミッション・インポッシブルのようです。私は数日間試していますが、解決策はありません。
次の関数を高速化する方法を誰か見てもらえますか??
alpha2f <- function(p0,t2,n1,n2){
alpha2=0
for (x1 in 12:n1){
for (x2 in 0:n2){
for (y2 in 0:n2){
a=dbinom(x1,n1,p0)
b=dbinom(x2,n2,p0)
alpha2=alpha2+a*b*dbinom(y2,n2,p0)*ifelse(ztest(x1+x2,y2,n1+n2,n2)>t2,1,0)
}}}
return(alpha2)
}