1、2012級(jí) 金融數(shù)學(xué)專業(yè)應(yīng)用統(tǒng)計(jì)課程期末結(jié)課題目：1950-2010中國(guó)GDP數(shù)據(jù)分析作者 ：，學(xué)號(hào)： 1222150284任課教師：王沛成績(jī)： 日期：2015年 6 月 30 日-1950-2010 年中國(guó) GDP 的數(shù)據(jù)分析摘要：GDP 不僅能夠反映一個(gè)國(guó)家（或地區(qū)）的生產(chǎn)規(guī)模，而且能夠反映這個(gè)國(guó)家的產(chǎn)業(yè)結(jié)構(gòu)，如三大產(chǎn)業(yè)在整個(gè)國(guó)民經(jīng)濟(jì)中所占的份額。本將運(yùn)用 SAS技術(shù)，采用正態(tài)分布檢驗(yàn)、方差分析、相關(guān)分析以及回歸分析等方法對(duì)第一、二、三產(chǎn)業(yè)進(jìn)行分析，從而更加深刻的理解第一、二、三產(chǎn)業(yè)之間的相互關(guān)系和影響以及其對(duì) GDP 的貢獻(xiàn)。本文基于 SAS 以我國(guó) 1950-2010 年的國(guó)內(nèi)生產(chǎn)總值

2、為基礎(chǔ)，對(duì)數(shù)據(jù)進(jìn)行繪圖分析、模型識(shí)別、模型估計(jì)及預(yù)計(jì)。關(guān)鍵字：三大產(chǎn)業(yè) GDP正態(tài)性檢驗(yàn) 方差分析 相關(guān)分析 回歸分析1.前言國(guó)內(nèi)生產(chǎn)總值(Gross Domestic Product,簡(jiǎn)稱 GDP)是指在一定時(shí)期內(nèi)(一個(gè)季度或一年)，一個(gè)國(guó)家或地區(qū)的經(jīng)濟(jì)中所生產(chǎn)出的全部最終產(chǎn)品和勞務(wù)的價(jià)值，常被公認(rèn)為衡量國(guó)家經(jīng)濟(jì)狀況的最佳指標(biāo)。它不但可反映一個(gè)國(guó)家的經(jīng)濟(jì)表現(xiàn)，更可以反映一國(guó)的國(guó)力與。GDP 是最重要的宏觀經(jīng)濟(jì)統(tǒng)計(jì)指標(biāo)之一，它是人們了解和把握一個(gè)國(guó)家(或地區(qū))的宏觀經(jīng)濟(jì)運(yùn)行狀況的有效工具，是制定經(jīng)濟(jì)政策的重要依據(jù)，也是檢驗(yàn)經(jīng)濟(jì)政策科學(xué)性和有效性的重要從生產(chǎn)者角度來(lái)說(shuō)，分析了解三大產(chǎn)業(yè)之間的結(jié)構(gòu)

3、是非常重要的。因此，正態(tài)性檢驗(yàn)因?yàn)椋疚牡臉颖緮?shù)量小于 2000，為了更好的了解自 1952 至 2012 年三大產(chǎn)業(yè)每年的發(fā)展情況，對(duì)三大產(chǎn)業(yè)進(jìn)行了正態(tài)性檢驗(yàn)。第一產(chǎn)業(yè)正態(tài)性檢驗(yàn).-從上圖的正態(tài)檢驗(yàn)看，Shapiro-Wilk 的 W 值為 0.695947 太小，其檢驗(yàn)的顯著性水平小于值 0.05，即原假設(shè) H0 假設(shè)，表明數(shù)據(jù)基本上不屬于正態(tài)分布。因此，可以得出第一產(chǎn)業(yè)自 1952-2012 年，其數(shù)據(jù)不屬于正態(tài)分布。（如果數(shù)據(jù)呈現(xiàn)正態(tài)性，正態(tài)概率圖中的*點(diǎn)應(yīng)該覆蓋+號(hào)，因?yàn)?號(hào)線表示數(shù)據(jù)實(shí)際的“分布線”，“+”號(hào)線表示理論上要求的“正態(tài)線”）因此，得出第一產(chǎn)業(yè)第二產(chǎn)業(yè)正態(tài)分布檢驗(yàn).正態(tài)

