CCRモデルのソースを表示

'''【しーしーあーるもでる (CCR (Charnes, Cooper and Rhodes) model)】'''
=== 概要 ===

最初のDEA(包絡分析法)のモデルで, チャーンズ・クーパー・ローズにより提案され, 3人の頭文字をとって名づけられたモデルである. 仮想的出力/仮想的入力が1以下の制約下で評価対象事業体にとって最も大きくなるようにする比率形式モデルをもとに, それを線形計画法で記述した各種のモデルがある. 線形計画法で記述したモデルでは主問題, 双対問題それぞれに解釈がなされ, また, 入力改善を指向するか, 出力改善を指向するかでも定式化は異なる.
=== 詳説 ===
　[[DEA]](包絡分析法)のモデルとしてCharnes, Cooper and Rhodesにより提案され, 3人の頭文字をとって名づけられたモデルである[1]. 

　<math>n\, </math> 個の事業体(DMU)に関する<math>m\, </math>個の入力データ<math>X=(x_{ij})\in \mathbf{\mathrm{R}}^{m\times n}\, </math>と<math>s\, </math>個の出力データ<math>Y=(y_{ij})\in \mathbf{\mathrm{R}}^{s\times n}\, </math>をもとに着目DMU <math>J (=1,2,...,n)\, </math>の効率性を測定するために仮想的出力/仮想的入力に着目したモデルに始まり, 多くのモデルが提案された. 

【[[比率形式モデル]]CCR-IR (Input-oriented Ratioform)】


<center><math> \mbox{max.} \;\; D_{J}=\sum_{r=1}^{s} u_{r}y_{rJ}/\sum_{i=1}^{m} v_{i}x_{iJ}\, </math></center>


＜目的関数の解釈：DMU <math>J\, </math> にとって最も有利となるようにウェィト<math>v_i\, </math>, <math>u_r\, </math>を決める. ＞


<table align="center">
<tr>
<td><math>\mbox{s. t.}\, </math></td>
<td></td>
<td></td>
</tr>
<tr>
<td></td>
<td colspan="2">
<math>\sum_{r=1}^{s} u_{r}y_{rj}/\sum_{i=1}^{m} v_{i}x_{ij}\leq 1 \; (j=1, 2, \ldots ,n),\, </math> </td>
</tr>
<tr>
<td colspan="3"><math>u_{r} \geq 0 \ (r=1,2, \ldots ,s) \ ;\; v_{i} \geq 0 \ (i=1, 2, \ldots ,m).\, </math>
</td>
</tr>
</table>


＜制約の解釈：どのDMUの効率値も1以下＞

このモデルの線形計画法(Linear Programming; LP)による定式化は以下のようになる. 

【同値なLP問題<math>CCR_D\, </math>-I：[[入力指向型モデル]]】


<table align="center">
<tr>
<td width="50"><math>\mbox{max.} \, </math></td>
	<td><math>\sum_{r=1}^{s} u_{r}y_{rJ} \, </math></td>
</tr>
<tr>
<td><math>\mbox{s. t.} \, </math></td>
	<td><math>\sum_{i=1}^{m} v_{i}x_{iJ}=1, \, </math></td>
</tr>
<tr>
<td></td>
	<td><math>\sum_{i=1}^{m} v_{i}x_{ij}-\sum_{r=1}^{s} u_{r}y_{rj}\geq 0 \ (j=1,2, \ldots ,n), \, </math></td>
</tr>
<tr>
<td></td>
	<td><math>u_{r} \geq 0 \ (r=1,2, \ldots ,s),\ v_{i} \geq 0 \ (i=1,2,\ldots ,m). \, </math></td>
</tr>
</table>



　このモデルは「乗数」と呼ばれるウェイト<math>v_{i}\, </math>, <math>u_{r}\, </math>を用いていることから乗数形式モデルとも呼ばれる. (乗数については [2] では無限小正数<math>\varepsilon\, </math>以上という制約を課しているが,  ここでは [6] の付録Bの主張に従い, 非負制約のみとした. )

【同値なLP双対問題<math>CCR_P\, </math>-I：入力指向型】


<table align="center">
<tr>
<td width="50"><math>\mbox{min.} \, </math></td>
	<td><math>\theta_{J} \, </math></td>
</tr>
<tr>
<td><math>\mbox{s. t.} \, </math></td>
	<td><math>\theta_{J}x_{iJ}-\sum_{j=1}^{n} \lambda_{j}x_{ij} \geq 0 \ (i=1,2,\ldots ,m), \, </math></td>
</tr>
<tr>
<td></td>
	<td><math>y_{rJ}-\sum_{j=1}^{n} \lambda_{j}y_{rj} \leq 0 \ (r=1,2, \ldots ,s), \, </math></td>
</tr>
<tr>
<td></td>
	<td><math>\lambda_{j} \geq 0 \; (j=1,2, \ldots ,n). \, </math> </td>
