私は現在、ベルヌーイ データのシミュレーション研究を実行して、サンプル サイズが大きい場合、サンプル比率 ˆp もほぼ正規分布することを示しています。
演習から、ベルヌーイ データは次のように生成されると言われています。
rbinom(n=50, size=1, prob=0.5)
サンプルサイズを 50 として、真の比率を 0.5 として使用しているためです。
指数データを使用したシミュレーション スタディのコードが提供されましたが、指数データ コードを使用せずに、上記のベルヌーイ データのコードを使用してコードを変更する必要があります。サンプル コードは次のとおりです。
n = 10 # sample size
set.seed(142981)
simreps = 1000 # simulation replicates
# (just needs to be a big number)
xbar = rep(0, simreps)
for(i in 1:simreps){
xbar[i] = mean(rexp(n, rate=0.001))
}
hist(xbar, xlim=c(0,3000))
qqnorm(xbar, ylim=c(0,3000));qqline(xbar)
しかし、発生する問題は、次のコードを追加する場合です。
set.seed(142981)
simreps = 1000 # simulation replicates
# (just needs to be a big number)
xbar = rep(0, simreps)
for(i in 1:simreps){
xbar[i] = mean(rbinom(n=50, size=1, prob=0.5))
}
hist(xbar)
qqnorm(xbar);qqline(xbar)
このような値を持つヒストグラムの結果が得られます
しかし、QQ プロットの場合、結果は のよう
になり、正しくないようです。何が悪かったのかわかりません。どんな助けでも大歓迎です。