4、性檢驗(yàn)檢驗(yàn)統(tǒng)計(jì)量p 值Shapiro-WilkW0.6959470.2576971.3213097.077739Pr W0.00010.01000.0050 DCramer-von MisesW-SqPr W-SqAnderson-DarlingA-SqPr A-Sq-從上圖的正態(tài)檢驗(yàn)看，Shapiro-Wilk 的 W 值為 0.606243 太小，其檢驗(yàn)的顯著性水平小于值 0.05，即原假設(shè) H0 假設(shè)，表明數(shù)據(jù)不屬于正態(tài)分布。（第二產(chǎn)業(yè)正態(tài)分布檢驗(yàn)圖）（如果數(shù)據(jù)呈現(xiàn)正態(tài)性，正態(tài)概率圖中的*點(diǎn)應(yīng)該覆蓋+號(hào)，因?yàn)?號(hào)線表示數(shù)據(jù)實(shí)際的“分布線”，“+”號(hào)線表示理論上要求的“正態(tài)線”）因此，第二

5、產(chǎn)業(yè)的自 1952-2012 年的數(shù)據(jù)為非正態(tài)分布。的自 1952-2012 年的數(shù)據(jù)為非正態(tài)分布。.正態(tài)性檢驗(yàn)檢驗(yàn)統(tǒng)計(jì)量p 值Shapiro-WilkW0.6062430.3055981.8588549.646728Pr W0.00010.01000.0050 DCramer-von MisesW-SqPr W-SqAnderson-DarlingA-SqPr A-Sq-第三產(chǎn)業(yè)正態(tài)分布檢驗(yàn)（第三產(chǎn)業(yè)正態(tài)檢驗(yàn)數(shù)據(jù)）從上圖的正態(tài)檢驗(yàn)看，Shapiro-Wilk 的 W 值為 0.575966 太小，其檢驗(yàn)的顯著性水平小于值 0.05，即原假設(shè) H0 假設(shè)，表明數(shù)據(jù)不屬于正態(tài)分布。（第三產(chǎn)業(yè)正態(tài)

6、檢驗(yàn)圖）（如果數(shù)據(jù)呈現(xiàn)正態(tài)性，正態(tài)概率圖中的*點(diǎn)應(yīng)該覆蓋+號(hào)，因?yàn)?號(hào)線表示數(shù)據(jù)實(shí)際的“分布線”，“+”號(hào)線表示理論上要求的“正態(tài)線”）因此，得出第三產(chǎn)業(yè)的自 1952-2012 年的數(shù)據(jù)為非正態(tài)分布。.正態(tài)性檢驗(yàn)檢驗(yàn)統(tǒng)計(jì)量p 值Shapiro-WilkW0.5769660.3050422.05968710.49023Pr W0.00010.01000.0050 DCramer-von MisesW-SqPr W-SqAnderson-DarlingA-SqPr A-Sq-相關(guān)系數(shù)分析為了了解第一、二、三產(chǎn)業(yè)與國(guó)內(nèi)生產(chǎn)的相關(guān)性，于是進(jìn)行了其進(jìn)行了相關(guān)系數(shù)分析，同時(shí)，通過(guò)對(duì)三大產(chǎn)業(yè)與國(guó)內(nèi)生產(chǎn)總值

7、的相關(guān)性的分析，為后面做關(guān)于國(guó)內(nèi)生產(chǎn)總值的回歸奠定基礎(chǔ)。分析結(jié)果如下：，國(guó)內(nèi)生產(chǎn)總值與第一產(chǎn)業(yè)的相關(guān)系數(shù)為：0.98746，且“當(dāng) H0： Pho=0 時(shí)，Prob |r|”的概率值 0.0001 |r|”的概率值 0.0001 |r|”的概率值 0.0001 |t|ercept1144.577887 237.93222034.81.00011.5461990.0509152 30.37 Fseco1 2051904493 2051904493 922.22 .00011315044549315044549 141.60 Fseco1 870902778558 870902778558 391

