各地区生产总值及构成项目的差异

目录TOC\o"1-5"\h\z1实训的题目与数据1实训的题目与目的1样本与变量1数据文件12统计软件操作3对我国收入生产总值及构成项目差异的描述性统计分析3对我国收入生产总值及构成项目差异的统计软件处理3对我国收入生产总值及构成项目差异的计算结果的说明4对我国收入生产总值及构成项目差异的多元线性分析5对我国收入生产总值及构成项目差异的多元线性软件处理5对我国收入生产总值及构成项目差异的多元线性计算结果的说明10对我国收入生产总值及构成项目差异的聚类分析152.3.1对我国收入生产总值及构成项目差异的聚类分析的软件处理152.3.1对我国收入生产总值及构成项目差异的聚类分析的计算结果的说明193对我国收入生产总值及构成项目差异的研究结论23我国收入生产总值差距过大23劳动者报酬占收入生产总值比重最高23我国东中西部收入生产总值差距过大23参考文献241实训的题目与数据1.1实训的题目与目的(1)实训的题目：2013我国生产总值及构成项目的差异分析(2)实训的目的：了解我国各地区生产总值之间的差异以及构成项目之间的差异。样本与变量样本为我国31个省、自治区、直辖市(不包括港澳台)，变量为y1营业盈余、y2固定资产折旧、y3生产税净额、y4劳动者报酬、y5收入生产总值。数据文件依据2013年《中国统计年鉴》选取我国31个省、自治区、直辖市(不包括港澳台)的营业盈余、固定资产折旧、生产税净额、劳动者报酬、收入生产总值。原始数据见表1。表12013我国各地区生产总值及构成项目的差异地区二|yiy2ny3V4irys匚1北京市4034392368623179S410217.9619B00.E1।2n天津市4639301615592431555755.571444201「3河北省6255363671.744196901431B.9528442954n山西省2587.742180962133.085763.471266525r5内蒙古5283.83174B232254.6976297516916.506n辽宁省6621Q933081150854411690.6827213.22匚7吉林省3695.092192562139255019501304640s-1黑龙江省4626541598162172.246057,9614454.91匚9上海市55003525383944623493163721B1815L二I江苏省18166.28770471809673257856559753.3711浙江省9483.104663635711.101789875377565812口安徽省4874132476.463066.8387919219229.3413福建省5257.542239.7430936311277572186849匚14江西皆40303819278522408662111014410.1915山东省16269327654.539208.6922097.805523032|16河南省7797.263649644667.76161767432191.30

