社会消费品零售总额的SAS报告分析.doc

基于sas分析湖南省社会消费品零售总额影响因素的研究【摘要】本文旨在分析1978-2008年改革开放以来，我国社会消费品零售总额变动情况，影响其变动的因素。首先，我们提出了关于收入和消费的主要理论观点，然后再引入其他有关变量，进而建立了理论模型。然后，收集了相关的数据，利用sas软件对计量模型进行了参数估计和检验，并加以修正。最后，我们对所得的分析结果作了经济意义的分析，详细剖析其成因，进一步进行预测，并相应提出一些政策建议。【关键词】社会消费品总额 财政支出 城镇居民家庭人均可支配收入 居民消费 一引言为促进中部地区（山西 江西 河南 湖北 湖南和安徽六省）经济快速发展，中共中央提出“中部崛起”的战略。中部六省对中部崛起战略非常拥护，并抱有很高的期望。中部六省希望中共加大对中部崛起的政策支持力度，支持中部建设全国粮食核心主产区，支持中部建立前进制造业基地，支持中部加快老工业基地改造，资源型城市转型和国有企业改革，支持中部解决交通设施的薄弱环节，支持中部治理生态和环境，支持中部教育卫生事业发展，支持中部减轻财政负担。二 理论背景为了加强在中部崛起竞争力湖南省加快了经济建设的步伐。社会消费品零售总额所计量的是各种经济类型的商业由于经济的发展和社会的进步，特别是社会主义市场经济的建立，商品生产和商品交换的领域进一步扩大，用已确立和描述各类消费品市场对居民和社会集团出售商品总和的商品零售额指标的口径范围也作了相应的调整。社会消费品零售总额的增长体现了该省在市场经济中的竞争力，因此对影响社会消费品零售总额的因素作出分析和验证有着重要的必要性。 社会消费品零售总额是指各种经济类型的批发零售贸易业、餐饮业、制造业和其他行业对城乡居民和社会集团的消费品零售额和农民对非农业居民零售额的总和。这个指标反映通过各种商品流通渠道向居民和社会集团供应的生活消费品来满足他们生活需要，是研究人民生活、社会消费品购买力、货币流通等问题的重要指标。 研究影响社会消费品零售总额的因素可以进一步证实其反映出了人民生活水平，生活质量的提高；也可以反映出湖南省经济的脉象，便于零售商指定营销战略和企划。三模型建立本模型是通过对湖南1978-2008年数据进行研究。1. 确定模型所包含的变量，因变量（被解释变量）为社会消费品零售总额（y单位：亿元）自变量（解释变量）有3个财政支出（x1单位：亿元），城镇居民家庭人均可支配收入（x2单位：元） ，居民消费（x3单位：亿元）。2. 确定模型的数学形式，我们先假设数学模型为： 四 实验方法本次研究主要是分析自变量和因变量的关系以及找出合适的模型，并且运用时间序列的方法对社会消费品零售总额进行预测分析。我们先用sas对我们假定的模型对变量x1，x2，x3分析，看看它们与y的关系（运用的方法的逐步回归和主成分法）然后在对社会消费品零售总额进行时间序列分析，最后预测未来3年的数据。五 数据来源与sas操作1 数据来源与cnki网2 数据引人年份社会消费品零售总额y（亿元）财政支出x1（亿元）城镇家庭可支配收入x2（元）居民消费x3（亿元）197854.524.46323.8893.93 197964.825.17365.07111.48198076.323.71475.92125.42 198186.721.39505.08150.27198294.823.26519169.051983106.725.31564190.191984123.630.04645213.961985156.540.09760.8240.71 1986179.554.29904.44265.54 1987212.555.931017.8298.85 198827664.891254.97364.37 1989297.974.231492.61393.49 1990299.180.081591.45454.90 1991339.788.581783.24500.84 1992398.799.12166.5582.91 1993493.5132.032816.5706.55 1994669151.493887.6883.46 1995849.23173.944699.21108.78 1996960.8217.745052.11409.44 19971056.4230.825209.71540.67 19981141.3273.645434.31601.83 19991246.6313.125815.41690.66 20001383.7347.836218.71782.77 20011531.8431.76780.61871.64 20021701.8533.026958.61992.05 20031885.6573.74537674.22016.88 20042149.6719.548617.482269.69 20052459.1873.429523.972553.41 20062834.221064.5210504.672872.58 20073356.491357.0312293.543217.29 20084222.61765.224913821.163619.45 3. sas操作（1） 逐步回归法逐步回归的基本思想是：首先在待选的m个变量中选择一个对因变量影响最大的自变量，这可以通过因变量与每一自变量进行回归得到的f值来判断。若最大的f值在给定的显著性水平下是显著的，则该变量被选中，否则选元结束。第二步，在剩下的m-1个变量中再选择一个变量加入到模型中，这可以通过偏f检验来判断。若最大的偏f统计量在给定的显著性水平下是显著的，则对应的自变量则被加入到模型中，否则不再加入其它变量，选元结束。第三步，对已在模型中的每个变量进行显著性检验。若检验不显著，则去掉该变量，再重复进行第二步和第三步，否则保留该变量，重复进行第二步和第三步，这一过程一直进行到待选的全部自变量根据给定的显著性水平没有一个再能被选入模型或排除出刚构成的回归模型为止。输入程序：data exp3; input y x1 x2 x3;cards;54.524.46323.8893.93 64.825.17365.07111.4876.323.71475.92125.42 86.721.39505.08150.2794.823.26519169.05106.725.31564190.19123.630.04645213.96156.540.09760.8240.71 179.554.29904.44265.54 212.555.931017.8298.85 27664.891254.97364.37 297.974.231492.61393.49 299.180.081591.45454.90 339.788.581783.24500.84 398.799.12166.5582.91 493.5132.032816.5706.55 669151.493887.6883.46 849.23173.944699.21108.78 960.8217.745052.11409.44 1056.4230.825209.71540.67 1141.3273.645434.31601.83 1246.6313.125815.41690.66 1383.7347.836218.71782.77 1531.8431.76780.61871.64 1701.8533.026958.61992.05 1885.6573.74537674.22016.88 2149.6719.548617.482269.69 2459.1873.429523.972553.41 2834.221064.5210504.672872.58 3356.491357.0312293.543217.29 4222.61765.224913821.163619.45 ;proc reg graphics;model y=x1 x2 x3/selection=stepwise sls=0.05 sle=0.2 r;plot student.*p.=*;run;该程序告诉我们先建立一个临时数据集（一般临时的数据集在逻辑库里的work中）名叫exp3；里面有4个变量y x1 x2 x3.其中proc reg graphics 是表示回归作图“selection=stepwise”表示逐步回归（选元，并给出回归结果），“sls=0.05”表示变量保留在模型中的显著性水平为0.05，而“sle=0.2”则表示变量选入到模型中的显著性水平为0.2，plot语句画标准残差与预测值的残差图。逐步回归结果图和残差图图1从图中对因变量y进行逐步回归结果表明，第一步，变量x2被加入到模型中，说明三个自变量分别关于y回归，其中由x2与y回归得到的f统计量值最大,为1087.10，对应的概率为0.0001，在0.05的显著性水平下通过检验。第二步，分别把变量x1和x3加入到刚建立的模型中，并分别计算其偏f检验值，结果表明x1的偏f统计值大于x3的偏f统计值，为638.91对应的概率为0.0001，小于进入模型的显著性水平0.20，所以x1被选入到模型中，此时模型中已有expend和ads两个变量。第三步，对x3和x1分别计算其偏f统计量值及其对应的概率(程序繁多就不一一输出)，结果两个变量的偏f检验对应的概率都小于保留在模型中的显著性水平0.05，因此两个变量都被保留下来。第四步，重复第二步和第三步，即把最后一个变量加入到由x2和x1对y进行回归的模型中，结果两个变量的偏f检验对应的概率也都小于保留在模型中的显著性水平0.05，因此两个变量也都可以被保留下来。从残差图中和逐步回归表中可以得到回归方程如下：统计量值：（5.87） （ 892.2） （8.72） (12.84) f=12165.4对应概率：（0.0223）（0.0001） （0.0064） （0.0013） (0.0001)r2= 0.9993，c（p）=4，mse=960.58764，ssr=25936.以上模型中，截据项的t在0.05的显著性水平下通过检验。通过以上分析我们可以得出：1. 财政支出，城镇家庭可支配收入以及居民消费都和社会消费品零售总额有显著的关系。2. 财政支出的系数表明财政支出每增加1亿元，社会消费品零售总额将增加1.25381亿元，而城镇家庭可支配收入系数表明城镇家庭可支配收入每增加1元，社会消费品零售总额将增加0.06827亿元，居民消费的系数表明居民消费每增加1亿元，社会消费品零售总额将增加0.28841亿元。（2） 多重共线性多重共线性是指在多元回归分析中自变量之间存在线性关系，用数学术语来表达就是系数矩阵的秩p+1，即r（）p+1，换句话说，某一个自变量可以被其它自变量线性表示，此时称自变量之间完全共线。在sas中输入代码：proc reg corr;model y=x1 x2 x3/ vif collin collinoint;run;在上面的程序中，proc reg调用回归分析过程，并按model语句给出的模型结构进行分析，其中选项corr要求计算变量之间的相关系数， vif要求计算方差膨胀因子，而collin和collinoint则要求进行多重共线诊断，两者的区别是前者没有对截据项进行调整，而后者进行了调整。当截据项经检验是显著时，应考察由collinoint输出的结果，否则参看由collin输出的结果(由于我们前面已经检验了截距项是显著的，所以我们这里考察collinoint输出的结果)。程序运行结果如下：变量之间线性回归及多重共线诊断结果从变量之间线性回归及多重共线诊断结果中最上面的相关系数阵可以看出，x1、x2、x3与y都是正相关，并且x1、x2,x3都与y高度正相关，相关系数分别为0.9840,0.9869和0.9823。在自变量之间，x2与y的相关系数为0.9869最高，而x2与x3的相关系数高达0.9974（x2为城镇家庭可支配收入，x3为居民消费）x1和x2以及x1和x3的相关系数也很好高其值分别为0.9445和0.9351；上述说明x2与y存在高度的线性关系，也说明x2和x3，x1存在高度相关性，且x1和x3也存在高度相关性(可以说明存在多重共线性)；这与我们从经济理论得到的判断是一致的。输出表中的方差膨胀因子(vaiance inflation)和最下面的共线性诊断（因为截据项检验显著）也说明了这一点。其中条件数为36.60466，而在变量x2和x3上的方差比率分别为0.99326和0.98746，远超过50%，说明两变量高度共线。如果不对多重共线进行处理而直接进行线性回归，根据输出分析表可以写出以下方程：检验统计量：（-2.42） (29.87) (2.95) (3.58) f=12165.4对应概率: （0.0223）(0.0001) (0.0064) (0.0013) (0.0001)r2=0.9993 r2-adj=0.9992 root mse=30.99335上述分析中我们发现的x2和x3存在高度相关性，虽然x2和x3的t检验显著。由于自变量x2与x3高度相关,为了消除这一影响，通常可以去掉一个不重要的因素，然后就剩下的两个自变量对因变量进行回归。考虑到x2和x3对y都很重要，为保留这两个变量，以下我们将reg过程就三个自变量对y做主成分回归。输入代码：proc reg data=exp3 outest=result1 outvif;model y=x1 x2 x3/pcomit=1;run;proc print data=result1;run;在上述程序中，model语句后的选项“pcomit=1”表示去除最后一个主成分，用剩下的p-1个主成分作为自变量关于因变量回归。本例中有三个自变量，所以主成分有三个，去除一个后还有两个。用这两个主成分对进口额进行回归，回归的结果存放到由proc reg语句规定的选项“outest=”的数据集result1中。程序运行结果如下：由图表的第三行可以写出主成分回归方程如下：从方程中我们可以得出如下结果：财政支出的系数表明财政支出每增加1亿元，社会消费品零售总额将增加1.24843亿元，而城镇家庭可支配收入系数表明城镇家庭可支配收入每增加1元，社会消费品零售总额将增加0.07575亿元，居民消费的系数表明居民消费每增加1亿元，社会消费品零售总额将增加0.2645亿元。（3）.预测我们预测y（预测未来三年的社会消费品零售总额）可以通过时间序列的方法进行预测。先通过sas的文件菜单点击导入数据，并且把导入的数据命名为exp4(一般导入的数据在逻辑库的work里面)。对数据进行平稳性检验。输入代码：proc gplot data=work.exp4;plot y*date;symbol c=red i=jion v=star;run;proc arima data=work.exp4;identify var=y nlag=12;run;第一个代码表示对该序列绘制时序图用于观测平稳性，其中plot y*date的含义是规定将输出以date为横坐标，以y为纵坐标的曲线图；symbol c=red i=jion v=star该语句的作用是规定所绘制图形的格式。此图形中各点的形状为星号，各点是通过直线连接的，曲线的颜色为红色。第二个代码表示对该序列输出自相关图，偏自相关图和随机性检验结果图；绘制图如下：从图中可以看出具有明显的趋势性（有点成曲线指数增加的趋势）为非平稳模型。从纯随机性结果图可以看到该模型p值检验显著，综上所述该序列为非平稳非白噪声序列；由于出现了明显的指数趋势性的缘故，故此我们对原序列指数趋势转化为线性趋势，然后在进行差分以消除线性趋势；data work.exp4;set work.exp4;y1=log(y);run; proc gplot data=work.exp4;plot y1*date;symbol c=red v=star i=join;run;该代码表示对y进行log运算产生新值为y1；并绘制图形。图形可以看出经过指数变换后该序列成线性趋势；然后在进行差分以消除线性趋势data work.exp4;set work.exp4;dif=dif(y1);proc gplot data=work.exp4;plot dif*date;symbol c=red v=star i=join;run;该代码表示对新序列（y1）进行一阶差分,并绘制图形。输出结果图为: 从图像可以看出做完一阶差分序列趋于平稳，为此我们来观察纯随机检验图。proc arima data=exp4;identify var=y1(1);run;对模型进行自相关图的产生和纯随机检验结果图如下；从纯随机性结果图可以看到该模型p值(0.0320.05)故检验显著，综上所述该序列为平稳非白噪声序列；我们下一步进行模型识别与优化，定阶：proc arima data=exp4;identify var=y1 minic p=(0:7) q=(0:7);run;该代码表示自相关延迟阶数p在0和7之间，骗自相关系数阶数q在0和7之间的所有arma(p,q)模型的bic信息量，并显示其中的最小值。可以看到arma（7,7）的bic信息量最小，其值为-9