質問: 閲覧ありがとうございます。 Rおよび統計を始めたての者です(RStudioを使用しています)。 Rの学習で使っているサンプルコードの中に、 何の処理なのかが分からない部分があるため、質問させていただきます。 「**下のコードの★行が何を求める処理なのか**」をぜひ教えていただきたいです。 (#~は、各行で行われていることを私が解釈してメモしたものです。) 何卒よろしくお願いいたします。 ``` set.seed(1) #乱数の種を指定 n=300 #データの個数 data=rgamma(n,3,2) #k=3,θ=2のガンマ分布に従う乱数(データ)をn個生成 fit0=fitdist(data,"gamma") #生成したデータをガンマ分布に当てはめる fit0 #最尤推定法によりパラメータ(k,θ(=shape,rate)が推定される fit0$estimate[1]+qnorm(0.975)*fit0$sd[1]*c(-1,1) #★ #fit0$estimate[1]でパラメータkの推定 #qnorm(0.975)は、-∞からの累積確率が97.5%になる座標(ヒストグラムの横軸の座標) #ガンマ分布なのにqnormとなっているのはサンプルコードの誤植? ```
回答: ★部分は、**「真のパラメータを含むであろう区間」を信頼係数95%で推定する処理**であることが分かりました。 また、「ガンマ分布なのにqnormとなっているのはサンプルコードの誤植?」...ではありませんでした。 ``` fit0=fitdist(data,"gamma") fit0 #最尤推定法によりパラメータが推定される(点推定) #↓区間推定 fit0$estimate[1]+qnorm(0.975)*fit0$sd[1]*c(-1,1) #★近似的な信頼区間:推定値±z*標準誤差 #fit0$estimate[1]:パラメータ(k)の推定値 #c(-1,1)で± #qnorm(0.975)は、標準正規分布において、-∞からの累積確率が97.5%になる(=下側97.5%の)点z ```