匚17湖北省6059.043064.343692.5312075.9224791.6313湖南省5615.8326497240461712309952462167匚19广东省15277.9477S3.879623.3929809.60624747920广西3033971538492857456920001444990匚21海南省523.04469.31427.751767.463177.5622重庆市3809.1215397822B1.3651530012783.26匚23四川省7640.853252.594149.96113466726392.07匚24贯州省1270631016271473D1432695808686匚25云南省2153.021255.362430.225993.7111832.31匚26西藏112.651166375.86510.5381567匚27陕西省4S19.712014.4S2912.5567587416205.4523甘肃省1131.491032.4798772312901633069匚29青海省474.92410.46279.60957.082122.0630宁夏516.8840360257757匚31新疆1442.821223.201197.1845S0648443.64IJ金资料来源：2013年《中国统计年鉴》2统计软件操作对我国收入生产总值及构成项目差异的描述性统计分析对我国收入生产总值及构成项目差异的统计软件处理调出主菜单(Analyze/DescriptiveStatistics/Descriptives)首先打开要进行描述统计分析的数据文件。其次在SPSS统计软件数据数据编辑窗口，依次点击Analyze/DescriptiveStatistics/Descriptives,调出描述统计主菜单并进入描述统计主对话框Descriptive(见图1)。(2)主对话框的设置在描述统计主对话框中，从左侧的源变量框中将研究变量y1、y2、y3、y4、y5等调入到右侧的Variable(s)框中。同时选择对话框左下角的Savestandardizedvaluesasvariables,输出各个变量的标准化值(见图2)。(3)二级对话框的设置(Optins)在主对话框中，点击Options按钮，进入二级对话框Descriptives:Optins。在二级对话框中，根据研究目的选择需要计算的指标，如均值Mean标准差Std.deviation、最小值Minimum最大值Maximum偏度Skewness峰度Kurtosis等(见图3)。上述设置完成后，点击Continue返回主对话框。(4)运行程序DESCRIPTIVESV^RlABLES=yly2y3y4y5/STATISTICS=MEANSTDDEVMINMAXKURTOSISSKEWNESS.FitEM'力DataTrans^nn邛面早MutliKMmm申力底5r.因1a胤H雪国威应醛49♦嗨堂区个法如产特3E1女■知租1jLBTil4D34393邓电31791M1U1T1H199W312天津市4EH-3Q圮243155ST555714442013耳北詈EZ5B961R1T44IW9CM3IB352M42骷❽:含口d山B91网TH2谢q。用3MSN3471第能2555WW1WI23弊1历脸94嫡16钟--*6口的6£对旧JE0B11丐值54411G3B5E272132271ES5M211ZKM3&2&印博州1S44E4Dnxna4GH1的*晡21空勤锄1船UJ643I零上即明匾消2例铸血A2婚»1$372恻8隹1b惘1跳宛77MT1ibHTl因拈玷SS753J71t近工*34AJ104663EJ17BKTB3TT&E5E•it"j以881ZttBW^^.1324?EW1国5助ITO3Z19229M13iWi1年57启2251Ta7姬总H2?T87曾獭啪1i・明算1K7IB22成鞋我11M«J1Q1915山开安7时51被郎弱切的整55230J21Ei同博*7用735JMB&I4G6T7&tflfG743ZW3DIT琳>的M1K4M5出螳24F»116■Mi描旺园2除?豫mt?1?W9952崎IS1Bf¥1翳77则7763B?46233*<534747520ra3HJB71理将理丁码E82D0D14445MZ1鼻即■5Z3M45931WTT51T5J4E3«FT5Eas■吒击3Kti2163ITi豁1然死3岫1271326二，一……frleu±PW占gidii^图1描述统计菜单

图2描述统计主对话框的设置埼Descriptives；OptionsMMean□SumDispersion“Std.deviationWMinimum匚VarianceVMaximumCLRange|__SEmeanDistribution“Kurtosis4SKewnessDisplayOrder©^VariablelistOAlphabetic<_：Ascendingmeans。DescendingmeansContinueCancelHeJp图3描述统计二级对话框的设置2.1.2对我国收入生产总值及构成项目差异的计算结果的说明运行以上程序，即可得到描述统计分析的计算结果，并据此给出研究结论

表2描述统计量N极小值极大值均值标准差营业盈余31112.6518166.285251.72944448.86106固定资产折旧31116.637763.872571.90772015.62725生产税净额3175.869623.383253.84132386.54399劳动者报酬31510.5329809.609385.27427055.82470收入生产总值31815.6762474.7920462.752315709.73731有效的N(列表状态)31表2的计算结果是我国31个省、自治区、直辖市的收入生产总值及构成项目的平均水平、差异程度和分布形态。在收入生产总值的均值中，劳动者报酬的均值最高，达到了9385.2742元；固定资产折旧最低，仅为2571.9077元。离散程度最大的是收入生产总值，最小的是固定资产折旧。2.2对我国收入生产总值及构成项目差异的多元线性分析对我国收入生产总值及构成项目差异的多元线性软件处理(1)调出主菜单(Analyze/Regression/Linear)在数据编辑窗口，依次点击Analyze/Regression/Linear,调出多元线性回归分析菜单并进入多元线性回归分析主对话框LinearRegression(见图4)

