統計ソフトRで回帰分析のestimateを対数から真数にするには?

にほんブログ村 科学ブログ 数学へ

統計ソフトRの回帰分析のestimateが自然対数の場合、 そのまま使うわけにもいかずどうすりゃいいのさと困る。

真数に変えてくれる方法は探しても見つからなかったので、 自前で関数を作った。

rr.95ci <- function(x){
 res <- exp(cbind(coef(x),confint(x)))
 colnames(res) <- c("RR","LL(2.5%)","UL(97.5%)")
 round(res,2)
}

 

rr.95ci()という関数。

()内にはglm()でlogistic modelやlog-linear modelで 計算した結果のobjectを入れる。

Objectからcoefficientsとconfidence intervalsを抜き出して、 exponentialして返す関数。

小数点以下の桁数は最後のround()内の数字で変えることができる。

これ一つあるだけで、だいぶ楽。