管理过程的数学模型

管理过程的数学模型主要内容虚拟解释变量虚拟因变量滞后变量趋势变量面板数据非线性模型第一页，共八十四页。3.1虚拟解释变量一.虚拟解释变量模型虚拟变量：某些因素对解释因变量是必须的，但它们是定性的、不可计量的；为了将这些变量引入所要研究的模型，必须将它们数量化：比如起作用时赋值1，或0；不起作用时赋值0，或1，这样的变量称为虚拟变量比如：天气、季节、性别等本节讨论虚拟变量作为解释变量的模型第二页，共八十四页。例1研究学生体重与身高的关系，随机抽样80名学生，其中48名男生，32名女生（W表示体重，单位磅；h表示身高，单位英寸）模型A：(-5.2066)(8.6246)模型B:(-2.5884)(4.0149)(5.1613)

其中第三页，共八十四页。虚拟变量的影响hW男生女生截距项位移第四页，共八十四页。例2研究工作经验对工资的影响，假设由于法律的原因，男性和女性参加工作时的起薪是一样的建立模型（Y是月薪，X是工龄）第五页，共八十四页。虚拟变量的影响斜率项位移hW男性女性第六页，共八十四页。例3对例2，如果男女参加工作时的起薪存在性别歧视，应如何修正模型？第七页，共八十四页。虚拟变量的影响截距项和斜率项均位移hW男性女性第八页，共八十四页。例4-A研究学历对工作的影响，建立模型（Y:参加工作的起薪）模型A:第九页，共八十四页。例4-B模型B:第十页，共八十四页。虚拟变量陷阱对模型B正规方程系数矩阵非满秩矩阵存在多重共线性第十一页，共八十四页。虚拟变量陷阱欲表征的状态数等于虚拟变量个数此时存在完全多重共线性因此，虚拟变量个数=欲表征的状态数-1第十二页，共八十四页。例4-C模型C：第十三页，共八十四页。季节性变动虚拟变量销售函数模型考虑销售量的季节性波动，应如何引入虚拟变量?i=1,2,3第十四页，共八十四页。二.结构稳定性检验与分段线性回归1.结构稳定性检验检验模型反映的经济结构是否由于受到某种因素的影响而发生变化例如:时间序列数据，1992年后解释变量的参数可能发生变化界面数据，东部省份和中西部省份的回归结果不一致第十五页，共八十四页。结构稳定性的Chow检验对总样本进行回归将总样本分成两个子样本，分别回归构造统计量第十六页，共八十四页。参数稳定性Chow检验判断样本扩大时，模型参数的估计是否具有稳定性对原样本（n1）进行回归增加n2个观测值，组成新样本（n1+n2）进行回归构成统计量第十七页，共八十四页。举例研究各省市旅游外汇收入，建立模型（Y：旅游外汇收入；X1旅行社职工人数；X2国际旅游人数）：得到回归方程：(3.067)(6.653)(3.378)分别对东部地区（东北、华北、华东、华南）15省市和西部地区（华中、西南、西北）16省市进行回归第十八页，共八十四页。2.分段线性回归模型结构发生变化，但回归函数保持连续在一个时刻结构发生变化：XYb1第十九页，共八十四页。第二十页，共八十四页。在两个时刻结构发生变化：第二十一页，共八十四页。3.2虚拟因变量

