0

こんにちは、次の 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)
}
4

1 に答える 1