Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

順序回帰(ordered regression)

概要

Ordered Regression は、順序尺度(ordinal scale)のカテゴリカルデータを目的変数にとる回帰モデル。

例えば目的変数 yy が:

評価数値
悪い1
普通2
良い3
とても良い4

のように

  • 大小関係はある(順序性あり)

  • 間隔は等しいとは限らない(距離不明)

というカテゴリ変数である場合、

  • 線形回帰(間隔尺度を仮定)

  • 多項ロジット(順序性を無視)

は不適切となるため、専用のモデルが必要となる。

色々なモデル

順序ロジットモデル(orderd logit moodel)順序プロビットモデル(ordered probit model) などがある。

ロジット関数(ロジスティック関数)を使うタイプであってもいくつかタイプがある様子

1️⃣ 累積ロジットモデル系 (まとめて 順序ロジットモデル(orderd logit moodel)順序ロジスティック回帰(ordered logistic regression) と呼ばれることも)

  • 比例オッズモデル(proportional odds model) :ロジスティック回帰モデルの差で確率を表現するモデル P(y=k)=P(yk)P(yk1)P(y = k) = P(y \le k) -P(y \le k-1)

  • 累積ロジットモデル(cumulative logit model) :比例オッズモデルの回帰係数がカテゴリごとに異なるよう一般化したモデル P(ykx)=11+exp[(τkxβk)]P(y \le k \mid \mathbf{x}) = \frac{1}{1+\exp\left[-(\tau_k - \mathbf{x}^\top \beta_k)\right]}

logP(ykx)P(y>kx)=τkxβ(k=1,,K1)\log \frac{P(y \le k \mid \mathbf{x})} {P(y > k \mid \mathbf{x})} = \tau_k - \mathbf{x}^\top \beta \quad (k=1,\dots,K-1)

2️⃣ 隣接カテゴリ・ロジット・モデル (adjacent category logit model)

logP(y=kx)P(y=k+1x)=αkxβ(k=1,,K1)\log \frac{P(y = k \mid \mathbf{x})}{P(y = k+1 \mid \mathbf{x})} = \alpha_k - \mathbf{x}^\top \boldsymbol{\beta} \quad (k = 1, \dots, K-1)

両者の比較

累積ロジットモデル隣接カテゴリロジットモデル
ロジットlogP(yk)P(y>k)\log \frac{P(y \le k)}{P(y > k)}logP(y=k)P(y=k+1)\log \frac{P(y = k)}{P(y = k+1)}
解釈kk以下」vs「kkより上」の累積オッズkk」vs「k+1k+1」の隣接オッズ
確率の導出累積確率の差再帰的に計算
比例オッズの意味累積オッズ比が xx に依らず一定隣接オッズ比が xx に依らず一定

隣接カテゴリロジットモデルの方が「直前のカテゴリから次のカテゴリへの移行」を直接モデル化するため、各カテゴリ間の遷移の難しさを表現する場合に直感的である。

関連文献

References
  1. Wurm, M. J., Rathouz, P. J., & Hanlon, B. M. (2021). Regularized Ordinal Regression and the ordinalNet R Package. Journal of Statistical Software, 99(6). 10.18637/jss.v099.i06