(受限因变量、分类选择模型)被解释变量是定性变量。如家庭是否拥有自己的住宅，企业是否在某个地区投资，成年男子是否在“参与劳动”等。常见模型：线性概率模型（LPM）对数单位模型（LogitModel）概率单位模型（ProbitModel）第二十二页，共八十四页。一.线性概率模型与一般线性模型相比，Y不再服从正态分布，而是一个二项分布第二十三页，共八十四页。第二十四页，共八十四页。线性概率模型的问题(1)随机项非正态性第二十五页，共八十四页。(2)异方差概率总和1第二十六页，共八十四页。异方差的处理随机项异方差，OLS估计量线性无偏但不是有效的用加权最小二乘法（WLS）处理方程两边同时除以第二十七页，共八十四页。(3)值域问题问题：条件期望的值可能超出[0,1]区间，这是线性概率模型的严重缺点处理：大于1的当作1，小于0的当作0采用Logit或Probit模型第二十八页，共八十四页。(4)拟合优度问题线性概率模型的拟合优度一般不高，在0.2到0.6之间LPMLPM对于受约束的LPM(b)一般不会大，大多数实例第二十九页，共八十四页。二.Logit模型思路：出于线性概率模型的缺陷，希望对它进行变换，使预测对于所有的X都落在（0,1）之间第三十页，共八十四页。基本思路随着X增加，因变量增加（或减少），因此用一个累计概率函数F来描述可有多种累计概率函数，得出不同的模型，一般只采用两种：累计logistic概率函数——Logit模型累计正态概率函数——Probit模型第三十一页，共八十四页。第三十二页，共八十四页。第三十三页，共八十四页。第三十四页，共八十四页。第三十五页，共八十四页。

