誰でもこれについて私を助けることができますか? 「p」が完全に外生的で、一様分布に従うとします。次に、TRUE(= 1)またはFALSE(= 0)のダミーであり、3つの要素(1-3、4-6、7-9、. .., 58-60) の "z" は 0 より大きくなければなりません。
たとえば、{1 0 0 1 1 0 0 0 0 0 1 0...} のような「z」を取得した場合、ループをもう一度繰り返したいと考えています ( sum(z[7:9])=0 であるため)。 {1 1 0 0 0 1 0 1 0 1 0 0...} のような新しい「z」を取得するまで、別の「エラー」を描画します。ここで、3 つの要素ごとのすべての合計は 0 より大きくなります。使用するコードは次のとおりです。次のように。どこが間違っていますか?
set.seed(005)
p<-runif(60, 0, 1)
for (i in 1:20) {
repeat {
error= -0.2*log((1/runif(60, 0, 1))-1) # a random component
z=(p<0.5+error) # TRUE/FALSE condition
z=replace(z, z==TRUE, 1) # replace z to 1 if z is true, else z=0
if (sum(z[(3*i-2):(3*i)])>0) {break}
}
}