计量经济学第十章定性选择模型ppt课件

1、第十章 定性选择模型 我们在第四章中曾引见解释变量为虚拟变量的模型，本章要讨论的是因变量为虚拟变量的情形。在这种模型中，因变量描画的是特征、选择或者种类等不能定量化的东西，如乘公交还是本人开车去上班、考不考研讨生等。在这些情况下，因变量是定性变量，我们可以用定义虚拟变量的方法来描写它们。这种因变量为虚拟变量的模型被称为定性选择模型Qualitative choice models或定性呼应模型Qualitative response models。 假设只需两个选择，我们可用0和1 分别表示它们，如乘公交为0，自驾车为1，这样的模型称为二元选择模型binary choice Models，多于

2、两个选择如上班方式加上一种骑自行车的定性选择模型称为多项选择模型Multinomial choice models。 第一节 线性概率模型 二元选择模型如何估计呢？由于它看上去象是一个典型的OLS回归模型，因此一个简单的想法是采用OLS法估计。当然，对结果的解释与常规线性回归模型不同，由于二元选择模型中因变量只能取两个预定的值。线性概率模型LPM普通方式如下： 这看上去与典型的OLS回归模型并无两样，但区别是这里Y只取0和1两个值，观测值可以是个人、公司、国家或任何其他横截面个体所作的决议。解释变量中可以包括正常变量和虚拟变量。01122(10.1)iiikkiiYXXXu 下面用一个关于能否

3、读研讨生的例子来阐明如何解释线性概率模型的结果。模型为：012(10.2)iiiiYGPAINCOMEu其中：01iY该生三年内未去读研三年内去读研个学生拿到学士学位后第iiiG P A 第 个 学 生 本 科 平 均 成 绩位：千美元）个学生家庭年收入（单第iiINCOME 设回归结果如下一切系数值均在10%程度统计上显著： 0.70.40.00210.3iiiYGPAINCOME （） 对每个观测值，我们可根据10.3式计算因变量的拟合值或预测值。在常规OLS回归中，因变量的拟合值或预测值的含义是，平均而言，我们可以预期的因变量的值。但在本例的情况下，这种解释就不适用了。假设学生甲的平均分

4、为3.5，家庭年收入为5万美圆，Y的拟合值为0.7 0.4 3.5 0.002 500.810.4Y （） 虽然因变量在这个二元选择模型中只能取两个值：0或1，可是该学生的的拟合值或预测值为0.8。我们将该拟合值解释为该生决议读研的概率的估计值。因此，该生决议读研的能够性或概率的估计值为0.8。需求留意的是，这种概率不是我们能观测到的数字，能观测的是读研还是不读研的决议。 对斜率系数的解释也不同了。在常规回归中，斜率系数代表的是其他解释变量不变的情况下，该解释变量的单位变动引起的因变量的变动。而在线性概率模型中，斜率系数表示其他解释变量不变的情况下，该解释变量的单位变动引起的因变量等于1的概率

5、的变动。 GPA的系数估计值0.4意味着家庭收入不变的情况下，一个学生的GPA添加一个点如从3.0到4.0，该生决议去读研的概率的估计值添加0.4。 INCOME的系数估计值0.002阐明，一个学生的成果不变，而家庭收入添加1000美圆，该生决议去读研的概率的估计值添加0.002。 LPM模型中，解释变量的变动与虚拟因变量值为1的概率线性相关，因此称为线性概率模型。线性概率模型存在的问题1线性概率模型假定自变量与Y=1的概率之间存在线性关系，而此关系往往不是线性的。2拟合值能够小于0或大于1，而概率值必需位于0和1的闭区间内。 回到有关读研的例子。假设学生乙的GPA为4.0，家庭收入为20万美

6、圆，那么代入10.3式，Y的拟合值为 从而得到一个不能够的结果概率值大于1。假设另有一个学生丙的GPA为1.0，家庭收入为5万元，那么其Y的拟合值为 -0.2，阐明读研的概率为负数，这也是一个不能够的结果。0.7 0.4 4.0 0.002 200 1.310.5Y （） 处理此问题的一种方法是，令一切负拟合值都等于0，一切大于1的拟合值都等于1。但也无法令人非常称心，由于在现实中很少会有决策前某人读研的概率就等于1的情况，同样，虽然某些人成果不是很好，但他去读研的时机仍会大于0。线性概率模型倾向于给出过多的极端结果：估计的概率等于0或1。 3 另一个问题是扰动项不是正态分布的。现实上，线性概

