cloglogロジスティック回帰の推定値の解釈


21

cloglogリンクを使用してロジスティック回帰からの推定値を解釈する方法について誰かにアドバイスしてもらえますか?

私は次のモデルを装着しましたlme4:

glm(cbind(dead, live) ~ time + factor(temp) * biomass,
    data=mussel, family=binomial(link=cloglog))

たとえば、時間の推定値は0.015です。単位時間あたりの死亡率にexp(0.015)= 1.015113(単位時間あたり〜1.5%増加)を掛けると言うのは正しいですか?
言い換えれば、loglogロジスティック回帰の場合と同様に、loglogで得られた推定値はlogオッズで表されますか?


コードを編集して、R構文規則に従ってください。持てない( '
— フランクハレル

元の投稿を編集し、コメントを削除します。
— フランクハレル

回答:


30

相補的なログとログのリンク関数では、ロジスティック回帰ではありません。「ロジスティック」という用語はロジットリンクを意味します。もちろん、まだ二項回帰です。

時間の推定値は0.015です。単位時間あたりの死亡率のオッズにexp(0.015)= 1.015113(単位時間あたり〜1.5%増加)を掛けると言うのは正しいですか?

いいえ、log-oddの観点からモデル化されていないためです。それはあなたがロジットリンクで持っているものです。log-oddsの観点から機能するモデルが必要な場合は、logit-linkを使用します。

補完ログログリンク機能は、

η(バツ)=ログ⁡(−ログ⁡(1−πバツ))=バツβ

πバツ=P(Y=1|バツ=バツ)

exp⁡(η)exp⁡(η)=−ログ⁡(1−πバツ)

exp⁡(−exp⁡(η))=(1−πバツ)1−exp⁡(−exp⁡(η))=πバツバツ

バツ


ベンがコメントで彼の質問をやさしくほのめかしたように:

単位時間あたりの死亡率(ハザード)が1.5%増加すると言うのは本当ですか?

相補対数-対数モデルのパラメーターには、ハザード比の観点からきちんと解釈されています。それがあります:

eη(バツ)=−ログ⁡(1−πバツ)=−ログ⁡(Sバツ)S

(例では、ログサバイバルは単位時間あたり約1.5%低下します。)

h(バツ)=−ddバツログ⁡(Sバツ)=ddバツeη(バツ)

P(Y=1)


7
単位時間あたりの死亡率(ハザード)が1.5%増加すると言うのは本当ですか?
— ベンボルカー
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.