应用回归分析第九章部分答案

1、第9章非线性回归9.1在非线性回归线性化时，对因变量作变换应注意什么问题？答：在对非线性回归模型线性化时，对因变量作变换时不仅要注意回归函数的形式，还要注意误差项的形式。如：（1）乘性误差项，模型形式为二AKL-e:加性误差项，模型形式为八AKL；对乘法误差项模型（1）可通过两边取对数转化成线性模型，（2）不能线性化。一般总是假定非线性模型误差项的形式就是能够使回归模型线性化的形式，为了方便通常省去误差项，仅考虑回归函数的形式。9.2为了研究生产率与废料率之间的关系，记录了如表9.14所示的数据，请画出散点图，根据散点图的趋势拟合适当的回归模型。表9.14生产率x（单位/周）100020003

2、0003500400045005000废品率y（%）5.26.56.88.110.210.313.0解:先画出散点图如下图:12.0010.008.006.00从散点图大致可以判断出x和y之间呈抛物线或指数曲线，由此采用二次方程式和指数函数进行曲线回归(1)二次曲线SPSS输出结果如下：ModelSummaryTheindependentvariableisx.RRSquareAdjustedRSquareStd.ErroroftheEstimate.981.962.942.651ANOVATheindependentvariableisx.SumofSquaresdfMeanSquareFS

3、ig.Regression42.571221.28650.160.001Residual1.6974.424Total44.2696CoefficientsUnstandardizedCoefficientsStandardizedCoefficientstSig.BStd.ErrorBetaX-.001.001-.449-.891.423x*24.47E-007.0001.4172.812.048(Constant)5.8431.3244.414.012从上表可以得到回归方程为："5.843-0.087x4.4710x2由X的系数检验P值大于0.05,得到X的系数未通过显著性检验。

4、由X2的系数检验P值小于0.05,得到X2的系数通过了显著性检验。(2)指数曲线ModelSummaryRRSquareAdjustedRSquareStd.ErroroftheEstimate.970.941.929.085Theindependentvariableisx.ANOVASumofSquaresdfMeanSquareFSig.Regression.5731.57379.538.000Residual.0365.007Total.6096Theindependentvariableisx.CoefficientsThedependentvariableisln(y).Unsta

5、ndardizedCoefficientsStandardizedCoefficientstSig.BStd.ErrorBetax.000.000.9708.918.000(Constant)4.003.34811.514.000从上表可以得到回归方程为：4.OO3e°.0002t由参数检验P值0<0.05,得到回归方程的参数都非常显著。OObseiveduadiatLcEspaneiita.从R2值，°的估计值和模型检验统计量F值、值及拟合图综合考虑,指数拟合效果更好一些。9.3已知变量x与y的样本数据如表9.15,画出散点图，试用ae®/x来拟合回归模型

6、，假设：表9.15(1) 乘性误差项，模型形式为y=ae卩/xe£加性误差项，模型形式为y=ae/x+£。序号xy序号xy序号xy14.200.08663.200.150112.200.35024.060.09073.000.170122.000.44033.800.10082.800.190131.800.62043.600.12092.600.220141.600.94053.400.130102.400.240151.401.620解：散点图:0.500-0.000-1.001.002.003.004.00(1)乘性误差项，模型形式为y=ae卩/xe£线性化

7、：lny=lna+®/x+e令yl=lny,a=lna,x1=1/x做y1与x1的线性回归，SPSS输出结果如下:ModelSummarPredictors:(Constant),x1ModelRRSquareAdjustedRSquareStd.ErroroftheEstimate1.999a.997.997.04783a. DependentVariable:y1ANOVAbModelSumofSquaresdfMeanSquareFSig.1Regression10.930110.9304778.305.000aResidual.03013.002Total10.96014aP

8、redictors:(Constant),x1b.DependentVariable:y1CoefficientsaModelUnstandardizedCoefficientsStandardizedCoefficientstSig.BStd.ErrorBeta1(Constant)-3.856.037-103.830.000x16.080.088.99969.125.000aDependentVariable:y1从以上结果可以得到回归方程为：y仁-3.856+6.08x1F检验和t检验的P值0<0.05,得到回归方程及其参数都非常显著。回代为原方程为：y=0021e6-08/x(2

9、)加性误差项，模型形式为y=a&/x+£不能线性化，直接非线性拟合。给初值a=0.021,B=6.08(线性化结果)，NLS结果如下：ParameterEstimatesParameterEstimateStd.Error95%ConfideneeIntervalLowerBoundUpperBounda.021.001.020.023b6.061.0445.9656.157IterationHistoryIterationNumber3ResidualSumofSquaresParameterab10001.0216.0801.10010216.0612.0.001.021

