离散选择模型.pdf

《离散选择模型.pdf》由会员分享，可在线阅读，更多相关《离散选择模型.pdf（17页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、 离散选择模型离散选择模型 郑安 郑安 南京审计学院 经济学院 南京审计学院 经济学院 在经典计量经济学模型中，被解释变量通常被假定为连续变量。但是，我们在数据分析与建模时，有时会遇到被解释变量在经典计量经济学模型中，被解释变量通常被假定为连续变量。但是，我们在数据分析与建模时，有时会遇到被解释变量 Y 不再是连续变量，而是离散变量。例如：公共交通工具和私人交通工具的选择问题，对某种商品的购买决策问题，求职者对某种职业的选择问题，劳动力迁移决策等，这时我们把因变量离散取值的计量模型称为离散选择模型（不再是连续变量，而是离散变量。例如：公共交通工具和私人交通工具的选择问题，对某种商品的购买决策问

2、题，求职者对某种职业的选择问题，劳动力迁移决策等，这时我们把因变量离散取值的计量模型称为离散选择模型（DCM，discrete choice model），如果因变量只能取两个值，则称之为二元选择模型（），如果因变量只能取两个值，则称之为二元选择模型（binary choice model），如果因变量可取多个数值，则称之为多元选择模型。），如果因变量可取多个数值，则称之为多元选择模型。离散选择模型起源于离散选择模型起源于Fechner于于1860年进行的动物条件二元反射研究。年进行的动物条件二元反射研究。1962年，年，Warner 首次应用于经济领域。首次应用于经济领域。20 世纪世纪 7

3、0 和和 80 年代，离散选择模型普遍应用于经济布局、交通问题、就业问题、购买决策问题等经济决策领域的研究。模型的估计方法主要发展于年代，离散选择模型普遍应用于经济布局、交通问题、就业问题、购买决策问题等经济决策领域的研究。模型的估计方法主要发展于 20 世纪世纪 80 年代初期，远远滞后于模型的应用，并且至今还在不断改进，它属于微观计量经济学年代初期，远远滞后于模型的应用，并且至今还在不断改进，它属于微观计量经济学即研究大量个人、家庭或企业的经济信息，即研究大量个人、家庭或企业的经济信息，McFadden 因为在微观计量经济学领域的贡献而获得因为在微观计量经济学领域的贡献而获得 2000 年

4、诺贝尔经济学奖。年诺贝尔经济学奖。一、离散选择模型内容提要：一、离散选择模型内容提要：1.1 定义：因变量取值是离散的，这类回归模型称为离散选择模型。根据因变量的个数以及因变量取值的不同，离散选择模型有如下若干情形：定义：因变量取值是离散的，这类回归模型称为离散选择模型。根据因变量的个数以及因变量取值的不同，离散选择模型有如下若干情形：（1）若因变量取值仅有两个，称为二元选择模型，包括重复观测值可以得到情况下的二元选择模型和重复观测值不可以得到情况下的二元选择模型。）若因变量取值仅有两个，称为二元选择模型，包括重复观测值可以得到情况下的二元选择模型和重复观测值不可以得到情况下的二元选择模型。（

5、2）若因变量取值两个以上，称为多元选择模型，包括一般多元选择模型和排序多元选择模型。）若因变量取值两个以上，称为多元选择模型，包括一般多元选择模型和排序多元选择模型。1.2 研究内容研究内容 （1）模型的设定）模型的设定（2）最大似然法估计参数；）最大似然法估计参数；（3）模型检验）模型检验（4）自变量的边际效应分析）自变量的边际效应分析 1.3 应用举例应用举例（1）买房决策：买与不买，解释变量假定为收入；对于一个确定收入的市民而言，买房与不买房都有可能，这是随机事件，自然而然的，我们所关心的就是买房的概率如何确定。）买房决策：买与不买，解释变量假定为收入；对于一个确定收入的市民而言，买房与

