ポアソン回帰GLMで９５％信頼区間を求める。

Question

ポアソン回帰、対数リンク関数のGLMで、ある病気の発症率について分析しています。

説明変数はage, area, sexでいずれもカテゴリカル変数。
ageは４水準、areaも４水準、sexは２水準です。
発症率=件数(number)/人口(population)が関心のある対象なので、
目的変数をnumberとし、populationをオフセット項として用いています。

やりたいことは、
年齢(age)のグループ”2”、
エリア(area)のグループ”3”、
性別(sex)のグループ”1”
に当てはまる人の発症率の推定値とその９５％信頼区間を求めることです。

ポアソン回帰GLMとして、

model = glm(number ~ offset(log(population)) + age + area + sex, 
                       family = poisson, data = lung)

とモデルを作りました。

上記のグループに当てはまる人の発症率の推定値は、

exp(model$coefficients[1]   #Intercept  
   + model$coefficients[2]    #age2      
   + model$coefficients[6]    #area3   
   + model$coefficients[8])   #sex1   

上のコードで求まると思うのですが、

この９５％信頼区間を求めるために、どうすればよいのかわかりません。
confintを使うのだとは思うのですが。
教えていただけるとありがたいです。

Answer 1

求めたい対象である
年齢(age)のグループ”2”、
エリア(area)のグループ”3”、
性別(sex)のグループ”1”　をベースラインに再設定しました。
そのインターセプトの信頼区間が、求めたい信頼区間となります。