7、率模型的扰动项服从二项分布。 4此外，线性概率模型存在异方差性。扰动项的方差是 ，这里 是因变量等于1的概率，此概率对于每个观测值不同，因此扰动项方差将不是常数，导致异方差性。可以运用WLS法，但不是很有效，并且将改动结果的含义。(1)ppp 5最后一个问题是在线性概率模型中， 以及 不再是适宜的拟合优度测度。现实上，此问题不仅是线性概率模型的问题，而是一切定性选择模型的问题。较好一点的测度是模型正确预测的观测值的百分比。首先，我们将每一预测归类为1或0。假设拟合值大于等于0.5，那么以为因变量的预测值为1。假设小于0.5，那么以为因变量的预测值为0。然后，将这些预测值与实践发生的情况相比较，

8、计算出正确预测的百分比： 2R2R100观测值总数正确预测的观测值数比正确预测观测值的百分 需求指出的是，这个测度也不是很理想，但预测结果的好坏，并非定性选择模型独一关怀的事，这类模型常被用于研讨影响人们进展某个决策的要素。让我们来看一个竞选的例子。假设候选人甲和乙二人竞选某市市长，我们可以用一个二元选择模型来研讨影响选民决策的要素，数据见表101，模型为：0123110.6iiiiiCANDINCOMEAGEMALEu（）其中：110iCAND票个选民不投候选人甲的如果第个选民投候选人甲的票如果第ii位：千美元）个选民的家庭收入（单第iiINCOME个选民的年龄第iiAGE01iMALE女性

9、男性VariableCoefficientStandard errort-Statisticp-ValueConstant-0.510.19-2.650.01 INCOME0.00980.0033.250.00 AGE0.0160.00533.080.00 MALE0.00310.130.020.98 表10-2 两候选人选举线性概率模型回归结果Dependent variable：CAND1Observations：30 = 0.58Adjusted = 0.53Residual Sum of Squares =3.15F-statistic = 11.872R2R 如表102所示，INCO

10、ME的斜率估计值为正，且在1%的程度上显著。年龄和性别不变的情况下，收入添加1000元，选择候选人甲的概率添加0.0098。 AGE的斜率估计值也在1%的程度上显著。在收入和性别不变的情况下，年龄添加1岁，选择候选人甲的概率添加0.016。MALE的斜率系数统计上不显著，因此没有证听阐明样本中男人和女人的选票不同。 我们可以得出如下结论：年老一些、富有一些的选民更喜欢投票给候选人甲。 表103给出CAND1的拟合值，每个大于等于0.5的拟合值计入CAND1为1的预测，而小于0.5的拟合值那么计入CAND1为0的预测。 从表103可看出，30个观测值中，27个或90%预测正确。选甲的14人中，1

11、2人或85.7%预测正确。选乙的16人中，15人或93.8%预测正确。 是0.58，阐明模型解释了因变量的58%的变动，这与90%的正确预测比例相比，低了不少。留意表103中有一些拟合值大于1或小于0。这是我们前面指出的这类模型的缺陷之一，这些拟合值是概率的估计值，而概率永远不能够大于1或小于0。2R第二节 Probit模型和Logit模型一Probit和Logit方法概要 估计二元选择模型的另一类方法假定回归模型为 这里 不可观测，通常称为潜变量latent variable。我们能观测到的是虚拟变量：*01(10.7)kijijijYXu01iY*iY010.8若（）其它*iY 这就是Pr

12、obit和Logit方法的思绪。Probit模型和Logit模型的区别在于对10.7式中扰动项u的分布的设定，前者设定为正态分布，后者设定为logistic分布。 10.7式与线性概率模型的区别是，这里假设潜变量的存在。例如，假设被观测的虚拟变量是某人买车还是不买车， 将被定义为“买车的愿望或才干，留意这里的提法是“愿望和“才干，因此10.7式中的解释变量是解释这些元素的。 从10.8式可看出， 乘上任何正数都不会改动 ，因此这里习惯上假设 Var(ui) = 1，从而固定 的规模。由10.7和10.8式，我们有*iY*iY*iYiY0101Pr(1)Pr()1 ()kiiijijjkjijj