Fil5-EdilView0^19Trgn§fqtfmAnalyzeDirectMarketingGnpjphs旦tilrt]/Aqld-qn&WindowHelp西『月即口弟►口电白仃巾如自S1atigtiC£►Tables，CompareMeanLSkGensiralLinearModel►1地区yiy5varvar恸1北京4034.399619800.Bl2天津4639.30GenefallzedLinearbedelsk5714442.DI——.河北6255.36MixedModels►9628442.064ill四2587.74Cotrelate卜IT136BS.265内蒙52B383RegressiCinkLoglln#ar，NewiNetworks►Classify-►DlmeniglcinRedudson,ScalekNafiparamEtricTssts►r6辽宁662109gUnear.7吉林369509E黑龙463654IZiCurveEsilimartJQn..9上海560035明ParllalLeastSquares...10江苏1616628QEHnaiyLogistic...11浙江546310।MuHihDmialLogistic...E3。整nalSftoWt12安演48741313福建5257.S4Forecas.tingkSurvivalr14江西4030.3B16山东16M9.32MultipleRemp口nsef由MissingValus.Analysis...MultipleImpul?tion►ComplexSaim制自专►QualityCoolrol>囱ROCCuvve^EigNonlinear..展世制ghtEstimallon...点2-SlageLeas-tSquares..16间隔7797.2E17湖北6059.0418湖喃5615.83OplimalScaling(CATREQ)...19f15277.94bU—20广西3083.97)014449.9021海南62304Jlutjf4b3177.5622里氏3MJ91215397B22B1356153.OQ1Z7032623四川76408532525941499611348672639207IlkD^taViewVariableLin网...图4多元线性回归分析菜单（2）主对话框的设置（Linear）设置因变量。在主对话框中，将因变量y5调入到主对话框右上角的因变量Dependent框中。设置自变量。设置完因变量后，在主对话框中，依次选择自变量y1,y2,y3,y4并调入到主对话框右侧中部的自变量框Independent（s）中。设置自变量进入模型的方法。在主对话框的右侧中部的Method框中,按照默认选项选择Enter（见图5）。（3）各个二级对话框的设置①设置模型参数及有效性的检验1。在主对话框中，点击Statistics,进入多元线性回归分析二级对话框LinearRegression:Statistics。在二级对话框中，主要进行以下设置。设置模型参数和模型参数有效性检验。在RegressionCoefficients框中选择Estimates,输出模型参数及模型参数有效性检验的结果（即t检验）。设置模型的拟合优度检验和模型有效性检验。选择右侧的Modelfit,输出模型的拟合优度检验和模型有效性检验的结果（即R检验与F检验）。

图5多元线性回归分析主对话框设置多重共线性检验。选择右侧的Collinearitydiagnostics,输出多重共线性检验的结果（即VIF检验）。设置序列自相关检验。在左侧的Residuals框中，选才DDurbin-Watson,输出序列自相关检验结果（即DW佥验）。设置异常值与强影响点检验。在左侧的Reaiduals框中，选择Casewisediagnostics,输出异常值与强影响点检验的结果。上述内容见图6。完成后，点击Continue返回主对话框。

图6多元线性回归分析模型参数与有效性检验设置模型参数有效性检验2。在主对话框中，点击Plots,进入多元线性回归分析二级对话框LinearRegression:Plots二级对话框。在二级对话框中，主要进行以下设置。正态性检验。在下方的StandardizedResidualPlots框中，选择Normalprobabilityplot,进行正态性检验。方差齐检验。在左侧框中，选择ZPRED故X轴，选择ZRESID故Y轴，进行方差齐检验。之后返回上级对话框（见图7）。图7多元线性回归分析模拟有效性检验设置输出内容与模型检验3。在主对话框中，点击Save进入LinearRegression:Save二级对话框。在这个对话框中，主要有以下设置。输出点预测值。在PredictedValues框中，选择Unstandardized,输出非标准化的点预测值。输出预测区间。在PredictedValues框中，选择Individual,输出概率为95%勺预测区间。因变量y的异常值和强影响点检验。在Residuals框中分别选择标准化残差Standardized、学生化残差Studentized、剔除残差Deleted,进行因变量y的异常值和强影响点的检验。自变量x的异常值和强影响点检验。在Distances框中，选择Cook's库克距离，进行自变量x的异常值和强影响点的检验（见图8）。上述设置完成后，返回到主对话框。④设置自变量进入模型的条件。在主对话框中，点击Options进入LinearRegressionOption二级对话框，查看自变量进入模型或被剔除模型的条件(也可修改软件指定的条件)。(4)运行程序(OKREGRESSION/MISSINGLISTWISE/STATISTICSCOEFFOUTSRANOVACOLLINTOL心RITERIA=PIN(.05)POUT(.10)CIN(95)/NOORIGIN/DEPENDENTy5/METHOD=ENTERy1y2y4/SCATTERPLOT=(*ZRESID,*ZPRED)/RESIDUALSDURBIN/CASEWISEPLOT(ZRESID)OUTLIERS(3)/SAVEPREDCOOKICINSRESIDDRESID

