北京联合大学统计学期末作业.docx

1、北京联合大学统计学期末作业北京联合大学题目：影响我国农村居民人均纯收入因素的分析姓名：胡锋学号：2013110444008学院：管理学院专业：信息管理与信息系统指导教师：康意2013年6月8日一、背景资料我国是农业大国，农村人口有9.4亿，占全国人口70%。农村居民收入问题关系着我国全面建设小康社会的大局 , 提高农村居民收入是建设社会主义新农村的重要内容之一 。改革开放以来， 我国农民人均纯收入从 1978 年的 133.57 元提高到 2013 年的8896元 ， 农民收入持续增长 , 但收入分配差距也在逐渐扩大 , 虽然目前分配差距仍处于警戒线之下 , 但已逼近警戒线 , 因而要高度重视

2、农村收入分配差距持续扩大的问题 。1978 的基尼系数虽然尚处于合理的区域内 , 但城乡收入差异仍在继续扩大 。1978 年 , 城乡收入差距为216 : 1 ,2002年以来，我国城乡收入比一直在“3”以上，并在2007年和2009年扩大到了“3.33”，城乡收入差距的居高不下一度成为我国经济发展的重要问题之一。近年来，随着我国农民收入增长迅速，城乡收入比正在回落，城乡差距正在逐步走向缩小。 城乡收入差距严重影响经济健康发展与和谐社会构建。农村居民收入过低，不仅影响农民生活水平提高，而且影响粮食生产和农产品供给；不仅制约农村经济发展，而且制约整个国民经济增长。因此，增加农村居民纯收入，缩小城

3、乡收入差距势在必行。1. 影响我国农村居民人均纯收入的主要因素有：（1）.财政支农支出。公共财政理论和实践表明，财政支农政策对于一国的农业增长至关重要。财政支农政策的积极作用主要表现在两个方面：其一，财政支农政策是国家调控农业生产进而影响农民收入 的一个基本工具；其二，财政支持能有效产品的外部性问题并具有规模经济的优势。(2).人均GDP值。它是重要的宏观经济指标之一，它是将一个国家核算期内（通常是一年）实现的国内生产总值与这个国家的常住人口（目前使用户籍人口）相比进行计算，得到人均国内生产总值。是衡量各国人民生活水平的一个标准。包括农村居民在内的人民生活水平提高了，可理解为农村居民纯收入增加

4、。(3).非农就业。由于农业产值在国内生产总值中的比重太低，农村中人地矛盾突出，进行农业生产成本高、收益低。农村居民投入第二第三产业有利于增加收入。城市化水平。它是衡量一个国家或地区的城市化最主要的指标，由城镇人口占总人口的比重表示。伴随着城市化，经济社会分工逐渐细化，更多的人从事经济、政治、文化事业，能转移大量的农村劳动力，从而提高农村居民纯收入。(4).城市化水平：又叫城市化率，是衡量城市化发展程度的数量指标，一般用一定地域内城市人口占总人口比例来表示2.数据选取数据选取1978到2013年的数据1.农村居民人均纯收入，在模型中用y表示2.影响农村居民纯收入的因素：1) 财政支农支出，用X