6、不买房都有可能，这是随机事件，自然而然的，我们所关心的就是买房的概率如何确定。假设假设01y=未买房买房，收入记为，收入记为 income。样本数据：特定地区家庭信息（样本数据：特定地区家庭信息（y，income）（2）研究生录用中二元选择模型：）研究生录用中二元选择模型：01y=未录取录取，解释变量虚拟变量，解释变量虚拟变量 011D=非应届生应届生和分数和分数 score。如何确定一个考生的录取概率？。如何确定一个考生的录取概率？样本采集：所有考生信息样本采集：所有考生信息（y,score,D1）。（3）医学中二元选择模型：心肌梗塞和人体内的两项指标密切相关，即高密度蛋白质含量（）医学中二

7、元选择模型：心肌梗塞和人体内的两项指标密切相关，即高密度蛋白质含量（HDL）和纤维蛋白质含量（）和纤维蛋白质含量（Fib）。）。假定假定 01y=未患心肌梗塞患心肌梗塞 样本数据：所关心的某个群体的信息样本数据：所关心的某个群体的信息(y，HDL，Fib)（4）政府支持率的民意测验）政府支持率的民意测验 对政府的态度对政府的态度012y=支持中立反对，解释变量为市民收入，解释变量为市民收入income。样本数据：某城市市民信息（样本数据：某城市市民信息（y，income）二、问题的提出二、问题的提出 例例1：分析选举中投票者的收入对选举的影响。以投票者的态度（：分析选举中投票者的收入对选举的影

8、响。以投票者的态度（y）为被解释变量，令同意者）为被解释变量，令同意者y=1，反对者，反对者y=0，以投票者的月收入（，以投票者的月收入（x）为解释变量，如下表：）为解释变量，如下表：Y X Y X Y X 0 0 100 0 0 1100 1 1 2100 0 0 200 0 0 1200 1 1 2200 0 0 300 1 1 1300 1 1 2300 0 0 400 0 0 1400 1 1 2400 0 0 500 1 1 1500 1 1 2500 0 0 600 0 0 1600 1 1 2600 0 0 700 1 1 1700 1 1 2700 0 0 800 0 0 18

9、00 1 1 2800 0 0 900 1 1 1900 1 1 2900 0 0 1000 1 1 2000 1 1 3000 能否用经典正态线性回归模型处理呢？能否用经典正态线性回归模型处理呢？12iiiyxu=+2122212min()()()()iiiiiyxxxyyxxyx=经典正态线性回归模型的问题：一般假设经典正态线性回归模型的问题：一般假设u服从正态分布，无异方差，无自相关，无多重共线性下，最小二乘法得到的参数估计值是服从正态分布，无异方差，无自相关，无多重共线性下，最小二乘法得到的参数估计值是BLUE的的（1）12()1(1)0(0)(1)iiE yxP yP yP y=+=

10、+=，而，而12ix+可能超出可能超出0，1范围范围（2）iu非正态分布非正态分布 iu 概率 11y=121ix12(1)iP yx=+10y=12ix 121(1)1iP yx=+（）（3）iu的异方差的异方差 221212()()(1)iiiiiuEuEuxx=+所以用线性回归模型不能得到最优线性无偏估计量所以用线性回归模型不能得到最优线性无偏估计量 三、二元选择模型概述三、二元选择模型概述 3.1 二元选择模型要表达什么？二元选择模型要表达什么？如果我们把个体的选择看作是一次由量变所导致的质变，则二元选择模型要表达的就是“量变与质变的关系”。如果我们把个体的选择看作是一次由量变所导致的

11、质变，则二元选择模型要表达的就是“量变与质变的关系”。具体来讲，假设个体的决策用用因变量为具体来讲，假设个体的决策用用因变量为iy表示，则表示，则iy取值取值01或，分别表示是与否，选择与不选择等，影响决策的指标变量记为，分别表示是与否，选择与不选择等，影响决策的指标变量记为1(,)kxx?。二元选择模型要研究的是指标变量如何影响个体的决策概率？比如，特定收入条件下购房决策概率；若干指标下患病概率等等。因此，二元选择模型实际上就是建立个体决策的概率模型。二元选择模型要研究的是指标变量如何影响个体的决策概率？比如，特定收入条件下购房决策概率；若干指标下患病概率等等。因此，二元选择模型实际上就是建