13、Pob Yob uXFX 其中F是u的累积分布函数。假设u的分布是对称的，那么 ，我们可以将上式写成1()( )FzF z01()(1 0 .9 )kijijjPFX我们可写出似然函数：10(1)(10.10)iiiiYYLPP 10.9式中F的函数方式取决于有关扰动项u的假设，假设 的累积分布是logistic分布，那么我们得到的是logit模型。在这种情况下，累积分布函数为：iuexp()()(10.11)1exp()iiizF zz因此()log1()iiiF zzF z这是由于，由10.11式，有：()log1()exp()exp()1exp()1exp()loglogexp()1ex

14、p()exp()11exp()1exp()exp()1exp()loglog exp()11exp()iiiiiiiiiiiiiiiiF zF zzzzzzzzzzzzzzz结合10.9式，对于logit模型，有：01log1kiijijjipXp 上式的左端是时机odds的对数，称为对数时机比率log-odds ratio，因此上式阐明对数时机比率是各解释变量的线性函数，而对于线性概率模型， 为各解释变量的线性函数。 假设10.9式中 服从正态分布，我们得到的是probit模型或normit模型，在这种情况下，累积分布函数为：ip/21( )exp()(10.12)22izitF zdtiu

15、 无论是probit模型还是logit模型，极大似然函数10.10都伴随着非线性估计方法，目前很多计量经济分析软件已可用于probit和logit分析，用起来很方便。 由于累积正态分布和累积logistic分布很接近，只是尾部有点区别，因此，我们无论用10.11还是10.12，也就是无论用logit法还是probit法，得到的结果都不会有很大不同。可是，两种方法得到的参数估计值不是直接可比的。由于logistic分布的方差为 ，因此，logit模型得到的的估计值必需乘以 ，才干与probit模型得到的估计值相比较正态分布规范差为1。233二Probit模型 Probit模型可以处理很多线性概率

16、模型中遇到的问题。如我们在前面指出的，线性概率模型会给出小于0或大于1 的这种不能够的概率估计值，Probit模型所根据的是累积正态概率分布，将防止这类问题的发生，同时它给出接近0或1的概率估计值的时机也要小于线性概率模型。与线性概率模型相比，Probit模型更准确地描画我们计划研讨的许多决策过程。如图10-1所示概率=F(Z)10ZProbit模型线性概率模型图10-1 线性概率模型和Probit模型 虽然Probit模型实践是非线性的，但它可以以一种类似于其他经济模型的方式写出。首先，我们需求将等式10.12略微改写一下，它代表由累积正态概率函数执行的变换：/21( )exp()(10.1

17、3)22iziitPF zdt在上式中，F是一个函数，即将正态概率函数的一个值转换成概率的累积正态概率函数。Probit模型运用其反函数，将概率值转换成 Z 的值。1(10.14)iZFPProbit模型为 101122(10.15)iiiikkiiZFPXXXu虽然乍看上去上式像一个典型的回归模型，但它是一个非线性模型，由于有 这一项。Probit模型不能用OLS法估计，应采用极大似然法估计。 1iFP Probit模型以及我们下面要讨论的Logit模型在大样本观测值数以百计时运转最好。假设样本中两种能够的选择都有足够的信息，那么效果更佳。例如，对于我们前面的读研讨生的例子，设观测值为200

18、，假设其中仅3%的人决议读研，也就是200人中仅有6人，那么在此样本中就没有足够的信息来给出好的估计值，选择读研的样本过少，使得回归结果的可信程度不高。 我们可以将两个候选人的选举模型用Probit模型估计，运用与前面一样的变量和数据，估计结果如表10-4所示。表104 两候选人选举模型的Probit回归结果Dependent variable：CAND1VariableCoefficientStandard errort-Statisticp-ValueConstant-5.191.70 -3.060.00 INCOME0.0710.0342.10 0.04 AGE0.0730.0342.1