5、1表示2) 人均GDP值，用X2表示3) 非农就业系数，即第二、三产业的就业人数与就业总人数的比值，用X3表示4) 城市化水平，即城镇人口与总人口的比值，用X4表示5) 其他因素，我们将由各种原因未考虑到的或者无法度量的因素归入随机误差项，如农村居民信贷19782013年农村居民人均纯收入及其影响因素的统计数据表1：年份农村居民纯收入（元）y财政支农支出（亿元）x1人均GDP值（元）x2非农就业系数（二三产业就业人数/就业总人数）x3城市化水平(城镇人口/总人口）x41978133.676.953810.2950.1791979160.289.974190.3020.191980191.382

6、.124630.3130.1491981223.473.684920.3190.2021982270.179.885280.3190.2111983309.886.665830.3290.2161984355.395.936950.360.231985397.6101.048580.3760.2371986423.8124.39630.3910.2451987462.6134.1611120.40.2531988544.9158.7413660.4060.2581989601.5197.1215190.40.2621990686.3221.7616440.3990.2641991708.624

7、3.5518930.4030.2691992784269.0423110.4150.2741993921.6323.4229980.4360.2819941221399.740440.4570.28519951577.7430.2250460.4780.2919961926.1510.0758460.4950.30519972090.1560.7764200.5010.31919982162626.0267960.5020.33419992210.3677.4671590.4990.34820002253.4766.8978580.50.36220012366.40917.9686220.50

8、.37720022475.601102.793980.50.39120032622.201134.86105420.5090.40520042936.401693.79123360.5310.41820053254.901792.4141850.5520.43200635872161.35165000.5740.43920074140.363404.7201690.5920.45200847615955.522674.90.6040.45920095176.97253.123111.80.6190.47201059198579.7295240.6330.4832011697710408.635

9、198.60.6520.4992012791715000383540.6640.522013889613000425380.6750.537（数据来自于中国统计年鉴、中国金融统计年鉴)二 预测1.创建工作文件：选择【File】【New】【Workfile】（或者在命令窗口中输入：create）如图1所示：图12.查看序列Y的图形图2 ：19782013年农村居民纯收入指数时序图观测序列的时序图，可知序列具有曲线趋势，先用低阶差分法提取确定性信息。图3：对序列Y进行一阶差分运算图4：一阶差分运算公式图5：一阶差分序列(dz)图6：一阶差分曲线图可以看出一阶差分序列仍然具有线性趋势，继续进行差分分

10、析；二阶差分的操作和上面一阶差分的操作一样。图7：二阶差分序列(dz1)曲线图从二阶差分序列曲线图可以看出二阶差分序列中没有曲线趋势，二阶差分提取了曲线趋势。以下图7、图8是对二阶差分序列进行自相关分析图7 图8图9：二阶差分序列相关图可以看出二阶差分序列具有短期相关性的特征，无确定性信息，为平稳序列。3.arma模型定阶由图9看出，偏自相关系数在k=4后很快趋于0即4阶截尾，尝试拟合AR（4）；自相关系数在k=2处显著不为0，当k=4时在2倍标准差的置信带边缘，可以考虑拟合MA（2）或MA（4）；同时可以考虑ARMA（4，2）模型等。4.arma模型建立arma（p,q）由模型定阶发现，p可

11、能等于4，q可能等于4或2，我们根据各种组合来选择最优模型，在工作文件主窗口点击Quick/Estimate Equation ，在Equation Specification 对话框中填入dz1 ar(1)ar(2)ar(3)ar(4)ma(1)按回车，即得到参数估计结果，如图10图10：ARMA（4，1）模型估计结果由参数估计结果看出，各系数均不显著，说明模型并不适合拟合ARMA(4，1)模型。经过进一步筛选，逐步剔除不显著的滞后项或移动平均项，最后得到如下ARMA(2，1)模型：如图11图11：ARMA（2，1）模型估计结果综上可见，我们可以对同一个平稳序列建立多个适合模型，但比较AIC

12、和SC的值，以及综合考虑其他检验统计量，考虑模型的简约原则，我们认为ARMA（2，1）模型是较优选择。5模型预测预测方法常用有两种：Dynamicforecast和Staticforecast，前者是根据所选择的一定的估计区间，进行多步向前预测；后者是只滚动的进行向前一步预测，即每预测一次，用真实值代替预测值，加入到估计区间，再进行向前一步预测。在这里我选择Dynamicforecast，点击ok，出现如图12预测对话框：图12图13序列动态预测图预测值存放在dz1f序列中，此时我们可以观察原序列dz1和dz1f之间的动态关系，同时选中dz1和dz1f，击右键，点open/asgroup，然后

13、点击view/graph/line，则出现图14，动态预测值几乎是一条直线，说明动态预测效果很不好。图14 动态预测效果图进行静态预测，见图15，预测值存放在dz2f中，做dz1和dz2f之间的静态分析如图16，可以看出静态预测效果不错。图15：静态预测图图16：静态预测效果图三回归分析 因为数据是时间序列数据，所以在回归分析之前需要做平稳性检验。为了确定没有随机趋势或确定趋势，否则将会产生“伪回归”问题。伪回归是说，有时数据的高度相关仅仅是因为二者同时随时间有向上或向下的变动趋势, 并没有真正联系。这样数据中的趋势项,季节项等无法消除, 从而在残差分析中无法准确进行分析。 平稳性检验的方法可