12、立个体决策的概率模型。3.2 个体按什么进行决策呢？个体按什么进行决策呢？个体决策基于如下思想：个体决策受指标变量的影响，有些指标变量给决策带来正倾向，而有些带来负倾向，如果综合倾向超过某个临界值，则个体做个体决策基于如下思想：个体决策受指标变量的影响，有些指标变量给决策带来正倾向，而有些带来负倾向，如果综合倾向超过某个临界值，则个体做 出“是”的决策，否则做出“否”的决策。概括地讲，就是量变导致质变的思想。具体而言，涉及到三个方面的问题：出“是”的决策，否则做出“否”的决策。概括地讲，就是量变导致质变的思想。具体而言，涉及到三个方面的问题：（1）综合倾向如何表达？）综合倾向如何表达？（2）量

13、变的临界值如何选取？）量变的临界值如何选取？（3）个体决策的概率如何计算？）个体决策的概率如何计算？3.3潜回归潜回归 我们假设存在一个不可观察的潜在变量称为决策倾向，是指标变量的连续性 函 数，记 为我们假设存在一个不可观察的潜在变量称为决策倾向，是指标变量的连续性 函 数，记 为*iy，它 与 指 标 变 量，它 与 指 标 变 量ix之 间 具 有 如 下 线 性 关 系之 间 具 有 如 下 线 性 关 系 *1ikkiiiyxxu=+=+?该方程称为潜回归方程，其中该方程称为潜回归方程，其中iu是随机扰动项，是随机扰动项，1ikixx=?，1k=?3.4 量变临界值选取量变临界值选取

14、 量变到多少时，个体才进行选择呢？不妨设为量变到多少时，个体才进行选择呢？不妨设为0。这样，。这样，iy与与*iy有如下关系：有如下关系：*1000iiiyyy=做出“是”决策做出“否”决策 3.5 随机误差项的分布随机误差项的分布 如果给定随机误差项的分布，即给定如果给定随机误差项的分布，即给定()iuF ，则个体决策概率可完全确定。推导如下：，则个体决策概率可完全确定。推导如下：*(1)(0)()1()iiiiiP yP yP uxFx=*(0)(0)()()iiiiiP yP yP uxFx=四、常见的二元选择模型四、常见的二元选择模型 Probit模型模型 若若(0,1)iuN，则称二

15、元选择模型为，则称二元选择模型为Probit模型，可得决策概率：模型，可得决策概率：(1)1()()iiiP yFxFx=2/21(1)()()2ixziiiP yFxxedz=Logit模型模型 若若1xixeue+，则称二元选择模型为，则称二元选择模型为Logit模型，可得决策概率：模型，可得决策概率：(1)1()()iiiP yFxFx=(1)()1iixiixeP yFxe=+Extreme模型模型 若若1 exp()xiue，则称二元选择模型为，则称二元选择模型为Extreme模型，可得决策概率：模型，可得决策概率：(1)1()exp()ixiiP yFxe=在经济计量应用中，主要使

16、用在经济计量应用中，主要使用Probit模型和模型和Logit模型，尽管二者形式和参数有很大的差异，但在分析解释变量对因变量的边际影响时，对于多数的应用问题，二者的实际差异并不是很大。因此，具体选择哪一个模型，并没有严格的标准。模型，尽管二者形式和参数有很大的差异，但在分析解释变量对因变量的边际影响时，对于多数的应用问题，二者的实际差异并不是很大。因此，具体选择哪一个模型，并没有严格的标准。五、二元选择模型的参数估计五、二元选择模型的参数估计 二元选择模型一般采用极大似然估计。二元选择模型一般采用极大似然估计。假设假设iy 的分布列为：的分布列为：若给定样本若给定样本1(,)iikiy xx?

