順序付きロジスティック回帰の負の係数


17

我々は序数応答があるとy:{Bad, Neutral, Good}{1,2,3}と変数の集合X:=[x1,x2,x3]我々は考えることを説明するy。次に、y(応答)に対してX(設計行列)の順序付きロジスティック回帰を行います。y

推定係数と仮定x1、それが呼び出しβ 1である、順序付けられたロジスティック回帰に- 0.5e - 0.5 = 0.607のオッズ比(OR)をどのように解釈しますか?β^10.5e0.5=0.607

私が1つの単位増加のために」と言うか、paribusをceteris、観察のオッズ良いがある0.607観察の倍のオッズ悪いニュートラルを、とで同じ変更のためのx 1、観察のオッズ中性良いがある0.607 Badを観察する確率は?x1Good0.607BadNeutralx1NeutralGood0.607Bad

私は教科書やグーグルで負の係数の解釈の例を見つけることができません。


2
はい、それは正しいです。正の係数を解釈する方法とほぼ同じです。
ピーターフロム-モニカの復職

2
NB:通常、「Xにを回帰」と言いますが、その逆ではありません。yX
GUNG -復活モニカ

回答:


25

順調に進んでいますが、使用しているソフトウェアのドキュメントを常に参照して、実際にどのモデルが適合するかを確認してください。順序付けられたカテゴリ1 g kおよび予測子X 1X jX pを持つカテゴリ従属変数状況を想定します。Y1,,g,,kX1,,Xj,,Xp

「インザワイルド」では、異なる暗黙のパラメーターの意味を持つ理論的な比例オッズモデルを記述するための3つの同等の選択肢があります。

  1. logit(p(Yg))=lnp(Yg)p(Y>g)=β0g+β1X1++βpXp(g=1,,k1)
  2. logit(p(Yg))=lnp(Yg)p(Y>g)=β0g(β1X1++βpXp)(g=1,,k1)
  3. logit(p(Yg))=lnp(Yg)p(Y<g)=β0g+β1X1++βpXp(g=2,,k)

(モデル1及び2は、内という制限を有する別個のバイナリロジスティック回帰、βのJで変化しないG、及びβは、0 1 < ... < β 0 G < ... < β 0 K - 1、モデル3を有しています約同じ制限β jを、その必要β 0 2 > ... > β 0 gで > ... > β 0 Kk1βjgβ01<<β0g<<β0k1βjβ02>>β0g>>β0k

  • モデル1において、正予測子の増加を意味Xのjはのための増加オッズに関連付けられているでカテゴリYβjXjY
  • モデル1はやや直感に反するため、ソフトウェアではモデル2または3が推奨されるようです。ここでは、正のの予測の増加という手段Xのjはのために増加したオッズに関連付けられているでカテゴリYβjXjY
  • モデルの同じ見積りの1と2のリードが、そのための見積りβ jは反対の符号を有します。β0gβj
  • モデルの同じ見積りの2と3のリード、しかし、そのための見積りβ 0 gが反対の符号を有します。βjβ0g

ソフトウェアがモデル2または3を使用していると仮定すると、「が1単位増加すると、セトリスパリバス、「Y = 良い」と「Y = ニュートラルOR悪い」を観測する予測オッズは、E β 1 = 0.607。「と同様に」中1つの単位増加とX 1、paribusをceteris、予測『観察のオッズY = 良いかニュートラルを『観察対』Y = 悪いの要因によって』変更を電子βX1Y=GoodY=Neutral OR Badeβ^1=0.607X1Y=Good OR NeutralY=Bad。 "経験的なケースでは、実際のオッズではなく、予測オッズしかありません。eβ^1=0.607

以下に、カテゴリーのモデル1の追加図を示します。まず、比例オッズの累積ロジットの線形モデルの仮定。第二に、ほとんどのカテゴリgで観測される暗黙の確率。確率は、同じ形状のロジスティック関数に従います。 k=4genter image description here

カテゴリ確率自体については、描かれたモデルは次の順序付けられた関数を意味します。 enter image description here

PS私の知る限り、モデル2はSPSS MASS::polr()およびR関数およびで使用されordinal::clm()ます。モデル3はR関数rms::lrm()とで使用されますVGAM::vglm()。残念ながら、SASとStataについては知りません。


@Harokittyバイナリロジスティック回帰モデルには、線形回帰モデルのようなエラー項がありません。従属変数自体ではなく、確率をモデリングしていることに注意してください。の誤差分布に関する仮定Yたとえば、Rでを指定する必要がありますglm(..., family=binomial)
カラカル

3つの選択肢のリストで仕様#2を表現する方法を扱ったリファレンスがありますか?

1
@Harokitty Agrestiの「順序カテゴリデータの分析」、セクション3.2.2、p49、式3.8で簡単に説明されています。または、Agrestiの「Categorical Data Analysis」、セクション9.4、p323、式9.12で。
カラカル

Hi, sorry to bother you, do you have a reference for the 3rd one? Agresti doesn't seem to talk about that.

2
@Jase Well, Agresti just uses logit(Y>g) in the section linked above. For logit(Yg), see Harrell's "Regression Modeling Strategies", section 13.3.1, p333, eqn 13.4.
caracal
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.