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.

離散選択モデル

概要

  • 離散選択問題:意思決定者(例えば消費者や企業)が取りうる行動の選択肢の中から自身の目的関数(例えば効用や利潤)を最大化するような行動を選ぶという意思決定問題

    • ある商品の需要を考えるために、他の製品との代替性を考える

  • 動学的な意思決定を考える場合は「今日この行動をとる」と「今日は見送る(明日以降に行動するという選択肢を残す)」の両者の利得の現在割引価値を比較することになる

  • ここでは、ある製品を購入するかどうかの例で考えていく。

効用関数

消費者iiが製品jjを購入した場合の間接効用UijU_{ij}を、製品から消費者が平均的に得る効用VjV_jと個々人の好みよる効用ϵij\epsilon_{ij}(選好ショック)とに分けて捉える。

Uij=Vj+ϵijU_{ij} = V_j + \epsilon_{ij}

ここで Vj=β0+k=1KβkxjkαpjV_j = \beta^0 + \sum_{k=1}^K \beta^k x^k_j - \alpha p_jで、α,β0,...,βk\alpha, \beta^0, ..., \beta^k はパラメータである。

「購入しない」という選択肢は outside goods と呼ばれ、j=0j=0で表す。

以上をまとめると効用は次のように表される。

効用関数

消費者iiが製品jjを購入した場合の間接効用UijU_{ij}