10、5.Q6121.001.0216.06130.001.0216.06131001.0216.061ANOVASourceSumofSquaresdfMeanSquaresRegression4.45822.229Residual.00113.000UncorrectedTotal4.45915CorrectedTotal2.46714Dependentvariable:ya.Rsquared=1-(ResidualSumofSquares)/(CorrectedSumofSquares)=1.000.从以上结果可以得到回归方程为：y=0021占°61%根据R21,参数的区间估计不包括

11、零点且较短，可知回归方程拟合非常好，且其参数都显著。PredictedValues9.4Logistic回归函数常用于拟合某种消费品的拥有率，表&仃（书上239页，此处略）是北京市每百户家庭平均拥有的照相机数，试针对以下两种情况拟合Logistic回归函数。1y=刁-boHu（1）已知u=100，用线性化方法拟合，（2）u未知，用非线性最小二乘法拟合。根据经济学的意义知道，u是拥有率的上限，初值可取100；b0>0,0vblvl初值请读者自己选择。解：（1），u=100时，的线性拟合。对y=函数线性化得到：uAAAAAAIn（）=1.851-0.264In（）=1nb。tlnbi

12、，令yln（），作yy100y100y100关于t的线性回归分析，SPSS输出结果如下：ModelSummarjbModelRRSquareAdjustedRSquareStd.ErroroftheEstimate1.994a.988.987.16820a. Predictors:(Constant),tb. DependentVariable:y3ANOVAbModelSumofSquaresdfMeanSquareFSig.1Regression39.839139.8391408.165.000aResidual.48117.028Total40.32018a. Predictors:(C

13、onstant),tb. DependentVariable:y3CoefficientsaModelUnstandardizedCoefficientsStandardizedCoefficientstSig.BStd.ErrorBeta1(Constant)-1.851.080-23.039.000t-.264.007-.994-37.526.000a.DependentVariable:y3由表ModelSummary得到，R=0.994趋于1,回归方程的拟合优度好，由表ANOVA得到回归方程显著，由Coefficients表得到，回归系数都是显著的，得到方程：ln(-)=-1.850.

14、264，进一步计算得到：y100b0=0.157，0=0.768(u=100)回代变量得到最终方程形式为：y?1t0.01+最后看拟合效果，通过sequenee画图：由图可知回归效果比较令人满意。(2)非线性最小二乘拟合，取初值u=100,b0=0.157,13!=0.768:一共循环迭代8次，得到回归分析结果为：ParameterEstimatesParameterEstimateStd.Error95%ConfidenceIntervalLowerBoundUpperBoundu91.0622.03586.74795.377b.211.028.152.271e.727.012.701.75

15、3ANOVAaSourceSumofSquaresdfMeanSquaresRegression60774.331320258.110Residual85.369165.336UncorrectedTotal60859.70019CorrectedTotal15690.38618Dependentvariable:ya.Rsquared=1-(ResidualSumofSquares)/(CorrectedSumofSquares)=.995.u=91.062，R=0.995>0.994,得到回归效果比线性拟合要好，且:b0=0.211，D=0.727，回归方程为：y=11一-0.211

16、*0.727t91.062最后看拟合效果，由sequenee画图:PredictedValues10D口一得到回归效果很好，而且较优于线性回归9.5表9.仃（书上233页，此处略）数据中GDP和投资额K都是用定基居民消费价格指数（CPI）缩减后的，以佃78年的价格指数为100。（1）用线性化乘性误差项模型拟合C-D生产函数；（2）用非线性最小二乘拟合加性误差项模型的C-D生产函数；（3）对线性化检验自相关，如果存在自相关则用自回归方法改进；（4）对线性化检验多重共线性，如果存在多重共线性则用岭回归方法改进；解：（1）对乘法误差项模型可通过两边取对数转化成线性模型。Iny=lnA+:lnK+In

17、L令y'=lny,Bo=lnA,xi=lnK,X2=lnL,则传化为线性冋归方秫:y'=B°+二xi+：L+；SPSS输出结果如下：模型综述表ModelSummara.Predictors:(Constant),lnL,lnKbDependentVariable:lnYModelRRSquareAdjustedRSquareStd.ErroroftheEstimate1.997a.994.993.04836从模型综述表中可以看到，调整后的R为0.993,说明C-D生产函数拟合效果很好，也说明GDP的增长是一个指数模型。方差分析表ANOVAbPredictors:(Co

18、nstant),lnL,lnKModelSumofSquaresdfMeanSquareFSig.1Regression8.44624.2231805.601.000aResidual.05122.002Total8.49724a. DependentVariable:lnY从方差分析表中可以看到，F值很大，P值为零，说明模型通过了检验，这与上述分析结果一致。系数表Coefficientsaa.DependentVariable:InYModelUnstandardizedCoefficientsStandardizedCoefficientstSig.BStd.ErrorBeta1(Cons

19、tant)-1.7851.438-1.241.228lnK.801.056.86114.370.000lnL.402.171.1412.354.028根据系数表显示，回归方程为：111Y=IL785+O.8O1111K+0,4021nL尽管模型通过了检验，但是也可以看到，常数项没有通过检验，但在这个模型里，当InK和InL都为零时，InY为-1.785，即当K和L都为1时，GDP为0.168,也就是说当投入资本和劳动力都为1个单位时，GDP将增加0.168个单位，这种解释在我们的承受范围内，可以认为模型可以用。最终方程结果为：y=0618K0'801L0-404（2）用非线性最小二乘法