14、以用Uint Root Test 检验,。消除趋势可以用差分法（比如一阶差分），模型也只有通过平稳性检验才有统计分析的意义。1. 新建工作文件2.计算变量间的相关系数3时间序列的平稳性检验（1）观察Y序列趋势图在eviews中得到时间序列趋势图，在quick菜单中单击graph，在serieslist对话框中输入序列名称Y，其他选择默认操作。图形表明序列随时间变化存在上升趋势。（2）对原序列进行ADF检验得到序列的ADF平稳性检验结果，检测值2.587946大于所有临界值，则表明序列不平稳。以此方法，对各时间序列依次进行ADF检验，将检验值与临界值比较，发现所有序列的检验值均大于临界值，表明各

15、原序列都是非平稳的。（3）时间序列数据的一阶差分的ADF检验得到序列的ADF平稳性检验结果，检测值0.234432大于所有临界值，则表明序列不平稳。以此方法，对各时间序列依次进行ADF检验，将检验值与临界值比较，发现所有序列的检验值均大于临界值，表明各原序列的一阶差分都是非平稳的。（4）时间序列数据的二阶差分的ADF检验得到序列的ADF平稳性检验结果，检测值-6.118138远小于所有临界值，则表明序列二阶差分平稳。以此方法，对各时间序列的一阶差分依次进行ADF检验，将检验值与临界值比较，发现所有序列的检验值均小于临界值，表明各序列二阶差分都是平稳的。由此可知，以上时间序列均为二阶单整时间序列

16、。2.异方差检验（原假设H0：不存在异方差）根据上图得到样本回归方程：lnY=-1058.087+（-0.000363）lnX1+0.166601X2+3654.220lnX3+56.33930lnX4t=(-4.304796） （-0.507434） （27.84190） （ 6.075560） （0.172339） R2=0.995392 F= 1674.272 D.W= 0.689172异方差检验根据上述的检验结果，Obs*R-squared=30.22279，P=0.0984，该P值相比较显著性水平5%要大，因此可以接受原假设 ，此时，基本不存在异方差性。3.残差自相关性检验检验得到D

17、W=0.689172，给定显著性水平=0.05，在n=36,k=4时，查Durbin-Waston表得有上限临界值dL= 1.30，下限临界值dU=1.65,此时0 D.W dL，故残差存在一阶正自相关。4消除一阶正自相关（运用迭代法）在命令窗口中输入：ls y c x1 x2 x3 x4 ar(1)，做一阶迭代输出结果如下图所示由图形可知，此时的DW=1.463638，为消除正自相关5。消除二阶正自相关，在命令窗口ls y c x1 x2 x3 x4 ar(1) ar(2)，输出结果如下图由图形可知DW=1.756138，因为dUDW4- dU,所以此时消除了正自相关性四实验总结在实验过程中

18、，首先运用arma模型分析法对农村居民纯收入进行预测。再对数据做回归分析，分为3 步。平稳性检验，运用ADF方法检验，如果数据不好，可以进行1阶或者2阶差分处理；异方差检验，利用参数方程，得到检验值理想，确定数据不存在异方差性；残差自相关性检验，通过分析得到残差存在正自相关性，然后通过迭代法消除正自相关性。通过对模型的诊断和调整，得出模型公式为lnY=-1058.087+（-0.000363）lnX1+0.166601X2+3654.220lnX3+56.33930lnX4说明农村居民人均纯收入与财政支农支出、非农就业系数人均GDP、城市化水平存在一定的关系和影响。做二元回归分析，在数据处理方面运用的方法不是很多，主要是取对数法，诸如权值参数法、广义加性法等不适合二元分析。涉及Beta系数、弹性系数、偏相关系数等还不清楚如何处理，要求更深的知识，还有待于提升。