8、425 .00011 315044549.26 315044549.26 141.60 F模型2 871217823107 435608911554 195783 .000158 129047342.19 2224954.175760 871346870449誤差校正合計(jì)-結(jié)果分析：（1）方差分析的統(tǒng)計(jì)量M的 Sum of Square ：回歸平方和（SSR）=871217823107，誤差平方和（SSE）= 129047342.19，總平方和（SST）= 871346870449，三者的度分別為：2、58、60?；貧w平方和得均值（MSR）=435608911554，誤差平方和均值（MSE）=

9、 2224954.1757。實(shí)際的 F 值= 195783，其概率0.0001 |t|:該為“回歸系數(shù) B0 0 ”的 t 值的檢驗(yàn)值，變量第二產(chǎn)業(yè)（seco）的 t 檢驗(yàn)的概率 P 值0.05=因此，變量第二產(chǎn)業(yè)（seco）的回歸系數(shù)有效。同理，變量第三產(chǎn)業(yè)（thir）的 t 檢驗(yàn)的概率 P 值0.05=，非常顯著。（8）（9）Standard Error：這是回歸系數(shù)估計(jì)的標(biāo)準(zhǔn)誤差。因此截距項(xiàng)的標(biāo)準(zhǔn)誤差 237.9322203，誤差大。而第二產(chǎn)業(yè)（seco）的截距為 0.0509152，第.thir0.6498090.0546085 11.90.0001-三產(chǎn)業(yè)（thir）的截距為 0.0

10、546085,誤差小，說(shuō)明效果好。（10） Estimate：回歸方程式及回歸系數(shù)的估計(jì)值，得到第二產(chǎn)業(yè)的回歸系數(shù)為：1.546199，第三產(chǎn)業(yè)（thir）的回歸系數(shù)為：0.649809 截距項(xiàng)為：1144.577887，由此，最終得到回歸方程為：Y 1144.577887 1.546199X1 0.649809X2四、總結(jié)我國(guó)的三大產(chǎn)業(yè)包括農(nóng)業(yè)（種植業(yè)、林業(yè)、牧業(yè)、副業(yè)和漁業(yè)），工業(yè)和建筑業(yè)，以及除了第一二產(chǎn)業(yè)以外的其他部門，通過(guò)對(duì)其均值的比較，發(fā)現(xiàn)第二產(chǎn)業(yè)，即工業(yè)和建筑業(yè)在整個(gè)國(guó)民經(jīng)濟(jì)中依然占主導(dǎo)地位，原因要?dú)w咎于我國(guó)的國(guó)情，我國(guó)自開(kāi)放以后，各個(gè)產(chǎn)業(yè)有了巨大的改變，但是，第二產(chǎn)業(yè)為主導(dǎo)的產(chǎn)

11、業(yè)結(jié)構(gòu)至 1998 年底，依舊沒(méi)有動(dòng)搖。通過(guò)對(duì)三大產(chǎn)業(yè)自 1952 年至 2012 年的數(shù)據(jù)的正態(tài)性檢驗(yàn)，發(fā)現(xiàn)，它并不是呈正態(tài)分布，而是隨著時(shí)間的遞增而遞增，拋除通貨膨脹，可以得出，三大產(chǎn)業(yè)正呈現(xiàn)出良好的發(fā)展趨勢(shì)，說(shuō)明的國(guó)家正在健康的成長(zhǎng)。最后，進(jìn)行了方差分析，相關(guān)系數(shù)分析以及回歸分析，得出，三大產(chǎn)業(yè)之間有很好的相關(guān)性，他們對(duì)國(guó)內(nèi)生產(chǎn)總值的貢獻(xiàn)是無(wú)可厚非的，第一產(chǎn)業(yè)是第二三產(chǎn)業(yè)的基礎(chǔ)，而一二產(chǎn)業(yè)為第三產(chǎn)業(yè)的條件，第三產(chǎn)業(yè)發(fā)展促進(jìn)第一二產(chǎn)業(yè)的進(jìn)步，第二三產(chǎn)業(yè)對(duì)第一產(chǎn)業(yè)有帶動(dòng)作用?？偟恼f(shuō)來(lái)，國(guó)內(nèi)生產(chǎn)總值反映了一個(gè)國(guó)家的生產(chǎn)規(guī)模，同時(shí)也反映了這個(gè)國(guó)家的產(chǎn)業(yè)結(jié)構(gòu)，因此，了解該國(guó)家的產(chǎn)業(yè)結(jié)構(gòu)現(xiàn)狀及其發(fā)展