17、(1,2,)in=?，则样本似然函数为：，则样本似然函数为：11()()1()iinyyiiiLF xF x=则对数似然函数为：则对数似然函数为：1ln()(ln()(1)(1()niiiiiLyF xyF x=+对数似然函数取最大值的一阶条件为：对数似然函数取最大值的一阶条件为：1()()ln(1)()0()1()niiiiiiiif xf xLyyxF xF x=+=上述一阶条件，相当于得到上述一阶条件，相当于得到k个方程组，通常是非线性的。用迭代法进一步计算出个方程组，通常是非线性的。用迭代法进一步计算出估计量，具体计算借助于计量软件处理。估计量，具体计算借助于计量软件处理。六、二元选择

18、模型的参数检验六、二元选择模型的参数检验 6.1单个系数的显著性检验单个系数的显著性检验 一个解释变量（对二元决策的概率）是否有显著性影响的检验，如同正态一个解释变量（对二元决策的概率）是否有显著性影响的检验，如同正态iy 0 1 P 1()iF x ()iF x 线性回归分析的单个系数的检验类似，根据模型中的待估系数与其方差计算线性回归分析的单个系数的检验类似，根据模型中的待估系数与其方差计算z统计量，并检验假设统计量，并检验假设0:0iH=。6.2总体显著性检验总体显著性检验 由于由于Logit模型、模型、Probit模型是非线性的，在同时检验多个系数是否为模型是非线性的，在同时检验多个系

19、数是否为0时，时，F检验不能使用。可以采用下述几种检验方法。检验不能使用。可以采用下述几种检验方法。（1）Wald检验（适用于线性和非线性约束）检验（适用于线性和非线性约束）023:0kH=?检验统计量：检验统计量：1()()()WRqRVRRq=其中，其中，1000kRI=，Rq=是零假设是零假设0H的矩阵表示，的矩阵表示，V是估计系数的协方差矩阵，是估计系数的协方差矩阵，是无约束模型得到的估计值。可以证明，是无约束模型得到的估计值。可以证明，W渐进服从渐进服从2(1)k分布。所以分布。所以W检验只需要估计无约束模型检验只需要估计无约束模型（2）对数似然比检验（只适用于线性约束）对数似然比检

20、验（只适用于线性约束）023:0kH=?检验统计量：检验统计量：2ln()ln()RLRLL=其中，其中，ln()RL是约束模型的最大对数似然函数值，是约束模型的最大对数似然函数值，ln()L是非约束模型的最大对数似然函数值。可以证明，在零假设下，是非约束模型的最大对数似然函数值。可以证明，在零假设下，LR渐进服从渐进服从2(1)k分布。所以分布。所以LR检验既需要估计有约束模型，又需要估计无约束模型检验既需要估计有约束模型，又需要估计无约束模型（3）拉格朗日乘子检验（适用于线性和非线性约束）拉格朗日乘子检验（适用于线性和非线性约束）023:0kH=?检验统计量：检验统计量：LMg Vg=其中

21、，其中，g是将有约束模型得到的参数估计值是将有约束模型得到的参数估计值R代入无约束对数似然函数的一阶导数向量中，代入无约束对数似然函数的一阶导数向量中，V是有约束模型得到的参数估计值是有约束模型得到的参数估计值R的渐进协方差估计。可以证明，的渐进协方差估计。可以证明，LM渐进服从渐进服从2(1)k分布。所以分布。所以LM检验只需要估计有约束模型检验只需要估计有约束模型 6.3拟合优度检验拟合优度检验（1）McFadden R2=0ln1lnLL L是估计模型的最大似然函数值，是估计模型的最大似然函数值，0L是是230k=?时最大似然函数值。时最大似然函数值。该指标随着模型拟合的改善而提高。该指

22、标随着模型拟合的改善而提高。（2）期望）期望-预期表预期表（3）H-L Statistic 6.4异方差问题异方差问题 编程计算。见后编程计算。见后 解决方法：稳健估计解决方法：稳健估计 七、自变量的边际效应分析七、自变量的边际效应分析 7.1 边际效应定义：反映自变量单位变化所引起的因变量的变化量。边际效应定义：反映自变量单位变化所引起的因变量的变化量。7.2 线性回归模型的边际效应：线性回归模型的边际效应：1 1kkiyxxu=+?ix的边际效应：的边际效应：iiyx=，因此回归系数直接解释为该变量的边际效应。，因此回归系数直接解释为该变量的边际效应。7.3 二元选择模型的边际效应：二元选

23、择模型的边际效应：二元选择模型：二元选择模型：(1)()iiEyP yFx=（假设（假设iu的分布是对称的）的分布是对称的）则则ix的边际效应：的边际效应：()()iiiiFxEyfxxx=其中其中()f是密度函数，可见系数不能解释为边际效应，是密度函数，可见系数不能解释为边际效应，ix的边际效应除受系数影响外，还受其它变量取值的影响。的边际效应除受系数影响外，还受其它变量取值的影响。八、八、EVIEWS应用举例（例应用举例（例1）：）：8.1模型输出结果模型输出结果 Dependent Variable:Y Method:ML-Binary Probit(Quadratic hill cli

24、mbing)Date:07/10/09 Time:22:20 Sample:1 30 Included observations:30 Convergence achieved after 5 iterations Covariance matrix computed using second derivatives Variable CoefficientStd.Errorz-StatisticProb.X 0.0030670.0011922.5731210.0101C-4.7538961.892117-2.5124750.0120McFadden R-squared 0.706837 Me

25、an dependent var 0.500000S.D.dependent var 0.508548 S.E.of regression 0.274450 Akaike info criterion 0.539743 Sum squared resid 2.109040Schwarz criterion 0.633156 Log likelihood-6.096147Hannan-Quinn criter.0.569627 Restr.log likelihood-20.79442LR statistic 29.39654 Avg.log likelihood-0.203205Prob(LR

26、 statistic)0.000000 Obs with Dep=0 15 Total obs 30Obs with Dep=1 15 所以，估计的潜回归模型为：所以，估计的潜回归模型为：*i-4.753896+0.003067xiy=8.2模型预测模型预测 Y X Yp Y*Y X Yp Y*0 0 100 4.35E-06-4.4471930 0 1600 0.560939 0.153351 0 0 200 1.73E-05-4.14049 1 1 1700 0.677261 0.460054 0 0 300 6.31E-05-3.8337870 0 1800 0.778387 0.766

27、757 0 0 400 0.000210-3.5270841 1 1900 0.858468 1.07346 0 0 500 0.000640-3.2203811 1 2000 0.916232 1.380163 0 0 600 0.001786-2.9136781 1 2100 0.954185 1.686866 0 0 700 0.004567-2.6069751 1 2200 0.976900 1.993569 0 0 800 0.010716-2.3002721 1 2300 0.989284 2.300272 0 0 900 0.023100-1.9935691 1 2400 0.9

28、95433 2.606975 0 0 1000 0.045815-1.6868661 1 2500 0.998214 2.913678 0 0 1100 0.083768-1.3801631 1 2600 0.999360 3.220381 0 0 1200 0.141532-1.07346 1 1 2700 0.999790 3.527084 1 1 1300 0.221613-0.7667571 1 2800 0.999937 3.833787 0 0 1400 0.322739-0.4600541 1 2900 0.999983 4.14049 1 1 1500 0.439061-0.1

29、533511 1 3000 0.999996 4.447193 8.3期望期望-预期表：预期表：Expectation-Prediction Evaluation for Binary Specification Equation:UNTITLED Date:07/10/09 Time:23:02 Success cutoff:C=0.5 Estimated Equation Constant Probability Dep=0 Dep=1 TotalDep=0Dep=1Total P(Dep=1)C 2 13 15000 Total 15 15 30151530 Correct 13 13

30、2615015%Correct 86.67 86.67 86.67100.000.0050.00%Incorrect 13.33 13.33 13.330.00100.0050.00 Total Gain*-13.33 86.67 36.67 Percent Gain*NA 86.67 73.33 Estimated Equation Constant Probability Dep=0 Dep=1 TotalDep=0Dep=1Total E(#of Dep=0)13.03 1.97 15.007.507.5015.00 E(#of Dep=1)1.97 13.03 15.007.507.5

31、015.00 Total 15.00 15.00 30.0015.0015.0030.00 Correct 13.03 13.03 26.057.507.5015.00%Correct 86.84 86.84 86.8450.0050.0050.00%Incorrect 13.16 13.16 13.1650.0050.0050.00 Total Gain*36.84 36.84 36.84 Percent Gain*73.68 73.68 73.68 *Change in%Correct from default(constant probability)specification*Perc

32、ent of incorrect(default)prediction corrected by equation 8.4 H-L Statistic，Andrews Statistic 当自变量数量增加，尤其是连续自变量纳入模型后，协变类型的数量很大，所以很多协变类型只有很少的观测案例，皮尔逊当自变量数量增加，尤其是连续自变量纳入模型后，协变类型的数量很大，所以很多协变类型只有很少的观测案例，皮尔逊2不再适用于拟合优度检验。所以采用不再适用于拟合优度检验。所以采用H-L Statistic，Andrews Statistic，零假设为拟合完全充分，零假设为拟合完全充分 Goodness-of

33、-Fit Evaluation for Binary Specification Andrews and Hosmer-Lemeshow Tests Equation:UNTITLED Date:07/11/09 Time:00:30 Grouping based upon predicted risk(randomize ties)Quantile of Risk Dep=0 Dep=1 Total H-L Low High ActualExpectActualExpectObs Value 1 4.E-06 6.E-05 32.9999208.5E-053 8.5E-05 2 0.0002

34、 0.0018 32.9973600.002643 0.00264 3 0.0046 0.0231 32.9616200.038383 0.03888 4 0.0458 0.1415 32.7288900.271113 0.29805 5 0.2216 0.4391 12.0165920.983413 1.56336 6 0.5609 0.7784 20.9834112.016593 1.56336 7 0.8585 0.9542 00.2711132.728893 0.29805 8 0.9769 0.9954 00.0383832.961623 0.03888 9 0.9982 0.999

35、8 00.0026432.997363 0.00264 10 0.9999 1.0000 08.5E-0532.999923 8.5E-05 Total 1515.00001515.000030 3.80603 H-L Statistic 3.8060 Prob.Chi-Sq(8)0.8742 Andrews Statistic 18.7041 Prob.Chi-Sq(10)0.0442 8.5异方差检验编程异方差检验编程(1)原理：原理：0H:同方差同方差 1H:有如下形式的的异方差：有如下形式的的异方差：()exp(2)iiVar uz=z是影响随机误差项方差的自变量，是影响随机误差项方差

36、的自变量，为参数。为参数。构造构造LM统计量是辅助回归式统计量是辅助回归式12(*)(*)(*)(1)(1)(1)iiiiiiiiiiypfyfyyx bz bpppppp=+的拟合值的平方和，该的拟合值的平方和，该LM统计量近似服从自由度为统计量近似服从自由度为z中变量个数的中变量个数的2分布。分布。ip为二元选择模型得到的拟合概率值，为二元选择模型得到的拟合概率值，*y为二元选择模型得到的潜回归的拟合值。为二元选择模型得到的潜回归的拟合值。z是影响随机误差项方差的自变量，是影响随机误差项方差的自变量，12,b b为待估参数，为待估参数，x是其他的自变量。是其他的自变量。(2)步骤：步骤：（

37、1）生成生成(1)iiiiyppp，命名为命名为se（2）生成生成ip，命名为命名为yp，生成生成*y，命名为命名为yf（3）构造序列构造序列(*)(1)iifypp，命名为命名为fac，series fac=dnorm(-yf)/sqrt(yp*(1-yp)（4）生成截距项的序列生成截距项的序列 series x0=1（5）用最小二乘法拟合辅助回归式用最小二乘法拟合辅助回归式 ls se(x0*fac)(x*fac*(-yf)（6）生成辅助回归的拟合值生成辅助回归的拟合值sef（7）构造构造LM统计量统计量 scalar lm_test=sumsq(sef)（8）计算计算p值值 scalar

38、 p_val=1-cchisq(lm_test,1)结果为结果为P=0.999，所以不存在异方差，所以不存在异方差 九、重复观测值可以得到情况下二元离散选择模型的参数估计（九、重复观测值可以得到情况下二元离散选择模型的参数估计（stata）举例：研究年龄和接受某辐射对人群致死的影响举例：研究年龄和接受某辐射对人群致死的影响 age exposed deaths pop 0 0 5 120 0 1 8 106 1 0 16 85 1 1 22 98 GLS估计：估计：gprobit deaths pop exposed age ML估计：估计：bprobit deaths pop exposed

39、 age 十、多重选择模型十、多重选择模型 10.1 背景问题背景问题（1）城市旅行交通工具选择：火车、汽车、飞机；）城市旅行交通工具选择：火车、汽车、飞机；（2）多职业选择；）多职业选择；（3）等级评定的多重选择。等等。）等级评定的多重选择。等等。10.2 定义定义 称因变量取值超过两个的离散选择模型为多重选择模型；若多重选择有序，则称之为有序反应模型；若多重选择无序，则称之为无序反应模型。称因变量取值超过两个的离散选择模型为多重选择模型；若多重选择有序，则称之为有序反应模型；若多重选择无序，则称之为无序反应模型。如如 排序选择模型（有序）：银行信誉等级、学生奖学金等级等，数字反映好坏程度，

40、数字之间的差有意义。排序选择模型（有序）：银行信誉等级、学生奖学金等级等，数字反映好坏程度，数字之间的差有意义。一般选择模型（无序）：交通工具选择、超市购物选择等，数字之间的差值没有任何意义。一般选择模型（无序）：交通工具选择、超市购物选择等，数字之间的差值没有任何意义。十一、排序选择模型十一、排序选择模型 11.1 潜回归假设潜回归假设 假 设 存 在 一 个 尚 未 观 察 到 的 连 续 性 随 机 变 量假 设 存 在 一 个 尚 未 观 察 到 的 连 续 性 随 机 变 量*iy，*11iikikiiyxxxu=+=+?，即为潜回归。，即为潜回归。*iiyy与 关系是，如果连续性随

41、机变量关系是，如果连续性随机变量*iy超过某个临界值，则对应超过某个临界值，则对应iy的一个确定性选择。如买房多重选择为：租房、小户型、大户型、别墅，回归变量为收入，则个体选择依赖于的一个确定性选择。如买房多重选择为：租房、小户型、大户型、别墅，回归变量为收入，则个体选择依赖于*iy与临界值与临界值123,c c c的关系：的关系：112233*iiiiyccyccycyc租房小户型大户型别墅 假设多重选择有假设多重选择有1m+个，则个，则*iiyy与 关系具体如下：关系具体如下：*1*12*23*012iiiiimyccycycycmyc=?其中临界值关系满足：其中临界值关系满足：12mcc

42、c?。11.2 排序选择模型排序选择模型 相应地，若假设干扰项相应地，若假设干扰项iu服从标准正态分布或服从标准正态分布或logit分布，则可得排序选择模型的概率形式：分布，则可得排序选择模型的概率形式：*111(0)()()()iiiiiP yP ycPxucF cx=+=*121221(1)()()()()iiiiiiP yP cycP cxucF cxF cx=+=11.3 模型的参数估计模型的参数估计 假设随机抽样样本为：假设随机抽样样本为：(,)iiy x 1,2,in=?，利用样本确定模型参数：，利用样本确定模型参数：（1）回归系数）回归系数01,k?；（2）临界值）临界值12,m

43、c cc?参数估计的方法：极大似然估计。参数估计的方法：极大似然估计。11.4 EVIEWS应用举例（政府支持率的民意测验）应用举例（政府支持率的民意测验）态度态度iy是离散因变量，是离散因变量，0,1,2iy=，分别表示第，分别表示第i个市民的支持、中立、不支持。序体现在个市民的支持、中立、不支持。序体现在iy数值越大，支持度越低。回归变量数值越大，支持度越低。回归变量ix为市民收入，样本数据如下表：为市民收入，样本数据如下表：收入收入ix 态度态度iy 收入收入ix 态度态度iy收入收入ix 态度态度iy 550 0 850 1 1250 2 600 0 950 1 1350 2 650

44、0 1050 1 1450 2 700 0 1100 1 1500 2 750 0 1150 1 1550 2 800 0 1200 1 1600 2 900 0 1300 1 1650 2 1000 0 1400 1 1700 2 若假设干扰项若假设干扰项iu服从标准正态分布，请计算排序选择的服从标准正态分布，请计算排序选择的Probit的模型，并解释。的模型，并解释。待计算的多元选择模型：待计算的多元选择模型：*1*12*2012iiiiycycyccy=待计算的潜回归模型：待计算的潜回归模型：*iiiyxu=+其中其中(0,1)iuN 待估参数：待估参数：12,c c 利用利用Eview

45、s软件，计算结果如下：软件，计算结果如下：Dependent Variable:Y Method:ML-Ordered Probit(Quadratic hill climbing)Date:07/11/09 Time:22:33 Sample:1 24 Included observations:24 Number of ordered indicator values:3 Convergence achieved after 6 iterations Covariance matrix computed using second derivatives Variable Coefficie

46、ntStd.Errorz-StatisticProb.X 0.0087510.0028703.0495890.0023 Limit Points LIMIT_1:C(2)8.0930752.7039272.9930810.0028LIMIT_2:C(3)11.597793.8282363.0295400.0024Pseudo R-squared 0.674560 Akaike info criterion 0.965066Schwarz criterion 1.112322 Log likelihood-8.580789Hannan-Quinn criter.1.004133 Restr.lo

47、g likelihood-26.36669LR statistic 35.57181 Avg.log likelihood-0.357533Prob(LR statistic)0.000000 可得潜回归模型为：可得潜回归模型为：*0.008751iiiyxu=+注：因为对政府支持等级与序数反方向变化，即序数越大，支持度越低，反之亦然，故注：因为对政府支持等级与序数反方向变化，即序数越大，支持度越低，反之亦然，故*iy可理解为对政府的不满倾向。可理解为对政府的不满倾向。临界值为：临界值为：18.093074c=，211.59779c=可得支持态度可得支持态度iy的分布如下：的分布如下：*208

48、.09374,18.0937411.597792iiiiyyyyc=表示支持表示中立表示反对 给定收入给定收入ix，可得确定的态度，可得确定的态度iy的分布列：的分布列：iy的取值的取值 iy的取值的概率的取值的概率 0 1 2(8.0930.0087)ix (11.60.0087)(8.0930.0087)iixx 1(11.60.0087)ix 模型预测：模型预测：Procs/make model:i_y =0.00875149689847*x y_0 =cnorm(8.09307479625 -i_y)y_1 =cnorm(11.5977932253-i_y)-cnorm(8.09307

49、479625-i_y)y_2=1-cnorm(11.5977932253-i_y)收入收入ix 态 度态 度iy*iy 收入收入ix态 度态 度iy*iy 收入收入ix态 度态 度iy*iy 550 0 4.813323 850 1 7.4387721250 2 10.93937600 0 5.250898 950 1 8.3139221350 2 11.81452650 0 5.688473 1050 1 9.1890721450 2 12.68967700 0 6.126048 1100 1 9.6266471500 2 13.12725750 0 6.563623 1150 1 10.0

50、64221550 2 13.56482800 0 7.001198 1200 1 10.501801600 2 14.00240900 0 7.876347 1300 1 11.376951650 2 14.439971000 0 8.751497 1400 1 12.252101700 2 14.87754 期望期望-预期表预期表 Prediction Evaluation for Ordered Specification Equation:UNTITLED Date:07/11/09 Time:22:58 Estimated Equation Dep.Value Obs.Correct