Uij={β0+k=1Kβkxjkαpj+ϵijj=1,...,Jϵijj=0U_{ij} = \begin{cases} \beta^0 + \sum_{k=1}^K \beta^k x^k_j - \alpha p_j + \epsilon_{ij} & j = 1, ...,J\\ \epsilon_{ij} & j=0 \end{cases}
  • α,β0,...,βk\alpha, \beta^0, ..., \beta^k :パラメータ

  • xjkx_j^k:製品属性

  • ϵij\epsilon_{ij}:その製品に対する個人特有の選好で、特定の確率分布に従う確率変数(誤差項)

離散選択問題

以上の効用関数のもとで、消費者iiは最も高い効用が得られる選択肢did_iを選ぶと考える。

di=argmaxj{1,,J}Uijd_i=\underset{j \in\{1, \ldots, J\}}{\operatorname{argmax}} U_{i j}

分析者の視点ではϵij\epsilon_{ij}は直接観測できず、誤差項のようなものである。

そのランダム性を考慮して確率で考えると、消費者iiが製品jjを選択する確率は、選択肢jjから得られる効用UijU_{ij}が他の選択肢から得られる効用UilU_{il}よりも高くなるような事象が発生する確率

Pr(di=j)=Pr({ϵij}j=0JVj+ϵij>Vl+ϵil,lj)\operatorname{Pr}\left(d_i=j\right)=\operatorname{Pr}\left(\left\{\epsilon_{i j}\right\}_{j=0}^J \mid V_j+\epsilon_{i j}>V_l+\epsilon_{i l}, \forall l \neq j\right)

となる。この確率はϵij\epsilon_{ij}が従うと仮定する確率分布の設定によって変わる。

二項選択モデル

最も単純なケースとして、選択肢の数J=2J=2の場合を考える。各財から得られる効用を、平均的な満足度βj\beta_j、価格pjp_j、ランダムな選好ショックϵij\epsilon_{ij}で表す:

Ui1=β1αp1+ϵi1Ui2=β2αp2+ϵi2\begin{aligned} & U_{i 1}=\beta_1-\alpha p_1+\epsilon_{i 1} \\ & U_{i 2}=\beta_2-\alpha p_2+\epsilon_{i 2} \end{aligned}

このとき、消費者iiがとる選択肢did_iは:

di={1 if β1αp1+ϵi1>β2αp2+ϵi22 if β1αp1+ϵi1β2αp2+ϵi2d_i= \begin{cases}1 & \text { if } \beta_1-\alpha p_1+\epsilon_{i 1}>\beta_2-\alpha p_2+\epsilon_{i 2} \\ 2 & \text { if } \beta_1-\alpha p_1+\epsilon_{i 1} \leq \beta_2-\alpha p_2+\epsilon_{i 2}\end{cases}

したがって購入確率は

Pr(di=1)=Pr(β1αp1+ϵi1>β2αp2+ϵi2)=Pr(ϵi1ϵi2>(β1β2α(p1p2)))\begin{aligned} \operatorname{Pr}\left(d_i=1\right) & =\operatorname{Pr}\left(\beta_1-\alpha p_1+\epsilon_{i 1}>\beta_2-\alpha p_2+\epsilon_{i 2}\right) \\ & =\operatorname{Pr}\big(\epsilon_{i 1}-\epsilon_{i 2}>-\left(\beta_1-\beta_2-\alpha\left(p_1-p_2\right)\right)\big) \end{aligned}

となる。

二項プロビットモデル

ϵij\epsilon_{ij}が独立かつ同一の正規分布N(0,0.5)\mathcal{N}(0, 0.5)に従うと仮定する場合をprobit modelという。
この場合、財1の購入確率は、標準正規分布の累積分布関数Φ()\Phi(\cdot)を使って次のように表すことができる。

Pr(di=1)=Φ(β1β2α(p1p2))\operatorname{Pr}\left(d_i=1\right) = \Phi\left(\beta_1-\beta_2-\alpha\left(p_1-p_2\right)\right)

二項ロジットモデル

選好ショックϵij\epsilon_{ij}が第I種極値分布 P(ϵijx)=F(x)=eexP(\epsilon_{i j} \leq x)=F(x)=e^{-e^{-x}} に従うと仮定する場合、言い換えると選好ショックの差分ϵi1ϵi2\epsilon_{i1} - \epsilon_{i2}が標準ロジスティック分布 F(t)=P(Xt)=11+exp(t)F(t)=P(X \leq t)=\frac{1}{1+\exp (-t)} に従うと仮定する場合、財1の購入確率は次のように表すことができる。

Pr(di=1)=11+exp({β1β2α(p1p2)})=exp(β1β2α(p1p2))1+exp(β1β2α(p1p2))\begin{aligned} \operatorname{Pr}\left(d_i=1\right) & =\frac{1}{1+\exp(-\{ \beta_1-\beta_2-\alpha\left(p_1-p_2\right)\})}\\ & =\frac{\exp \left(\beta_1-\beta_2-\alpha\left(p_1-p_2\right)\right)}{1+\exp \left(\beta_1-\beta_2-\alpha\left(p_1-p_2\right)\right)} \end{aligned}

多項ロジットモデル

二項選択モデルを一般化して多項にした場合、多項プロビットモデルは複雑な数値計算が必要になりあまり用いられない。ここでは多項ロジットモデルについて書く。

選好ショックϵij\epsilon_{ij}が第I種極値分布 P(ϵijx)=F(x)=eexP(\epsilon_{i j} \leq x)=F(x)=e^{-e^{-x}} に従うと仮定する場合、消費者にとっての製品jjの購入確率を以下のように書くことができる。

Pr(di=j)=exp(β0+k=1Kβkxjkαpj)1+l=1Jexp(β0+k=1Kβkxlkαpl)Pr(d_i=j) = \frac {\exp( \beta^0 + \sum_{k=1}^K \beta^k x^k_j - \alpha p_j )} {1 + \sum^J_{l=1} \exp( \beta^0 + \sum_{k=1}^K \beta^k x^k_l - \alpha p_l )}

jjの選択確率は自身の価格pjp_jや製品属性xjkx_j^kのみならず、他の選択肢の価格や製品属性(pl,xlk)(p_l, x_l^k)にも依存する点に注意。

最尤推定

顕示選好(「ある製品の選択肢のもとで、人々が実際にどの選択肢を選んだのか」が選好を反映している、とする考え方)の考え方を反映したものが、データを用いた最尤推定である。

製品の価格や性能などのデータ{pj,xj1,,xjK}j=1J\{p_j, x^1_j,\cdots,x^K_j\}^J_{j=1}と、各消費者がどれを購入したかというデータ{di}i=1N\{d_i\}^N_{i=1}が利用可能であるとする。このとき、多項ロジットモデルから尤度関数を以下のように書くことができる

L(θ)=i=1Nj=0J(Pr(di=j))1(di=j)L(\theta) = \prod_{i=1}^N \prod_{j=0}^J (Pr(d_i=j))^{1(d_i=j)}

この尤度関数をパラメータθ=(α,β0,,βK)\theta=(\alpha, \beta^0, \cdots, \beta^K)について最大化することで消費者の効用関数内のパラメータを得ることができる。

反実仮想シミュレーション

一度関数が推定できれば、例えば値下げ時の効用の変化は効用関数に代入する製品価格を下げればシミュレーションができ、売上高の予測もできる。

期待効用

多項ロジットモデルにおける期待効用は

E{ϵij}j[maxj{0,1,,J}Uij]E_{\{\epsilon_{ij}\}_j} [\max_{j\in \{0,1,\cdots,J\}} U_{ij}]

で、これは「選好へのショックが発生する前の時点で最も高い効用をエられるように選択した場合に得られる効用の期待値」となる。

ロジットモデルでは期待効用は以下のように書くことができる。

E{ϵij}j[maxj{0,1,,J}Uij]=log(j=0Jexp(Vij))+CE_{\{\epsilon_{ij}\}_j} [\max_{j\in \{0,1,\cdots,J\}} U_{ij}] = \log \left( \sum^J_{j=0} \exp(V_{ij}) \right) + C

ここでVij=β0+k=1KβkxjkαpjV_{ij} = \beta^0 + \sum_{k=1}^K \beta^k x^k_j - \alpha p_jである。この右辺の形はログ・サム公式(log-sum formula)と呼ばれる。

ランダム効用モデル

こうした確率的なショック項を入れた離散選択モデルはランダム効用モデルとも呼ばれる

文献

直接効用と間接効用

https://www.i-repository.net/contents/outemon/ir/102/102830312.pdf

直接効用関数 u(x)u(x)

数量と効用の関係を表す。

間接効用関数 v(p,m)v(p, m)

予算制約線の元での効用最大化問題を解いた後の、価格との関係性

maxxu(x1,,xj)s.t. p1x1++pjxj\max_x u(x_1, \cdots, x_j)\\ \text{s.t.} \ p_1x_1 + \cdots + p_jx_j\\

から決定される需要関数

xi=Di(p1,,pj,m)x_i = D_i(p_1, \cdots, p_j, m)

を効用関数に代入すると

u(x1,,xj)=u[D1(p,m),,Dj(p,m)]=v(p,m)u(x_1, \cdots, x_j)\\ =u[D_1(\boldsymbol{p}, m), \cdots, D_j(\boldsymbol{p}, m)]\\ =v(\boldsymbol{p}, m)

となる。これを間接効用関数と呼ぶ。