Rにおいてオッズ比の95％信頼区間の求め方教えてください。

Question

odds比の95%信頼区間が出せずにいます。
以下の結果に対して95%信頼区間の出し方を教えていただけると嬉しいです。

　　　　coef　　　 exp(coef) se(coef) z Pr(>|z|)
mean 0.59526718 1.8135154 0.007839948 75.927439 0.000000e+00

パッケージvcdを入れたのですが
表が膨大なのでoddsratio(x, log=FALSE)のxに何を入れたらいいのかがわからず
どなたかわかる方お願いします！！

Answer 1

ロジスティック回帰であれば
log <- glm (formula = y ~ x1 + x2 + x3, family = binomial(link = "logit"), data = df)
summary (log) #p value
exp(coef(log)) #odds ratio
exp(confint(log)) #95% CI

重回帰であれば
mra <- glm(y ~ x1 + x2 + x3, family = gaussian(link = "identity"), data = df)
summary (mra) #p value
coef(mra)) #odds ratio
confint(mra)) #95% CI

でそれぞれ出ました。

Comments

[miono]

ありがとうございます。やっている内容は条件付きロジスティック回帰なので最後にくっつけてみます！

Answer 2

https://oku.edu.mie-u.ac.jp/~okumura/stat/2by2.html
にある通りではないでしょうか。

表が膨大なのでoddsratio(x, log=FALSE)のxに何を入れたらいいのかがわからず

オッズ比を求めたい値を入れるだけだと思います。
「以下の結果」は逆に何を求めたのでしょうか？

Comments

[miono]

ありがとうございます。上記の結果からオッズ比1.81が出ていますが、その９５％信頼区間の出し方がわからないのです。方法あったら教えてください。