順調に進んでいますが、使用しているソフトウェアのドキュメントを常に参照して、実際にどのモデルが適合するかを確認してください。順序付けられたカテゴリ1 、… 、g 、… 、kおよび予測子X 1、… 、X j、… 、X pを持つカテゴリ従属変数状況を想定します。Y1,…,g,…,kX1,…,Xj,…,Xp
「インザワイルド」では、異なる暗黙のパラメーターの意味を持つ理論的な比例オッズモデルを記述するための3つの同等の選択肢があります。
- logit(p(Y⩽g))=lnp(Y⩽g)p(Y>g)=β0g+β1X1+⋯+βpXp(g=1,…,k−1)
- logit(p(Y⩽g))=lnp(Y⩽g)p(Y>g)=β0g−(β1X1+⋯+βpXp)(g=1,…,k−1)
- logit(p(Y⩾g))=lnp(Y⩾g)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 K)k−1βjgβ01<…<β0g<…<β0k−1β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=4g
カテゴリ確率自体については、描かれたモデルは次の順序付けられた関数を意味します。
PS私の知る限り、モデル2はSPSS MASS::polr()
およびR関数およびで使用されordinal::clm()
ます。モデル3はR関数rms::lrm()
とで使用されますVGAM::vglm()
。残念ながら、SASとStataについては知りません。