19、80.03 MALE-0.70 0.90 -0.780.44 Observations：30McFadden pseudo-R2 = 0.61Residual Sum of Squares = 2.62 采用Probit模型估计的结果与前面用线性概率模型估计的结果有所不同。采用Probit模型的情况下，INCOME和AGE的系数估计值在5%的误差程度上显著，而在线性概率模型的情况下，在1%的程度上显著。 由于我们知道线性概率模型存在严重的问题，因此Probit结果能够更准确一些。可是，假设是实践研讨的话，要有一个大得多的样本。Probit模型的系数估计值不能像线性概率模型那样，解释成概率的变动

20、。运用Probit模型的一种有意思的方式是求出拟合值进展预测，如我们用线性概率模型所做的一样表10-3。 Probit模型中用McFadden的pseudo-R2作为拟合优度的测度。pseudo-R2是用于虚拟因变量模型的拟合优度的测度的名字。pseudo-原意是伪假，这里采用它，意思是与常规R2类似但不一样，而不是说它是假的。 对于定性选择模型，曾经开发了几种有用的pseudo-R2测度，这里所用的是McFadden开发的。很多估计Probit或Logit模型的计量经济程序计算pseudo-R2。本例中给出的0.61的含义是，Probit模型解释了因变量61%的变动。三. Logit模型 L

21、ogit模型基于累积logistic分布，而不是probit模型所用的累积正态分布。对于任何一个回归，probit和logit估计方法的结果往往从统计显著性的角度看是类似的。 Logit模型给出的概率估计值限制在0和1之间，与probit一样，而且logit模型也防止了接近 0或1的极端概率值。这两个模型都抑制了线性概率模型遇到的主要问题。 Logit模型的方式如下：01122ln()(10 16)1iiikkiiiYXXXuY 在这里，因变量的拟合值代表 的能够性的对数。术语概率probability和时机odds不是一回事。假设一个事件的概率是0.25，那么时机将是：0.251110.25

22、3概率概率 我们通常将其写为1:3，读作1对3。假设概率是0.5或50%，那么相应为0.5/1-0.5=1/1，或1:1。我们可以给logit模型中斜率系数一个特别的解释：某个解释变量的变动对Y等于1的时机的影响。准确地说，logit模型的斜率系数通知我们，在其它解释变量坚持不变的情况下，该解释变量变动一个单位所引起的时机的对数的变动。 1iY 与probit模型一样，logit模型也不能用OLS法估计，而要用极大似然法估计。采用表10-1中的同样数据估计logit模型，回归结果如表10-5所示。表10-5 两候选人选举模型的Logit回归结果Dependent variable：CAND1V

23、ariableCoefficientStandard errort-Statisticp-ValueConstant-8.963.23-2.770.01INCOME0.120.061.980.05AGE0.130.062.030.04MALE-1.031.54-0.670.51Observations：30McFadden pseudo-R2 = 0.60Residual Sum of Squares = 2.59 McFadden pseudo-R2和统计显著性与probit模型的结果类似。INCOME和AGE的系数估计值亦在5%误差程度上显著。而MALE那么在两种模型回归中均不显著。而斜率

24、系数估计值那么不同，这是由于它们的意义不一样。例如，AGE的系数估计值0.13意味着收入和性别不变的情况下，年龄增大一岁，选举候选人甲的时机的对数添加0.13。实践上，除了斜率系数的解释不同，运用probit模型和logit模型并没有多大区别。 第三节 多项选择模型 我们能够遇到多于两个能够的选择的情况，如在选举模型例子中，有能够不止两个候选人，我们前面讨论的估计方法无法处置多于两项选择的情况。假设第三个候选人丙加进来了，我们就必需调整以前的估计方法，来思索加上第三项选择的情况。012345671(10.17)3(10.18)iiiiiiiiiiCANDINCOMEAGEMALEuCANDIN

25、COMEAGEMALEv其中， 11,031,0iiCANDiCANDi第 个选民选甲； ，其它第 个选民选丙； ，其它 两式的系数下标不一样，阐明两方程的系数可以取不同的值。我们用OLS法估计这两个方程，存在的问题与两个选择的情况一样。一. 线性概率模型 线性概率模型经过修正，可用于多于两项选择的非定序的情况。要将第三个候选人加到我们的选举模型，我们需求用两个方程普通而言，方程的数目是选择数目减1。 对于任何一个观测值，估计出的概率之和必需等于1。第i个选民选甲的概率的估计值由10.17式中因变量CAND1的拟合值给出，比如说0.5，与此类似，该选民选丙的概率的估计值由10.18式中因变量C

