に私が1.48を超えると0.67を下回っている拒絶領域にあるどのように私の F値の数を確認するために1000回のシミュレーションをしようとしています。ループ内でインクリメントしますか? R
私はこれを持っているが、彼らが必要として変数は増加しません。最終目標は
は、私も試してみましたBとの合計は見つけることです
for (k in 1:1000){
Adata = rnorm(100, mean = 30, sd = 10)
Bdata = rnorm(100, mean = 45, sd = 10)
f = (sd(Bdata)^2)/(sd(Adata)^2)
if (f > 1.48){
a = 0
a <- a + 1}
if (f < .67){
b = 0
b <- b + 1}
}
a
[1] 1
b
[1] 1
:
for (k in 1000){
Adata = rnorm(100, mean = 30, sd = 10)
Bdata = rnorm(100, mean = 45, sd = 10)
f = (sd(Bdata)^2)/(sd(Adata)^2)
a = f > 1.48
b = f < .67
}
y = sum(a)+sum(b)
y
[1] 0
f 'のは拒絶領域にありますか?
代わりにグローバル変数を作成してみてください。汚いですが、問題を解決します。 –