12、變化規(guī)律，制定正確的產(chǎn)業(yè)發(fā)展政策，引導(dǎo)產(chǎn)業(yè)結(jié)構(gòu)健康地發(fā)展具有重要的意義。.-附錄： SAS 程序：.-3204.35251.63506.63831 6587.24510.1方差分析5017 7717.45813.55288.69102.272275800 11699.5 9138.66882.116428.5 11323.89457.222372.2 1493011993 28537.9 17947.213844.2 33612.9 20427.514211.2 37222.7 23028.714599.6 38691.8 26104.3cards; 342.9 141.8 194.33781

13、92.5 253.5392211.7 255.3421222.2 266.8443.9 280.7 303.4430317321445.9 483.5 377.6383.8 615.5 439.7340.7 648.2 468.1441.1 388.9 390453.1 359.3 336.9497.5 407.6 328.2559513.5 381.5651.1 602.2 462.8702.2 709.5 456.3714.2 602.8 456.9726.3 537.3 459.5736.2 689.1 512.6793.3 912.2 547.2826.3 1022.8577.3827

14、.4 1084.2606.5907.5 1173 640.4945.2 1192 652.7971.1 1370.5655.79671337.2639.5942.1 1509.1750.71018.41745.2860.51258.91913.5865.81359.42192 966.41545.62255.51061.31761.62383 1150.11960.82646.21327.52295.53105.71769.82541.63866.62556.22763.94492.72945.624040 1037208855528627 12583111135233702 14900315

15、7639.6 187383.222059220326052377 235319231626;title 雙方差分析;proc anova; class i g; m y= i g; run;相關(guān)系數(shù)分析四data gdp;input id cta firs seco thir ;.14770.02847 41033.5815933873.4446914944.7225 45555.8779638713.9538515781.26905 49512.2909744361.6105416537.01966 53896.7677949898.9018217381.7177 62436.3121456

16、004.7263121412.73404 73904.3118764561.2920123070.43756 87364.5791173432.866054228 7278 5403.2data gdp;=1952 to 2012;do g=1 to 3;input y;output;end;end;-19762943.79671337.2label id=年份 cta=國(guó)內(nèi)生產(chǎn)總值639.519773201.9942.1 1509.1 firs=第一產(chǎn)業(yè) seco=第二產(chǎn)750.719783624.11018.41745.2業(yè) thir=第三產(chǎn)業(yè); 860.5cards;19794038.2

17、1258.91913.5865.819804517.81359.42192966.419814862.41545.62255.51061.319825294.71761.623831150.119835934.51960.82646.21327.5198471712295.53105.71769.819858964.42541.63866.619631233.3497.5407.62556.2328.2198610202.22763.94492.72945.6198711962.53204.35251.63506.6198814928.338316587.24510.1198916909.24

18、22872785403.219681723.1726.3537.3199018547.950177717.45813.5459.5199121617.85288.69102.2722719691937.9736.2689.1199226638.1580011699.59138.6512.6199334634.46882.116428.519702252.7793.3912.211323.8547.2199446759.49457.222372.219712426.4826.31022.814930577.3199558478.111993 28537.919722518.1827.41084.

19、217947.2606.5199667884.613844.233612.919732720.9907.5117320427.5640.4199774462.614211.237222.719742789.9945.2119223028.7652.7199879395.714599.638691.819752997.3971.11370.526104.3655.7199989677.05475 14770.0284741033.58159 33873.44469.2 709.5 456.31967 1773.9714.2 602.8456.9513.5 381.51965 1716.1651.

20、1 602.2462.81950 679342.9 141.8 194.31951 824378192.5 253.51952 859392211.7 255.31953 910421222.2 266.8.9 280.7 303.4317321.9 483.5 377.6.8 615.5 439.7.7 648.2 468.1.1 388.9 3901962 1149.3453.1 359.3336.9-2000 99214.55431 14944.7225 45555.87796 38713.953852001 109655.1706 15781.26905 49512.29097 443