26、AND3的拟合值给出，如0.3，那么我们知道，该选民选乙的概率估计值为0.2，这三个估计的概率之和必需等于1。 因此，我们无需为候选人乙回归第三个方程。现实上，三个候选人截距的估计值之和等于1，各斜率的估计值之和为0，因此我们估计两个方程后，第三个方程的斜率就可以算出来了。 对线性概率模型进展的这种修正只适用于各个方程中的解释变量都一样的情况。否那么，就必需用较复杂的GLS法。 表10-1中没有包括支持第三个候选人丙的选民的有关数据，表106列出了这些数据。这最后10个观测值都支持候选人丙并非巧合，它们未必是原样本中最后10个观测值，只不过是表10-1中省略了一切支持丙的观测值。将这些数据加到

27、表10-1的数据中，我们就得到一个包含三种选择的数据集，观测值数目为40。 要留意的是，在将表106的数据加到原来的30个观测值中的同时，CAND3变量代表候选人丙也应该加到原来的30个观测值中，CAND3在前30个观测值中取值为0。用这个新的数据集估计10.17、10.18式，估计结果如表10-7和表10-8所示。表10-6 选举模型增补观测值：支持候选人丙的个体观测值观测序号CAND1INCOMEAGEMALECAND331022191132024201133030221134021241135026211136030340137029241138033251139028271140032

28、3011表10-7 三候选人选举线性概率模型回归结果Dependent variable：CAND1VariableCoefficientStandard errort-Statisticp-ValueConstant-0.580.16-3.710.00 INCOME0.010 0.00273.740.00 AGE0.0170.00434.050.00 MALE-0.0350.099-0.350.73 Observations：40 = 0.62Adjusted = 0.59Residual Sum of Squares =3.41F-statistic = 19.992R2R表10-8 三候

29、选人选举线性概率模型回归结果Dependent variable：CAND3VariableCoefficientStandard errort-Statisticp-ValueConstant0.480.192.50 0.02 INCOME-0.000850.0033-0.26 0.80 AGE-0.0110.0053-2.06 0.05 MALE0.330.122.69 0.01 Observations：40 = 0.30Adjusted = 0.25Residual Sum of Squares =5.19F-statistic = 5.352R2R 表10-7阐明候选人甲作为因变量的

30、方程的结果与二元选择线性概率模型的结果类似与表10-2比较，对斜率系数的阐明也可沿用二元选择模型同样的方式。例如，AGE的斜率系数0.017意味着，INCOME和MALE坚持不变的情况下，选民的年龄大一岁，选甲的概率上升0.017。 表10-8中候选人丙的结果那么与甲的结果大不一样。INCOME的斜率估计值在甲的方程中显著，但在丙的方程中那么不显著。高收入者倾向于选甲，低收入者倾向于选丙或选乙，但收入似乎不怎样影响对丙的选择。AGE的斜率估计值在5%误差程度显著，其值为负，阐明年轻选民倾向于选丙，与候选人甲的情况刚好相反，甲的方程阐明，年龄较大的选民倾向于选甲。具有同样年龄和收入的男选民选丙的

31、估计概率比女选民高0.33，这是一个很大的差距。在三个候选人中，妇女最不接受的人是丙。二. 多项logit模型方法多项Logit模型Mutinomial logit用于估计多于两项选择的定性选择模型这些选择没有先后次序，该方法防止了线性概率模型出现的问题。与线性概率模型一样，所需求的方程的个数是选择的数目减1，其中一个选择被用作基准选择，该选择没有本人的方程。将多项logit模型运用于三候选人的选举模型，我们用候选人乙作为基准选择，给出下面两个方程：01234567CAND1ln()INCOME(13)CAND2CAND3ln()INCOME(14)CAND2iiiiiiiiiiiiAGEMALEuAGEMALEv其中CAND21,0ii若第 个选民选乙， ，其它。 多项logit模型中的方程必需用极大似然法联立地估计。大多数计量经济软件的学生版甚至全版不支持这类估计。 采用多项logit模型估计出的斜率系数的解释与二元logit模型不一样。在这里，