对财政支出数据的回归

1、应用回归分析论文题目：影响财政支出因素的应用回归分析姓 名： 孙天阳 学 号： 20112981 摘 要截取了中国国家统计局统计年鉴中的财政支出数据，主要是对影响财政支出的因素进行分析，因为，财政支出的情况也在一定程度上反映了我国当前时段内的经济和社会情况。选取了各地区财政2008年的支出情况数据，对其做多元线性分析，建立回归模型，并通过对回归系数做显著性检验与逐步回归来分析数据。最后的结果就体现了各地财政支出与教育，国防，等的关系。关键词;逐步回归、显著性检验一、 模型1、多元线性回归的模型表示矩阵表达2、估计最小二乘估计要寻找经过整理得到回归参数的最小二乘估计就是3.显著性检验（1）、F检

2、验H0:1=2=p=0SST = SSR + SSE 当H0成立时服从（2）、回归系数的显著性检验H0j:j=0, j=1,2,p（,（X）-1）记 (X)-1=（cij) i,j=0,1,2, ,p构造t统计量 其中4、逐步回归(1)前进法前进法的思想是变量由少到多，每次增加一个，直至没有可引入的变量为止。首先分别对因变量y建立m个一元线性回归方程，并分别计算这m个一元回归方程的m个回归系数的F检验值，记为，选其最大者记为：给定显著性水平，若F（1,n-2），则首先将xj引入回归方程，为方便，设xj就是x1。接下来因变量y分别与 (x1，x2)，(x1，x3)，(x1，xm)建立m-1个二元

3、线性回归方程，对这m-1个回归方程中x2,x3,，xm的回归系数进行F检验，计算F值，记为，选其最大的记为：若F（1，n-3），则接着将xj引入回归方程。（2）后退法后退法与前进法相反，首先用全部m个变量建立一个回归方程，然后在这m个变量中选择一个最不重要的变量，将它从方程中剔除。设对m个回归系数进行F检验，记求得的F值为，选其最小者记为：给定显著性水平，若F（1,n-m-1），则首先将xj从回归方程中剔除，为方便，设xj就是xm。接着对剩下的m-1个自变量重新建立回归方程，进行回归系数的显著性检验，像上面那样计算出，如果又有F（1,n-（m-1）-1），则剔除xj，重新建立y关于m-2个自变

4、量的回归方程，依此下去，直至回归方程中所剩余的p个自变量的F检验值均大于临界值F(1,n-p-1)，没有可剔除的自变量为止。这时，得到的回归方程就是最终确定的方程。（3）逐步回归法逐步回归的基本思想是“有进有出”。具体做法是将变量一个一个引入，当每引入一个自变量后，对已选入的变量要进行逐个检验，当原引入的变量由于后面变量的引入而变得不再显著时，要将其剔除。这个过程反复进行，直到既无显著的自变量选入回归方程，也无不显著自变量从回归方程中剔除为止。这样就避免了前进法和后退法各自的缺陷，保证了最后所得的回归子集是“最优”回归子集。在逐步回归中需要注意的一个问题是引入自变量和剔除自变量的显著性水平值是

