Q&Aポアソン回帰、対数リンク関数のGLMで、ある病気の発症率について分析しています。
説明変数はage, area, sexでいずれもカテゴリカル変数。
ageは４水準、areaも４水準、sexは２水準です。
発症率=件数(number)/人口(population)が関心のある対象なので、
目的変数をnumberとし、populationをオフセット項として用いています。
やりたいことは、
年齢(age)のグループ”2”、
エリア(area)のグループ”3”、
性別(sex)のグループ”1”
に当てはまる人の発症率の推定値とその９５％信頼区間を求めることです。
ポアソン回帰GLMとして、
R
1model = glm(number ~ offset(log(population)) + age + area + sex, 2 family = poisson, data = lung)
とモデルを作りました。
上記のグループに当てはまる人の発症率の推定値は、
R
1exp(model$coefficients[1] #Intercept 2 + model$coefficients[2] #age2 3 + model$coefficients[6] #area3 4 + model$coefficients[8]) #sex1
上のコードで求まると思うのですが、
この９５％信頼区間を求めるために、どうすればよいのかわかりません。
confintを使うのだとは思うのですが。
教えていただけるとありがたいです。
回答1件
0
自己解決
求めたい対象である
年齢(age)のグループ”2”、
エリア(area)のグループ”3”、
性別(sex)のグループ”1” をベースラインに再設定しました。
そのインターセプトの信頼区間が、求めたい信頼区間となります。
投稿2019/11/25 23:52
総合スコア21
あなたの回答
tips
プレビュー
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。