图8多元线性回归分析模型输出项设置2.2.2对我国收入生产总值及构成项目差异的多元线性计算结果的说明(1)模型的分析建立运用SPS筑计软件，计算得到多元线性回归分析的模型参数，见表3表3系数a模型非标准化系数标准系数tSig.共线性统计量B标准误差试用版容差VIF1（常量）营业盈余207.2311.034145.225.094.2931.42710.965.165.000.04224.007固定资产折旧1.694.253.2176.682.000.02835.639劳动者报酬1.116.049.50122.666.000.06116.469a.因变量：收入生产总值根据计算结果，可以得到多元线性回归模型:lny5=207.231+1.034y1+1.694y2+1.116y4该模型反映了三个构成项目对收入生产总值的影响方向和影响程度，表明三个构成项目的变动对收入生产总值的影响。根据模型，我国2013年三个构成项目影响收入生产总值的弹性系数分别为1.034,1.694和1.116。说明营业盈余和劳动者报酬的增加对收入生产总值增加的影响没有固定资产折旧的增加对收入生产总值增加的影响大。固定资产折旧每增加1个百分点，收入生产总值九增加1.694个百分点。说明目前情况下增加固定资产折旧收入是增加收入生产总值的最佳途径。（2）模型有效性检验拟合优度检验（R检验）。拟合优度检验专门用来研究模型自变量与因变量之间的拟合程度。表4模型汇总b模型RR方调整R方标准估计的误差Durbin-Watson1a1.000.999.999468.774862.599a.预测变量：（常量），劳动者报酬I,营业盈余，固定资产折旧b.因变量：收入生产总值根据检验的计算结果，R,R方和调整咽分别为1.0,0.999和0.999,都属于高度相关，所以该多元线性回归模型的拟合优度较高，自变量对因变量的解释程度较高，模型显著有效。序列自相关检验（D惭验）。根据上面的拟合优度计算表，可以看到模型序列自相关的检验结果。D怖计量为2.599。根据序列自相关检验理论，D惭计量离DWJ中心值较远，所以接受零假设，即基本认定模型存在序列自相关。模型有效性检验（F检验）。模型有效性检验是从整体上检验模型的有效性，这种检验对于多元线性回归模型非常重要。检验方法一般采用F检验c检验数据见表5。表5Anovaa模型平方和df均方FSig.回归7397942140.7632465980713.5811221.762A”b.000581残差5933246.45527219749.869总计7403875387.22300a.因变量：收入生产总值b.预测变量：（常量），劳动者报酬I,营业盈余，固定资产折旧由于相伴，概率Sig小于显著性水平，所以拒绝零假设而接受备选假设即模型所有参数不同时为零，认为模型整体有效。多重根据共线性检验（VIF检验）。根据上面的模型建立的计算表（表3）,可以看到模型多重共线性的检验结果。根据容忍度Tol的计算结果，三个构成项目的Tol值分别是0.042,0.028,0.061,均接近普通临界值0,所以可以认定该模型存在严重的多重共线性。若根据方差膨胀系数VIF的计算结果，三个产业的VIF值分别是24.007,35.639,16.469均大于普通临界值10,也可认为该模型存在严重的多重共线性。

回归标准化残基的标潴P-P图回归标准化残基的标潴P-P图图9P-P概率图期望的累枳概率正态性检验。根据P-P概率图，可以看到残差的实际累计概率与理论累计概率之间比较接近，可以认为残差的正态性假设基本成立，可以对此进行多元线性回归分析

