信頼区間

二項

信頼区間というのがある 母分布からデータセットをとって、母分布の平均が知りたいときに、昨日は点推定をしたわけだが、今日は区間推定 確率pで生起するときN回観測して、k回起きたという観測結果を基にpを区間推定する Rのbinomパッケージではそのヘルプ文…

線形回帰の信頼区間

適当に点をとって、信頼区間を計算する関数predict()の出力を確認してみる 線形回帰の結果の中身を覗く lm.out <- lm(y~x) summary(lm.out) > summary(lm.out) Call: lm(formula = y ~ x) Residuals: Min 1Q Median 3Q Max -2.6354 -0.2825 0.1845 0.5036 1.…

サンプル数が多くなるときの漸近的収束

こちらで、二項分布での正確な信頼区間推定を行っている サンプル数の増大により、信頼区間幅がサンプル数の平方根に関して収束していくことを書いている p<- 0.8 trial<- 10*seq(from=1,to=100,by=1) apply(cbind(trial*p,trial*(1-p)),1,binom.test,0.5) a…