5、不相同的，要求进<出否则可能产生“死循环”。也就是当进出时，如果某个自变量的显著性P值在进与出之间，那末这个自变量将被引入、剔除、再引入、再剔除、，循环往复，以至无穷。二、数据分析1、指标的选取Y 一般预算支出 X1 一般公共服务 X2 国防X3 公共安全X4 教育X5 科学技术X6文化体育与传媒X7社会保障与就业X8医疗卫生X9环境保护X10城乡社区事务X11农林水事务X12交通运输2、结果数据及分析表一以上是参数估计的结果，建立回归方程(1方程的系数检验:表二方差分析的结果。P<0.05，证明高度显著，即变量全体x1-x12对y有显著影响，但是我们不能确定没一个变量x都对y有显

6、著性。由表一中我们可以看出诸如x2，x3，x4，x5,等都没有通过检验的t统计量及相应的P值，得到，并不是每一个单独的对因变量y都有显著性影响，这说明尽管回归方程通过了显著性检验，但也出现某些单个变量对y并不显著的情况。表三表三是样本相关阵的输出数据，x7等对y的贡献较大，而x2等比较小，所以我们需要剔除一些变量。（2）异方差和自相关性检验表四我们可以看到在其中的值都小于.05，所以不存在异方差性。我们看到dw=1.951，故认为不存在自相关性，我们在上表中可以看到有些变量的方差扩大因子较大，vif超过了10，所以我们认为变量间可能存在着多重共线性，我们在后面将进行变量的选择来剔除那些共线性变

7、量。（3后退法取显著性水平a=0.1000.变量剔除由此回归子集的最优方程就是：复决定系数是0.9910，c（p）是3.4128（4前进法，我们可以看到依次引入了x4 x7 x10 x12，此时我们就能写出最优模型的回归方程：复决定系数是0.9874.c（p）位10.4241.以上就是用后退法，前进法选择了变量建立了模型，建立了相对最优的回归方程，但是这两种方法有一定的缺陷，后退法中一旦某个变量被剔除了，就再也不能回到模型中，这样可能使找到的模型不是最优的。前进法相反，一旦变量进入，就不会再被剔除，但是这个变量可能在某个模型中不再显著，这样的模型同样不是最优的。所以我们再运用逐步回归来筛选，建

8、立一个最优模型。（5逐步回归显著性水平a=0.05这里我们先后由六步选择得到了这个最优子集模型，回归方程为复决定系数是0.9874.c（p）是10.4241.，最后发现vif都小于10，也都通过了显著性检验，所以这个模型是合理的。（6）主成分分析和岭回归以上结果我们可以得到选出的最优子集是x1 x4 x7 x8 x10 x11 x12.下面我们来做主成分回归；结果如下：以上我们看到三个主成分就能够反映91.88%的信息，所以我们取三个主成分就足够了。写出print1=0.422138x1+0.424429x4+0.30694x7+0.416599x8+0.30683x10+0.384503x1

9、1+0.332418x12,其他的表达式类似。做岭回归如下结果：我们可以写出K=0.2时的岭回归方程;y=1.0053x1+0.9349x4+1.40079x7+3.5135x8+1.89463x1+0.90681x11+3.80314x12也可以写出取三个主成分的回归方程。Y=0.8437x1+0.8956x4+1.53278x7+3.4658x8+2.35182x10+0.6661x11+4.3448x12。三、结果总结我们得到了回归方程：Y=-286388+0.67x1+1.18x3+0.59x4+1.99x5+2.80x6+1.31x7+2.22x8+1.38x9+1.23x10+1.

10、24x11+2.07x12，因此我们可以说，科学，医疗，交通运输与一般性支出的关系是显著的。为了降低财政支出，我们可以在加大对一些方面的支出，以便产生效益，来降低对于一般事项的支出。例如在全国社会保障支出最高的辽宁省（4699675）降低了一般性支支出较大的省往往也在教育，社会保障，社区事务，交通运输，医疗等方面支出较大，说明这些地区已经意识到了需要在这些方面加大支出，搞好建设，才能进一步削减财政支出。四、结语在目前，经济发展较快，而国家财政支出进一步加大，那么在社会建设中确实存在这一些问题，比如医疗卫生，教育，基础设施建设等。通过建立这样的一个回归模型，我们发现在经济快速增长的下，一些发面的

11、发展已经滞后了，那么如果我们不改善这种情况，财政将会越来越困难，所以在当前我们必须加强教育，社会保障和基础设施的建设，然后我们才能创造一个良好的财政收支环境，以利于经济和社会发展。总之，我们必须着眼未来，加强相关方面的建设，实现财政的良好循环。附录：程序；data sty;input y x1-x12;cards;1959285719626642210016046163162957112188661113820932851450513354688199838312177368034611160451867724595905556316584811416985286530180120105716

12、441915610978918501213853971415607610331881669631761852584513678193769819216674290023271372612023797635561324967151897728920714855751315017522348332301990753623498681763812718922183811714959642863101176310969343360379015071454573224339141590976447820640171536343162431915179598205796815170439916071774

13、930871253567215343483076503537711475833306360049022629940946996758390124817982172211149293337251322128261180122317424811386071950418803461341162821091998556595176456067853743107335729297615987261542300424044742765392498325661992008852396382286253716984485105985375148151452328122241082593916119878079

14、267149574132606281202743495223334966912228402507584686252789711170818359860132474927516860153224252884559260329151736674462315241148611695179536179042761599117212028865382208575637245565004819962464539902867928637593141522314286904651831939499177417559423417203241647125323415312002786894928625572377

15、883274652282005103844254736714787461367525518196186285711377159189197431468915985233292325628122425110929147427411402647509368042684560798801651210073018066692158775292720685781114061877951793603769195318377707820147015637603112185172704661346823514045917326125509929571333552178285052514041845860022

16、894598235300030045425436232281609340658961928013918624440264304389414616330234314547267585101353822209585943172720057861650276326580311065312564862841940230584252487281300395083840919485220217670063830691254733176522493010646405011127951311260126589325346231031208759874170661200158176375456957513608

17、303778568162902317769639320647033269132515566715436283452011467470878283294019260251058404195777812971100224536631057956563251221016214929246712897697876832797407230331393970584781106657435797085332258474242533556331449816817649434918637768058225139554118114304242897101601121390183448046143541534951

18、151279165439172266551636252929712461387672424792597174072948826936378234884215314363692812258150345730448955414356067914511364752242609310535861952922105379221842076158246807062297665129878178360107456867437540439146509412171074937685889911470238821712323373810651262419508176695279819224720110458725

19、845826751361777748608413911948380658964063247553015304708002903392134279000163548570682599566287013510341842071428520822613641324376392026490551714483180742455574783906587158867798146291854226891714396843361271358964245270318292569474319446515369845831504684694478351073358347976463132363595069415551

20、78203753488084396649891265567324661519548420483842440711045518197432460644234703319192769540553432657092537049017107317506733377845199610079522053710593638178443617200683001199213214835827657110905985864033046718888051431573338486589434;run;proc print;run;proc corr data=sty; var y x1 x2 x3 x4 x5 x6

21、x7 x8 x9 x10 x11 x12; run; Proc reg data=sty; model y=x1-x12/r dw vif; plot y*x1 y*x2 y*x3 y*x4 y*x5 y*x6 y*x7 y*x8 y*x9 y*x10 y*x11 y*x12;run; proc reg data =sty; model y = x1-x12 / selection = backward sle = 0.05 r; run; model y = x1-x12 / selection = forward sle =0.05 r;run;proc reg data =sty; mo

22、del y = x1-x12 / selection = stepwise sls = 0.05 sle = 0.2 r vif; run; proc standard data=sty mean=0 std=1; var y x1-x12; run; proc print; run;proc reg data=sty outtest=rrr graphics;model y=x1-x12/ridge=0 to 1 by 0.05;plot /ridgeplot;run;proc print data=rrr;run;proc princomp data=sty;var x1 x4 x7 x8

23、 x10-x12;run;proc princomp data=sty out=score n=3;var x1 x4 x7 x8 x10-x12;run;proc print data=score;run;proc reg data=sty outvif outest=a ridge=0 to 0.5 by 0.10;model y=x1 x4 x7 x8 x10-x12/pcomit=3 outvif;plot/ridgeplot;run;proc print data=a;run; 各地区财政支出 (2008年)一 般一 般国 防公共安全教 育科学技术预算支出公共服务4924849497

24、451372279807434111312851857711051864819592857196266422100160461631629571121886867724595905556316584811416985286530188166963176185258451367819376981921667413150175223483323019907536234986817638114545732243391415909764478206401715363421534348307650353771147583330636004902261180122317424811386071950418

25、803461341161542300424044742765392498325661992008852593916119878079267149574132606281202743324749275168601532242528845592603291517322085756372455650048199624645399028679281647125323415312002786894928625572377881137715918919743146891598523329232562811210073018066692158775292720685781114062704661346823

26、514045917326125509929571333228160934065896192801391862444026430438916502763265803110653125648628419402305841765224930106464050111279513112601265893377856816290231776963932064703326913251551297110022453663105795656325122101621493579708533225847424253355633144981101601121390183448046143541534951151279

27、294882693637823488421531436369281225815010537922184207615824680706229766512987814702388217123233738106512624195081766953806589640632475530153047080029033142852082261364132437639202649055171448968433612713589642452703182925694743363595069415551782037534880843966432460644234703319192769540553432651059

28、36381784436172006830011992132148358文化体育社会保障医疗卫生环境保护城 乡农林水事务交通运输与传媒和就业社区事务95513146460006527102580138515114191814142356334144079546111382093285145051335468819983831217736803461180120105716441915610978918501213853971415602900232713726120237976355613249671518977289207271892218381171495964286310117631096934336037316243191517959820579681517043991607177493087299409469967583901248179821722111492933372513282109199855659517645606785374310733572929762396382286253716984485105985375148151452328149522333496691222840250758468625278971117081866744623152411