3-2-1-0-散点图因变量：收入生产总值回归标潴化残差-3111ttr3-2-1-0-散点图因变量：收入生产总值回归标潴化残差*2*10123回归标淮化预计值图10残差散点图⑥方差齐检验。根据图10所示的散点图，可以看到所有的散点都落在横轴、纵轴土3的范围之内，可以认为残差的方差齐假设基本成立，可以对此进行多元线性回归分析。综上所述，多元线性回归分析模型有效性检验的基本结果是：模型拟合优度较高,模拟整体有效，模型的多重共线性严重，模型的正态性和方差齐检验也都通过。所以，整体模型有效，可以使用。2.3对我国收入生产总值及构成项目差异的聚类分析2.3.1对我国收入生产总值及构成项目差异的聚类分析的软件处理(1)调出主菜单(Analyze/Classify/K-MeansCluster)在数据编辑窗口，依次点击Analyze/Classify/K-MeansCluster,调出快速聚类分析菜单并进入快速聚类分析主对话框K-MeansClusterAnalysis(见图11)。"ilgTranstorrn夕峥依金M^rke'lingGraphsW加/ftcM-onswin^W宅HRegoris►DescnoliveStam别les，Tables►CompareMeanskGah&ralLinearModel►GenieralizeULmearfMo-Clelak■唯;期A局SRE_1UICI_1PRE]芟量变量13>6471.6323JB.4759221S684930721S68.4B45014)0330818201S14814410197901441018969151701-2.9994?MixedModels►ggifir日Ms►,艮电第自力自2n►IL^glmeair►NemralNelwoirks►D.318125523033790S5230.32B01t6M481.29M032191.31237321911.3036917)154H82316即用13397Q2479183H3】1羽-,4432JM战湖246215801319,ai66fr4.77M662J74.7970862J74.767fi7201066-2.5703(r.najirm-hjjJidn414449.90939c。”的►Dlmens^nReduction►Sc^le►Nonparametricmisk部1Two^epCiustef21)210-5157431775619222»81137W囹Clu&i€r施|HierardnicsiClu划E「.127B3.2492623)121.2987(26392.0688824)088^.21848同小.&08666032鸿W671Fdrecasting►MultipheResponse►目MissingvaimeMal号吗,，Mu可口陋Imputallon*CamplesSamples.►QualityComlrol►MDiscriminanl...NflamstNeighlHir...11B32梵白9826)240-.5845：S15.6721627)179.44211162Q6.448M―-,~1V.J，M4wrWTSTTr28)047-1136(06820263306988263306904229)144-.350632052992122W96921220612930)1711-.4194-7.663202S77.579912577.6716631)163-.393633833013443849鸵844384144320rocCurve...3334■产HI图11K-Means聚类分析分析菜单(2)主对话框的设置在K-Means聚类分析主对话框中，分别进行以下设置：①聚类分析变量的设置在主对话框中的左侧，分别选择5个变量“y1至y5"，调入到右侧的Variables框中。②聚类对象(标签)的设置选择字符型变量“地区”，调入到LabelCasesby框中。③设置聚类数在NumberofClusters框中指定分类数，此处设置为5(这个过程非常重要，一般要经过多次的尝试才可完成，一定要十分慎重)。

④聚类类中心的设置聚类类中心的设置有2种方式：查找已有的类中心点击左侧的Centers,可以调出查找已有类中心的对话框，通过查找已有的SPSS数据文件（类中心）作为本次聚类分析的类中心。其前提是以前做过类似的聚类分析而有类中心；以前有过类似的聚类分析而又类中心。计算指定类中心若以前没有做过本问题的聚类分析，也没有过类似的聚类分析，则由计算机（SPSSB计软件）代替用户而自动指定类中心。此时，用户不用进行人户操作。⑤设置聚类分析的方式在Method框下选中Iterateandclassify,表示以初始类中心开始，经过一定的迭代过程，直到得到最后的类中心，聚类分析完毕（见图12）0图12K-Means聚类分析主对话框的设置（3）各个二级对话框的设计设置聚类分析终止的条件在主对话框中，点击Iterate按钮，进入K-MeansClusterAnalysis:Iterate二级对话框。按照默认的选项进行，用来设置程序终止的条件。终止的条件有两个：计算次数系统默认为10次，当计算达到10次时，则聚类分析自动终止。但实际山，计算的次数在2-5次左右。移动距离系统默认的类中心的移动距离为0.02（SPSSa计软件上看到的是0.000）若移动距离达不到0.02,则聚类分析自动终止。这两个条件达到一个，聚类就结束（见图13）。上述设置完成后，返回主对话框。图13K-Means聚类分析终止条件的设置设置保存内容在主对话框中，点击Save按钮，进入K-MeansCluster:SaveNewVariables二级对话框。选择2个备选项ClusterMembership和Distancefromclustercenter,即保存聚类分析的类号和该个体距离其所属的类的类中心的距离（见图14）。上述设置完成后，返回主对话框。十@K-MeansClusterSaveNe...囹Clustermembershipinmaihi■■■»■■inraiaimamrainhi■imamrairiraiaima^“^Distancefromclustw8nt4rfContinueCancelHelp图14K-Means聚类分析保存内容的设置③检验聚类分析变量的有效性在主对话框中，点击Options按钮，进入K-MeansClusterAnalysis:Options二级对话框Statistics框选中3个选项Initialclustercenters,即输出初始类中心，ANOVAable以聚类分析产生的类为控制变量，对参与聚类分析的所有变量进行的单因素方差分析，目的是判断该变量是否有“资格”作为聚类分析的依据。Clusterinformationforeachcase表示输出个案聚类信息及其距离所属类中心的距离MissingValues框选才¥Exclues,表示去除所有含缺失值的个案后再进行聚类分析（见图15）。上述设置完成后，返回主对话框图15K-Means聚类分析保存内容的设置(4)运行程序QUICKCLUSTERy1y2y3y4y5/MISSING=LISTWISE/CRITERIA=CLUSTER(5)MXITER(10)CONVERGE(0)/METHOD=KMEANS(NOUPDATE)/SAVECLUSTERDISTANCE/PRINTID(地区)INITIALANOVACLUSTERDISTAN.

