面板数据的统计分析方法-冯国双课件

面板数据的统计分析方法冯国双面板数据的统计分析方法冯国双主要内容1、面板数据简介2、面板数据的固定效应模型3、面板数据的随机效应模型4、动态面板数据模型5、空间面板数据模型主要内容1、面板数据简介面板数据横断面数据时间序列数据城市手足口报告发病率北京0.78上海6.82重庆2.50哈尔滨0.06郑州5.62广州14.22乌鲁木齐0.24月份手足口报告发病率1月0.782月0.763月6.144月24.715月53.776月67.737月47.268月19.379月15.2810月9.8811月9.1112月3.96面板数据横断面数据面板数据城市1月2月3月4月5月6月7月8月9月10月11月12月北京0.780.766.1424.7153.7767.7347.2619.3715.289.889.113.96上海6.823.1315.4234.3740.9741.3929.948.9210.487.117.197.20重庆2.501.254.5514.2414.2110.505.911.242.063.865.014.52哈尔滨0.060.010.221.133.746.4810.123.701.630.980.520.39郑州5.623.7018.5515.0123.1718.1011.416.963.893.002.321.76广州14.225.1918.0551.1664.4950.3433.6117.1224.3225.578.766.32乌鲁木齐0.240.040.242.0513.5827.1622.519.016.187.325.632.44面板数据（paneldata）面板数据城市1月2月3月4月5月6月7月8月9月10月11月面板数据面板数据：同时含有横断面和时间序列的数据对一组固定个体的多次观测数据其它称谓：重复追踪数据纵贯调查数据…………面板数据面板数据：面板数据微观面板数据：主要指个体的重复调查数据，一般个体数N较大，观察时间较短宏观面板数据通常涉及一段时间内国家、省等单位的重复调查数据，观察时间较长面板数据微观面板数据：面板数据模型面板数据用双下标变量yit表示其中，i=1,2,…,N

表示第i个个体，共N个个体t=1,2,…,T表示第t个观测时间，共观测T次yit表示第i个个体在第t时间点的观测值面板数据模型面板数据用双下标变量yit表示面板数据模型面板数据模型的一般形式αi是截距项i=1,2,…,N（N为个体数）t=1,2,…,T（T为每一个体的观察时点数）K为解释变量个数uit是随机误差项面板数据模型面板数据模型的一般形式混合模型（一）混合模型（pooledmodel）α表示对所有个体都有一个相同的截距项。混合模型（一）混合模型（pooledmodel）混合模型城市城市代码月份手足口发病率气温北京110.78-4.8北京120.76-1北京136.144.1北京1424.7111.2北京1553.7721.7北京1667.7324.7上海216.825.3上海223.137.5上海2315.429.3上海2434.3712.8上海2540.9720.8上海2641.3923.8哈尔滨310.06-17.1哈尔滨320.01-15.5哈尔滨330.22-7.2哈尔滨341.134哈尔滨353.7416哈尔滨366.4825.5乌鲁木齐410.24-10.6乌鲁木齐420.04-12.4乌鲁木齐430.24-4.4乌鲁木齐442.058.7乌鲁木齐4513.5816.2乌鲁木齐4627.1622.5混合模型城市城市代码月份手足口发病率气温北京110.78-4混合模型混合模型混合模型混合模型混合模型混合模型的参数估计：混合模型将所有数据混合在一起作为样本数据，可采用普通最小二乘法（ordinaryleastsquare，OLS）进行参数估计。混合模型混合模型的参数估计：混合模型混合模型的缺点：假定所有个体（如不同省、市、区县、家庭）都是同质的，没有区分不同个体，也无法说明是否不同个体之间存在差异。掩盖了个体间可能存在的异质性，把个体特征全部放到误差项中。混合模型混合模型的缺点：混合模型针对混合模型缺点的解决方法：需要体现出个体的差异，主要基于两种方法：1、个体固定效应模型2、个体随机效应模型混合模型针对混合模型缺点的解决方法：固定效应模型（二）个体固定效应模型（fixedeffectmodel）把各个个体的值设定为一套固定的参数，采用一定的方法估计出来个体固定效应模型根据不同情况还可分为：1、变截距的个体固定效应模型2、变斜率/系数的个体固定效应模型固定效应模型（二）个体固定效应模型（fixedeffect固定效应模型1、变截距的个体固定效应模型：αi反映了个体对总均值的偏离，表示不同个体有不同的截距项。固定效应模型1、变截距的个体固定效应模型：固定效应模型变截距的固定效应模型固定效应模型变截距的固定效应模型固定效应模型变截距固定效应模型的参数估计：最小二乘虚拟变量（leastsquaredummyvariable，LSDV）利用虚拟变量的方式，估计出所有个体的截距值如个体固定效应模型可以表示为：共需估计N个虚拟变量系数和K个自变量系数固定效应模型变截距固定效应模型的参数估计：固定效应模型2、变系数的个体固定效应模型也称似不相关回归（SeeminglyUnrelatedRegression）αi反映了个体对总均值的偏离，表示不同个体有不同的截距项。βi反映了个体对总斜率的偏离，表示不同个体有不同的斜率。固定效应模型2、变系数的个体固定效应模型固定效应模型变系数的固定效应模型固定效应模型变系数的固定效应模型随机效应模型（三）个体随机效应模型（randomeffectmodel）把各个个体的值设为一套具有特定概率分布的随机变量，假定这些个体仅是从一个具有特定分布的总体中随机抽取出来的随机效应模型也可分为两大类：1、变截距的个体随机效应模型2、变斜率/系数的个体随机效应模型随机效应模型（三）个体随机效应模型（randomeffec随机效应模型1、变截距的个体随机效应模型：2、变斜率的个体随机效应模型vi是独立于解释变量的一个具有特定概率分布的随机变量，反映了个体成员的随机影响。随机效应模型1、变截距的个体随机效应模型：随机效应模型个体随机效应模型的参数估计：由于个体之间可能存在一定的相关，普通的最小二乘法不再是最有效的估计通常采用可行广义最小二乘法（feasiblegeneralizedleastsquare，FGLS）。随机效应模型个体随机效应模型的参数估计：面板数据模型固定效应模型与随机效应模型：固定效应模型中，假定每一个体有自己（固定）的截距；模型允许未被观测到的变量与任何已观测的自变量之间存在相关。随机效应模型中，假定这些个体是从一个更大的总体中随机抽取的。未被观测到的变量与所有自变量之间不相关。面板数据模型固定效应模型与随机效应模型：面板数据模型如何选择混合、固定或随机效应模型1、约束的F检验（比较混合模型和固定效应模型）判断：如果P大于0.05，可以认为约束模型（混合模型）成立；否则认为应采用无约束模型（固定效应模型）面板数据模型如何选择混合、固定或随机效应模型面板数据模型如何选择混合、固定或随机效应模型2、BP检验（拉格朗日乘数检验，比较混合模型和随机效应模型）原假设（H0）为：接受混合模型BP检验的统计量判断：如果P大于0.05，可以接受混合模型；