21、61.610542002 120332.6893 16537.01966 53896.76779 49898.901822003 135822.7561 17381.7177 62436.31214 56004.726312004 159878.3379 21412.73404 73904.31187 64561.292012005 183867.8827 23070.43756 87364.57911 73432.866052006 21631424040 103720885552007 26581028627 1258311113522008 31404533702/p>

22、02007 34090335226 1576391480382008 401512.8 40533.6 187383.2 1735962009 47156447712 2205922032602010 51932252377 235319231626;title 國(guó)內(nèi)生產(chǎn)總值相關(guān)系數(shù)分析 oc corr; var cta firs seco thir; run;回歸分析.-label id=年份 cta=國(guó)內(nèi)生產(chǎn)總值 seco=第二產(chǎn)業(yè) thir=第三產(chǎn)業(yè);cards;1950 679141.8 194.31951 824192.5 253.51952 859211.7 255.31953 9

23、10222.2 266.8.7 303.4321.5 377.6.5 439.7.2 468.1.9 3901962 1149.3359.3 336.91963 1233.3407.6 328.2.5 381.51965 1716.1602.2 462.8.5 456.31967 1773.9602.8 456.91968 1723.1537.3 459.51969 1937.9689.1 512.61970 2252.7912.2 547.21971 2426.41022.8577.31972 2518.11084.2606.51973 2720.91173 640.41974 2789.9

24、1192 652.71975 2997.31370.5655.71976 2943.71337.2639.51977 3201.91509.1750.71978 3624.11745.2860.51979 4038.21913.5865.81980 4517.82192 966.41981 4862.42255.51061.31982 5294.72383 1150.11983 5934.52646.21327.51984 7171 3105.71769.81985 8964.43866.62556.21986 10202.2 4492.72945.61987 11962.5 5251.635

25、06.61988 14928.3 6587.24510.11989 16909.2 7278 5403.2.data gdp;input ideco thir ;-1990 18547.9 7717.45813.51991 21617.8 9102.272271992 26638.1 11699.5 9138.61993 34634.4 16428.5 11323.81994 46759.4 22372.2 149301995 58478.1 28537.9 17947.21996 67884.6 33612.9 20427.51997 74462.6 37222.7 23028.71998

26、79395.7 38691.8 26104.31999 89677.05475 41033.58159 33873.444692000 99214.55431 45555.87796 38713.953852001 109655.1706 49512.29097 44361.610542002 120332.6893 53896.76779 49898.901822003 135822.7561 62436.31214 56004.726312004 159878.3379 73904.31187 64561.292012005 183867.8827 87364.57911 73432.86

27、6052004 216314103720885552005 2658101258311113522006 3140451490031313402007340903 1576391480382008 401512.8 187383.2 1735962009 4715642205922032602010 519322235319231626；title 國(guó)內(nèi)生產(chǎn)總值回歸分析分析;proc glm;m cta= seco thir;run;參考書：SAS 統(tǒng)計(jì)分析與應(yīng)用實(shí)例著.-.國(guó)內(nèi)生產(chǎn)總值數(shù)據(jù)本表按當(dāng)年價(jià)格計(jì)算：億元年 份國(guó)民生產(chǎn)總值第一產(chǎn)業(yè)第二產(chǎn)業(yè)第三產(chǎn)業(yè)年 份國(guó)民生產(chǎn)總值第一產(chǎn)業(yè)第二產(chǎn)業(yè)第三

28、產(chǎn)業(yè)1950679342.9141.8194.319835934.51960.82646.21327.51951824378192.5253.5198471712295.53105.71769.81952859392211.7255.319858964.42541.63866.62556.21953910421222.2266.8198610202.22763.94492.72945.6.9280.7303.4198711962.53204.35251.63506.6317321198814928.338316587.24510.1.9483.5377.6198916909.242287278

29、5403.2.8615.5439.7199018547.950177717.45813.5.7648.2468.1199121617.85288.69102.27227.1388.9390199226638.1580011699.59138.619601149.3453.1359.3336.9199334634.46882.116428.511323.819611233.3497.5407.6328.2199446759.49457.222372.214930513.5381.5199558478.11199328537.917947.219631716.1651.1602.2462.8199667884.613844.233612.920427.5.2709.5456.3199774462.614211.237222.7230