2.3.1对我国收入生产总值及构成项目差异的聚类分析的计算结果的说明表6初始聚类中心聚类12345营业盈余3809.127640.8515277.94112.659483.10固定资产折旧1539.783252.597763.87116.634663.63生产税净额2281.364149.969623.3875.865711.10劳动者报酬5153.0011348.6729809.60510.5317898.75收入生产总值12783.2626392.0762474.79815.6737756.58由于我们以前没有做过类似的研究，所以缺少用来做聚类分析的类中心，故初始类中心是由计算机自动收索的类中心。在这份数据中，很明显是第3类的类中心的数据最优，第5类、第2类、第1类直到第4类依次降低。表7迭代历史记录a迭代聚类中心内的更改1234511281.3872718.6185332.2522508.4293126.8782122.873668.006.000993.429.0003589.730.000.000770.980.0004.000.000.000.000.000a.由于聚类中心内没有改动或改动较小而达到收敛。任何中心的最大绝对坐标更改为.000。当前迭代为4。初始中心间的最小距离为13411.850。它反映了5个类中心每次迭代是聚类偏移的过程。经过第1次的迭代，5个类的类中心分别偏移了1281.387、2718.618、5332.252、2508.429、3126.878,经过第2次的迭代，5个类的类中心分别偏移了122.873、668.006、.000、993.429、.000,经过第3次的迭代，5个类的类中心分别偏移了589.730、000、000、770.980、.000,到第4次计算的时候，位置不再移动。所以到第4次迭代时，5个类的类中心（几乎）没有发生变化，所以聚类分析到此结束。表8最终聚类中心聚类12345营业盈余3842.875762.0616571.18788.928640.18

固定资产折旧1766.142896.667707.70665.934106.59生产税净额2385.333877.138976.27687.725189.43劳动者报酬6126.2811261.7725897.682365.1817037.75收入生产总值14120.6223797.6159152.834507.7534973.94与最初类中心类似，表中数据显示的是聚类分析的最终类中心。数据显示，第3类各项指标的数值最好，第5类次之，第2类、第1类依次降低，第4类最差表9ANOVA聚类误差FSig.均方df均方df营业盈余142240667.2114954164.30626149.074.000固定资产折旧29177650.6414198922.83626146.678.000生产税净额40716482.9844307762.87126132.298.000劳动者报酬35449322835226121.969.000收入生产总值1799041953.41747988752.82926225.197.000F检验应仅用于描述性目的，因为选中的聚类将被用来最大化不同聚类中的案例间的差别。观测到的显著性水平并未据此进行更正，因此无法将其解释为是对聚类均值相等这一假设的检验。表中数据反映各个变量在不同类的均值的比较情况。5个变量的F统计量都较大，分别是149.074、146.678、132.298、121.969、225.197,同时其对应的相伴概率sig值都很小，分别是.000、.000、.000、.000、.000,都小于显著性水平0.05,故拒绝零假设，认为用于聚类分析的5个变量对于本次K-mea课类分析都是有效的，可以据此进行聚类分析。表10聚类成员案例号地区聚类距离1北京市24562.2222天津市1948.5903河北省25645.5184山西省12015.1375内蒙古13488.7156辽宁省23858.220

7吉林省11625.7708黑龙江省1896.7619上海市22870.99510江苏省31921.23411浙江省