(收入等于的家庭个数）

（其中拥有住房的家庭数）640885012106018………402520第三十六页，共八十四页。处理异方差第三十七页，共八十四页。三.Probit模型第三十八页，共八十四页。第三十九页，共八十四页。LPMLogitProbit比较Logit模型最常用第四十页，共八十四页。3.3滞后变量模型一.滞后变量模型对采用时间序列数据的模型，相关变量的滞后值作为解释变量外生滞后变量模型内生滞后变量模型滞后变量模型实质上考虑的是经济系统的动态变化，也可称为动态模型第四十一页，共八十四页。1.外生滞后变量模型（分布滞后模型）举例：利率的动态模型r：利率M：货币供给D：国库券弥补的预算赤字第四十二页，共八十四页。(1)考伊克滞后（几何滞后）间隔时间越长，对因变量的影响越小，类似的模型都可用考伊克变换第四十三页，共八十四页。考伊克滞后的权重123465i(滞后)第四十四页，共八十四页。(2)阿尔蒙滞后（多项式滞后）一般取3次多项式，4期滞后第四十五页，共八十四页。根据这个模型可以求出各个参数第四十六页，共八十四页。考伊克变换只能表示影响递减的情形阿尔蒙模型的好处在于，它可以用多项式近似获得连续函数，反映各种复杂的影响，比如在两三个季度或两三年后才产生的影响123465i(滞后)第四十七页，共八十四页。2.内生滞后变量模型（自回归模型）(1)部分调整模型例：某家公司的理想库存水平为，实际库存水平，假设理想的库存水平由销售量决定：由于市场摩擦，实际水平和理想水平之间的差距不能迅速合拢，只能部分弥补第四十八页，共八十四页。(2)适应性预期模型例：假设居民的消费决定于期望收入根据早前实现的期望值对期望值进行修改期望方程(1)(2)(3)(3)代入(1)(1)滞后一期，并乘1-r(4)(5)第四十九页，共八十四页。(4)-(5)得：与前面考伊克变换得出的模型形式一样，不过参数的意义有所不同第五十页，共八十四页。二.因果关系检验检验两个变量是否存在因果关系，由Granger和Sims提出，称为Granger因果关系检验基本思想：如果X的变化（因）引起Y的变化（果），则X的变化应当发生在Y的变化之前X是引起Y变化的原因，则必须同时满足：X有助于预测YY有助于预测X第五十一页，共八十四页。对两变量Y与X，格兰杰因果关系检验要求估计:(*)(**)可能存在有四种检验结果：（1）X对Y有单向影响，表现为（*）式X各滞后项前的参数整体为零，而Y各滞后项前的参数整体不为零；（2）Y对X有单向影响，表现为（**）式Y各滞后项前的参数整体为零，而X各滞后项前的参数整体不为零；（3）Y与X间存在双向影响，表现为Y与X各滞后项前的参数整体不为零；（4）Y与X间不存在影响，表现为Y与X各滞后项前的参数整体为零。第五十二页，共八十四页。无限制条件回归有限制条件回归得残差平方和ESSUR得残差平方和ESSRn为观测个数k为无限制条件回归待估参数个数如果:F>F(m,n-k)，则拒绝原假设，认为X是Y的格兰杰原因。第五十三页，共八十四页。常见问题格兰杰因果关系检验对于滞后期长度的选择有时很敏感。不同的滞后期可能会得到完全不同的检验结果。因此，一般而言，常进行不同滞后期长度的检验，以检验模型中随机误差项不存在序列相关的滞后期长度来选取滞后期。第五十四页，共八十四页。举例：检验1978~2000年间中国当年价GDP与居民消费CONS的因果关系

中国GDP与消费支出（亿元）

年份

人均居民消费

CONSP

人均GDP

GDPP

年份

人均居民消费

CONSP

人均GDP

GDPP

1978

1759.1

3605.6

1990

9113.2

18319.5

1979

2005.4

4074.0

1991

10315.9

21280.4

1980

2317.1

4551.3

1992

12459.8

25863.7

1981

2604.1

4901.4

1993

15682.4

34500.7

1982

2867.9

5489.2

1994

20809.8

46690.7

1983

3182.5

6076.3

1995

26944.5

58510.5

1984

3674.5

7164.4

1996

32152.3

68330.4

1985

4589

8792.1

1997

34854.6

74894.2

1986

5175

10132.8

1998

36921.1

79003.3

1987

5961.2

11784.7

1999

39334.4

82673.1

1988

7633.1

14704.0

2000

42911.9

89112.5

1989

8523.5

16466.0

第五十五页，共八十四页。取两阶滞后，Eviews给出的估计结果为：

判断：=5%，临界值F0.05(2,17)=3.59拒绝“GDP不是CONS的格兰杰原因”的假设，不拒绝“CONS不是GDP的格兰杰原因”的假设。因此，从2阶滞后的情况看，GDP的增长是居民消费增长的原因，而不是相反。

但在2阶滞后时，检验的模型存在1阶自相关性。第五十六页，共八十四页。第五十七页，共八十四页。随着滞后阶数的增加，拒绝“GDP是居民消费CONS的原因”的概率变大，而拒绝“居民消费CONS是GDP的原因”的概率变小。如果同时考虑检验模型的序列相关性以及赤池信息准则，发现：滞后4阶或5阶的检验模型不具有1阶自相关性，而且也拥有较小的AIC值判断结果是：GDP与CONS有双向的格兰杰因果关系，即相互影响。

分析：第五十八页，共八十四页。3.4趋势变量(时间变量)用时间序列的观测时期所代表的时间作为模型的解释变量，用来解释被解释变量随时间推移的自发变化趋势。这种变量称为时间变量或趋势变量。只采用时间变量作为解释变量的模型称为增长曲线模型引入时间变量只是为了更好地解释被解释变量，但时间并不是经济活动变化的原因，因此本届讨论的都是非因果关系模型t的取值：t=1980，1981，1982，……t=1，2，3，……t=……-3，-2，-1，0，1，2，3，……(n为奇数)t=……，-5，-3，-1，1，3，5，……(n为偶数)第五十九页，共八十四页。常见增长曲线模型多项式增长曲线模型简单指数型增长曲线模型修正指数型增长曲线模型逻辑增长曲线模型龚珀兹增长曲线模型第六十页，共八十四页。逻辑(Logistic)增长曲线由Verlulst于1845年提出，用于模拟人口增长，俗称“S曲线”一般形式：常见的简化形式（狭义逻辑增长曲线模型）：第六十一页，共八十四页。逻辑(Logistic)增长曲线K第六十二页，共八十四页。龚珀兹(Gompertz)曲线由B.Gompertz于1825年提出上限逼近值K，下限逼近值0，与逻辑增长曲线相似，但二者的拐点位置不同第六十三页，共八十四页。含趋势变量的一般模型一个例子：研究我国资本外逃对经济增长的影响，建立如下模型（贺力平等，2004）：CF是资本外逃量，时间变量t用来概括影响实际GDP的一般因素第六十四页，共八十四页。3.5面板数据PanelData，翻译为面板数据或平行数据，是指包含若干个体在一个时间区域内(若干时点)的样本。样本中的每一个个体都具有很多观测(构成时间序列)在每个确定的时点也具有由各个个体数据组成的观测(构成截面数据)其中N是截面的个体数量，T是时间序列的时段个数第六十五页，共八十四页。优点和问题平行数据很有用，它可以使研究人员得到单用截面数据或单用时间序列数据都无法获得的经济信息。其它的好处还有：平行数据通常含有很多的数据点，样本具有较大的自由度；截面变量和时间变量的结合信息能够显著地减少缺省变量所带来的问题。另一方面，平行数据的使用也使模型的确认变得更加困难。平行数据的干扰可能包含时间序列干扰、截面数据干扰，以及时间序列与截面的混合干扰。第六十六页，共八十四页。估计方法融合方法：将所有的时间序列和截面数据互相融合(或者说混合在一起)，然后用LS估计可能的模型。固定效应模型：添加虚拟变量以便允许截距变化，这主要基于缺省变量可能引起截面截距和时间序列截距的变化。随机效应模型：考虑截面和时间序列的干扰(误差)改进第一种方法中LS估计的有效性。第六十七页，共八十四页。一.融合方法(普通最小二乘法)对前述面板数据模型，如果误差项满足古典线性模型假设，我们可以对截面数据逐个回归，比如对于t=1：共有T个这样的模型。类似地，我们还可以对时间序列数据逐个回归，比如对于i=1:共有N个这样的模型。如果α,β的真值对于时间序列和截面个体来说都是一样的常数，我们就可以混合所有数据，用NT个观测进行一个大的融合回归：第六十八页，共八十四页。二.固定效应模型最小二乘融合方法的问题在于常数截距和常数斜率的假设可能不合理。如果每个截面都是不同的模型，那么融合就不合适了。处理截距问题的最好办法，是引进允许截距项随时间和截面个体变化的虚拟变量，这就是固定效应模型：

其中第六十九页，共八十四页。是否添加虚拟变量可以通过统计的假设检验决定。检验就是比较两种平行数据运用方法的误差平方和。因为第一种方法比固定效应模型包含更多的参数限制条件(不同时间和不同个体的截距相等)，通常拟合程度更差，误差平方和会大些。如果添加的限制条件引起的误差平方和增加的不显著，就认为添加的限制条件合适，采用第一种方法(可以融合回归)；如果误差平方和的变化过大，我们就选择固定效应模型。检验统计量为第七十页，共八十四页。固定效应模型的优点和问题能够分析对任一给定的截面单位，因变量与整个截面均值之间的差异程度虚拟变量的使用不直接确认回归直线随时间和个体而变动的原因虚拟变量的引入会大大减少模型的自由度第七十一页，共八十四页。三.随机效应模型在第二种方法中采用虚拟变量，是考虑到第一种方法对信息的利用可能不够充分。另一种改善的方法是通过误差项来描述这种信息的不完整性，这就是随机效应模型其中它们彼此不相关，且不存在自相关。第七十二页，共八十四页。四.似不相关回归(SUR)一个更常用的方法是“似不相关回归”，由泽尔尼1962年提出，并可在一般的计量经济学软件（如Eviews）轻松实现假设误差项在各个横截面组中不同，但并不随时间而不同，例如有三个横截面误差项之间存在同期相关性第七十三页，共八十四页。似不相关回归的步骤Step1:使用OLS估计每个方程并求残差Step2:使用残差估计方差和协方差Step3:使用Step2中的估计值求所有参数的可行广义最小二乘估计值（FG