</tr>
</table>


　これらのモデルは入力の改善に着目しているので[[入力指向型モデル]]と呼ばれるが, [[出力指向型モデル]]も同様に考えられる. 

　モデルCCR-IRまたは<math>CCR_D\, </math>-IはDMU <math>J\, </math>にとって最も有利な乗数<math>v_{i}\, </math>, <math>u_r\, </math>を求めることを意味する. そのため1項目でも誰にも負けない項目があれば, その項目だけで評価すれば効率値を1にできるので, 一芸入試的評価も可能となる. 

　モデル<math>CCR_P\, </math>-Iの制約は <math>\theta_{J}x_{iJ}, y_{rJ}\, </math> が効率的フロンティアに包みこまれることを意味し, これがData Envelopment Analysis(包絡分析法)の由来となっている (すなわち, [2] ではモデル<math>CCR_P\, </math>-Iを主問題と捕らえている：<math>CCR_P\, </math>-I の下付きのp).  DMU <math>J \;\;  (=1,2,\ldots n)\, </math> が効率的となるためには<math>\theta_{J}=1\, </math> であるばかりでなく,  モデル<math>CCR_P\, </math>-Iの制約におけるスラック変数で表現される入力の余剰<math>\boldsymbol {s}_x=   (s_{x1},s_{x2},\ldots , s_{xm})\, </math>,出力の不足<math>\boldsymbol {s}_y= (s_{y1},s_{y2},\ldots , s_{ym})\, </math> :


<center><math>s_{xi}=\theta_{J}x_{iJ}-\sum_{j=1}^{n} \lambda_{j}x_{ij} \, </math></center>

<center><math>s_{yr}=\sum_{j=1}^{n} \lambda_{j}y_{rj} - y_{rJ}\, </math></center> 


を解消しなければならない. そこで, 


<center><math>\theta_{J}=1 \ ; \boldsymbol {s}_x=\mathbf{0} \ ; \boldsymbol {s}_y=\mathbf{0}
\, </math></center>


のとき, DMU <math>J\, </math> は効率的といわれる. 

　DMU <math>J\, </math> が非効率的なときに, その参照集合の要素数が<math>(s+m-1)\, </math>ならば, DMU<math> J\,</math>は「自然に包絡されている」といわれる. そのときには各入力を<math>\theta_{J}\; (<1)\, </math>倍すればDMU <math>J\, </math> は効率的になる. このことは, CCR入力指向型モデルで考えている効率性が原点とDMU間の距離([[動径距離]] radial metric)で測られていること, すなわち比例的効率性を測られていることを示している. DMU <math>J\, </math> が自然に包絡されていないときには入力の余剰や出力の不足があるにも拘らず<math>\theta_{J}=1\, </math>となっている場合があり, <math>\theta_{J}\, </math>は効率性の適切な尺度になっていない. そこで[[CFA]] (constrained facet analysis) [4] ではDMU <math>J\, </math>に最も近い包絡面上のファセットを延長して1でない効率性尺度が求められるよう工夫している. 

　[[比率形式モデル]] CCR-IRから得られる[[乗数形式モデル (DEAの)|乗数形式モデル]] <math>CCR_D\, </math>-Iの入出力間の制約を凸錐(convex cone)に拡張した [[コーンレシオモデル]] [3] や乗数<math>v_{i}\, </math>, <math>u_r\, </math>の値域に制約を課す[[領域限定法]] [5] なども提案されている. 



----
'''参考文献'''

[1] A. Charnes, W. W. Cooper and E. Rhodes, "Measuring Efficiency of Decision Making Units," ''European Journal of Operational Research'', '''2''' (1978), 429-444.

[2] A. Charnes, W. W. Cooper, A. Y. Lewin and L. M. Seiford, ''Data Envelopment Analysis'' : ''Theory, Methodology and Applications'', Kluwer Academic Publishers, 1994. 刀根薫, 上田徹 監訳, 『経営効率評価ハンドブック』, 朝倉書店, 2000. 

[3] A. Charnes, W. W. Cooper, Q. L. Wei and Z. M. Huang, "Cone Ratio Data Envelopment Analysis and Multi-Objective Programming," ''International Journal of Systems Sciences'', '''20''' (1989), 1099-1118. 

[4] A. Bessent, W. Bessent, J. Elam and T. Clark, "Efficiency Frontier Determination by Constrained Facet Analysis," ''Operations Research'', '''36''' (1988), 785-796.  

[5] R. Allen, A. Athanassopoulos, R. G. Dyson and E. Thanassoulis, "Weights Restrictions and Value Judgements in Data Envelopment Analysis: Evolution, Development and Future Directions," ''Annals of Operations Research'', '''73''' (1997), 13-34.  

[6] 刀根薫, 『経営効率性の測定と改善―包絡分析法DEAによる』, 日科技連, 1993.

[[category:DEA(包絡分析法)|しーしーあーるもでる]]