Rの、★部分のコードの意味（＝何の処理をしているのか）を知りたい

閲覧ありがとうございます。
Rおよび統計を始めたての者です（RStudioを使用しています）。

Rの学習で使っているサンプルコードの中に、
何の処理なのかが分からない部分があるため、質問させていただきます。

「下のコードの★行が何を求める処理なのか」をぜひ教えていただきたいです。
（#～は、各行で行われていることを私が解釈してメモしたものです。）

何卒よろしくお願いいたします。

R
1set.seed(1)　#乱数の種を指定
2n=300　#データの個数
3data=rgamma(n,3,2)　#k=3,θ=2のガンマ分布に従う乱数（データ）をn個生成
4
5fit0=fitdist(data,"gamma")　#生成したデータをガンマ分布に当てはめる
6fit0 #最尤推定法によりパラメータ(k,θ(=shape,rate)が推定される
7
8fit0$estimate[1]+qnorm(0.975)*fit0$sd[1]*c(-1,1) #★
9　#fit0$estimate[1]でパラメータkの推定
10　#qnorm(0.975)は、-∞からの累積確率が97.5%になる座標（ヒストグラムの横軸の座標）
11　#ガンマ分布なのにqnormとなっているのはサンプルコードの誤植？
12

行動規範の内容に同意します

回答1件

自己解決

★部分は、「真のパラメータを含むであろう区間」を信頼係数95%で推定する処理であることが分かりました。
また、「ガンマ分布なのにqnormとなっているのはサンプルコードの誤植？」…ではありませんでした。

R
1fit0=fitdist(data,"gamma")
2fit0 #最尤推定法によりパラメータが推定される（点推定）
3
4#↓区間推定
5fit0$estimate[1]+qnorm(0.975)*fit0$sd[1]*c(-1,1) #★近似的な信頼区間：推定値±z*標準誤差
6　#fit0$estimate[1]：パラメータ(k)の推定値
7　#c(-1,1)で±
8　#qnorm(0.975)は、標準正規分布において、-∞からの累積確率が97.5%になる（＝下側97.5%の）点z