P小于0.05，可认为应采用随机效应模型面板数据模型如何选择混合、固定或随机效应模型面板数据模型如何选择混合、固定或随机效应模型3、Hausman检验（比较随机效应模型和固定效应模型）原假设（H0）为：接受随机效应模型Hausman检验的统计量判断：如果P大于0.05，可以接受随机效应模型；

P值小于0.05，可认为应采用固定效应模型面板数据模型如何选择混合、固定或随机效应模型面板数据模型如何选择混合、固定或随机效应模型4、更实际的考虑流行病学中个体追踪随访数据通常考虑随机效应模型更为合适省、市、区等的观察可考虑固定效应模型结合实际研究目的而定面板数据模型如何选择混合、固定或随机效应模型面板数据分析案例citycitymonthhfmtemperature北京110.78-4.8北京120.76-1北京136.144.1北京1424.7111.2北京1553.7721.7北京1667.7324.7上海216.825.3上海223.137.5上海2315.429.3上海2434.3712.8上海2540.9720.8上海2641.3923.8哈尔滨310.06-17.1哈尔滨320.01-15.5哈尔滨330.22-7.2哈尔滨341.134哈尔滨353.7416哈尔滨366.4825.5乌鲁木齐410.24-10.6乌鲁木齐420.04-12.4乌鲁木齐430.24-4.4乌鲁木齐442.058.7乌鲁木齐4513.5816.2乌鲁木齐4627.1622.5例1：四个城市1-6月份手足口病报告发病率与气温的关系面板数据分析案例citycitymonthhfmtemper面板数据分析案例首先进行因变量的变换原变量对数变换面板数据分析案例首先进行因变量的变换面板数据分析案例混合模型（pooled

model）的最小二乘估计：SAS程序：proc

reg;modellhfm=temp;run;面板数据分析案例混合模型（pooledmodel）的最小二面板数据分析案例最小二乘估计（pooled

model）结果：面板数据分析案例最小二乘估计（pooledmodel）结果面板数据分析案例自相关的Durbin-Watson检验：SAS程序：proc

reg;modellhfm=temp/dwdwprob;run;面板数据分析案例自相关的Durbin-Watson检验：面板数据分析案例Durbin-Watson检验结果：提示：可能存在空间或时间自相关面板数据分析案例Durbin-Watson检验结果：面板数据分析案例残差（Residual）分析：面板数据分析案例残差（Residual）分析：面板数据分析案例残差（Residual）分析面板数据分析案例残差（Residual）分析面板数据分析案例各城市的变化趋势图面板数据分析案例各城市的变化趋势图面板数据分析案例个体固定效应模型（最小二乘虚拟变量回归）SAS程序proc