20、拟合加性误差项模型的C-D生产函数；上述假设误差是乘性的，现假设误差是加性的情况下使用非线性最小二乘法估计。初值采用（1）中参数的结果，SPSS输出结果如下：参数估计表ParameterEstimatesParameterEstimateStd.Error95%ConfidenceIntervalLowerBoundUpperBoundP.407.885-1.4292.243a.868.066.7311.006b.270.243-.234.774SPSS经过多步迭代，最终得到的稳定参数值为P=0.407，a=0.868,b=0.270y=0407K0'868L0-270为了比较这两个方

21、程，我们观察下面两个图线性回归估计拟合曲线图DPnGDPDPnGDP非线性最小二乘估计拟合曲线图非线性最小二乘估计拟合曲线图DPPredictedValuesDPPredictedValues我们知道，乘性误差相当于是异方差的，做了对数变换后，乘性误差转为加性误差，这种情况下认为方差是相等的，那么第一种情况(对数变换线性化)就大大低估了GDP数值大的项，因此，它对GDP前期拟合的很好，而在后期偏差就变大了，同时也会受到自变量之间的自相关和多重共线性的综合影响；非线性最小二乘法完全依赖数据，如果自变量之间存在比较严重的异方差、自相关以及多重共线性，将对拟合结果造成很大的影响。因此，不排除异方差、

22、自相关以及多重共线性的存在。(3) 对线性化回归模型采用DW检验自相关，结果如下:模型综述表ModelSummarfModelRRSquareAdjustedRSquareStd.ErroroftheEstimateDurbin-Watson1.997a.994.993.04836.715a. Predictors:(Constant),InL,InKDependentVariable:lnYDW=0.715<1.27,落在自相关的区间，所以采用迭代法改进DW=11=064252Xf=Xt-pXt将得到的数据再取对数，而后用普通最小二乘法估计，保留DW值模型综述表ModelSummarM

23、odelRRSquareAdjustedRSquareStd.ErroroftheEstimateDurbin-Watson1.983a.967.964478.902711.618Predictors:(Constant),Ltt,Ktta. DependentVariable:Ytt方差分析表ANOVAbModelSumofSquaresdfMeanSquareFSig.1Regression7.55423.777601.286.000aResidual.13221.006Total7.68623Predictors:(Constant),InLtt,InKtta. DependentVar

24、iable:lnYtt系数表Coefficientsaa.DependentVariable:lnYttModelUnstandardizedCoefficientsStandardizedCoefficientstSig.BStd.ErrorBeta1(Constant)-1.8591.470-1.265.220lnKtt.755.054.85214.098.000lnLtt.465.180.1562.577.018从模型综述表中可以看到，DW=1.618>1.45,认为消除了自相关；方差分析表中可以看到F值很大，P值为零，说明模型通过了检验。从系数表可得回归方程：】口丫；=-1,65

25、9+OT7551nKt+OT4651nL't再迭代回去，最终得方程为：LnytLnyti=-1.859+0.755(LnKtLnKti)+0.465(LnLtLnLti)(4) 对线性化回归方程通过VIF检验多重共线性:方差分析表ANOVAbModelSumofSquaresdfMeanSquareFSig.1Regression8.44624.2231805.601.000aResidual.05122.002Total8.49724a.Predictors:(Constant),lnL,lnKbDependentVariable:lnY系数表CoefficientsaaDepend

26、entVariable:InYModelUnstandardizedCoefficientsStandardizedCoefficientstSig.CollinearityStatisticsBStd.ErrorBetaToleranceVIF1(Constant)-1.7851.438-1.241.228lnK.801.056.86114.370.000.07713.034lnL.402.171.1412.354.028.07713.034多重共线性诊断表CollinearityDiagnosticsaModelDimensionEigenvalueConditionIndexVarianceProportions(Constant)lnKlnL11232.997.0031.63E-0051.00030.539429.012.00.001.00.00.09.91.00.001.00a.DependentVariable:InY直观法：从模型综述表上可以看到，F值很大，而t值很小，这是多重共线性造成的影响；VIF检验法：从系数表上可以看到，VIF=13&