ネストされたfor-loops in R

Question

ですから、私はRにまったく新しく、ネストされたfor-loopをシミュレーションに実装するのに問題があります。このコードの目的は、2つのベータ版のディストリビューションでは、確率が1/2より大きいすべての結果をカウンターがインクリメントされると仮定しています。次に、このカウンタに基づいて、確率をNで除算することによって計算する（母集団サイズ）。しかし、私はこれをn個の値の範囲に対して実行し、それぞれのn個の値の確率（p）の分布を作成したいと思います。私はこれを実装するのが難しいです。どんな助けでも大歓迎です!! forループの内側のための

コード：

n = 5000 
    counter = 0 
    x = rbeta(n,3,1) 
    y = rbeta(n,3,1) 

    for (i in 1:n) 
    { 
     if((y[i]-x[i] >= 1/2)) 
     { 
     counter = counter +1 
     } 
    } 
    p = counter/n 

Answer 1

あなたがpを計算したい場合は、あなただけのp <- mean(y - x > 1/2)を行うことができます。これは私がやって探していたまさにですnの値の範囲のためにこれを行うには、値n.rangeのベクトルを作成し、

n.range <- seq(100, 5000, 100) 
p <- rep(NA, length(n.range)) 
for(n in n.range){ 
    x <- rbeta(n,3,1) 
    y <- rbeta(n,3,1) 
    p[n == n.range] <- mean(y - x > 1/2) 
} 
hist(p)