panelprintfixed;idcitymonth;modellhfm=temp/fixone;run;面板数据分析案例个体固定效应模型（最小二乘虚拟变量回归）面板数据分析案例虚拟变量及自变量的估计结果：面板数据分析案例虚拟变量及自变量的估计结果：面板数据分析案例各城市截距：乌鲁木齐：-0.4091哈尔滨：-0.4091-0.67097=-1.08007北京：-0.4091+1.123651=0.714551上海：

-0.4091+1.210117=0.801017面板数据分析案例各城市截距：面板数据分析案例固定效应的约束F检验：模型拟合效果：与最小二乘回归相比，R2增加，MSE降低。面板数据分析案例固定效应的约束F检验：面板数据分析案例固定效应模型的残差分析面板数据分析案例固定效应模型的残差分析面板数据分析案例Pooled模型预测值固定效应模型预测值面板数据分析案例Pooled模型预测值面板数据分析案例个体随机效应模型SAS程序proc

panel;idcitymonth;modellhfm=temp/ranone;run;面板数据分析案例个体随机效应模型面板数据分析案例随机效应模型估计结果：随机效应模型的拟合结果面板数据分析案例随机效应模型估计结果：面板数据分析案例

Hausman检验和BP检验结果：P值仅略大于0.05，考虑选择固定效应模型面板数据分析案例Hausman检验和BP检验结果：面板数据分析案例拟合变斜率面板数据模型——固定效应模型似不相关回归（SUR）proc

syslinsur;modely1=x1;modely2=x2;modely3=x3;modely4=x4;run;面板数据分析案例拟合变斜率面板数据模型——固定效应模型面板数据分析案例似不相关模型（SUR）面板数据分析案例似不相关模型（SUR）面板数据分析案例普通最小二乘回归模型（OLS）面板数据分析案例普通最小二乘回归模型（OLS）面板数据分析案例SUR模型OLS模型面板数据分析案例面板数据分析案例拟合变斜率面板数据模型——随机效应模型proc

mixedcovtestmethod=MIVQUE0;classcity;modellhfm=temp/solution;randominttemp/subject=citytype=vc;run;面板数据分析案例拟合变斜率面板数据模型——随机效应模型面板数据分析案例随机系数分析参数估计结果面板数据分析案例随机系数分析动态面板数据模型如果现在的结果不仅受到当前因素的影响，还可能受到以往结果或因素的影响，可考虑动态面板数据模型。动态面板数据模型的参数估计通常采用广义矩估计（GMM）动态面板数据模型如果现在的结果不仅受到当前因素的影响，还可能动态面板数据模型例2：续例1数据，考虑到手足口病发病有可能受到上个月份手足口病例数的影响，采用动态面板数据进行分析分析：只考虑前一个月手足口发病率的影响，将因变量的一阶滞后项纳入模型动态面板数据模型例2：续例1数据，考虑到手足口病发病有可能受动态面板数据模型SAS程序

proc

panel;idcitymonth;claghfm(1)/out=hfmlag;run;proc

paneldata=hfmlag;instdepvar;modelhfm=hfm_1temp/gmmtwostepmaxband=5;idcitymonth;run;动态面板数据模型SAS程序动态面板数据模型主要结果：动态面板数据模型主要结果：空间面板数据模型常见空间模型1、空间滞后模型（SpatialLagModel，SLM）该模型包括解释变量X和空间滞后项Wy式中，W是空间权重矩阵。ρ是空间自回归系数，反映样本数据的空间依赖性，测量y对周边或相邻空间观测的影响，如果ρ≠0，提示存在区域之间的相互影响。ε表示随机误差。空间面板数据模型常见空间模型空间面板数据模型常见空间模型2、空间误差模型（SpatialErrorModel，SEM）该模型结合了一个标准的回归模型和一个随机误差ε的空间自回归模型。式中，W是空间权重矩阵，λ是空间误差相关的参数，ε表示随机误差。空间面板数据模型常见空间模型空间面板数据模型两种空间模型的选择如果我们确信每个个体的y值受到周围值的直接影响，可以考虑空间滞后模型。如果我们相信y并没有直接受到周围值的影响，而是因为某些我们没有考虑到的因素而导致模型的误差出现空间相关，此时可以考虑空间误差模型。空间面板数据模型两种空间模型的选择空间面板数据模型

空间面板数据模型

空间面板数据模型两种空间面板数据模型（SPDM）1、空间面板滞后模型2、空间面板误差模型空间面板数据模型两种空间面板数据模型（SPDM）空间面板数