第四版应用回归分析课后习题第八章

  • 格式:docx
  • 大小:240.69 KB
  • 文档页数:17

第8章 非线性回归 思考与练习参考答案 8.1 在非线性回归线性化时,对因变量作变换应注意什么问题? 答:在对非线性回归模型线性化时,对因变量作变换时不仅要注意回归函数的形式, 还要注意误差项的形式。如:

(1) 乘性误差项,模型形式为, (2) 加性误差项,模型形式为。 对乘法误差项模型(1)可通过两边取对数转化成线性模型,(2)不能线性化。 一般总是假定非线性模型误差项的形式就是能够使回归模型线性化的形式,为了方便通常省去误差项,仅考虑回归函数的形式。 8.2为了研究生产率与废料率之间的关系,记录了如表8.15所示的数据,请画出散点图,根据散点图的趋势拟合适当的回归模型。 表8.15 生产率x(单位/周) 1000 2000 3000 3500 4000 4500 5000 废品率y(%) 5.2 6.5 6.8 8.1 10.2 10.3 13.0

解:先画出散点图如下图:

eyAKLyAKL 从散点图大致可以判断出x和y之间呈抛物线或指数曲线,由此采用二次方程式和指数函数进行曲线回归。 (1)二次曲线 SPSS输出结果如下:

5000.004000.003000.002000.001000.00x

12.0010.008.006.00

y

Model Summary.981.962.942.651RR SquareAdjustedR SquareStd. Error ofthe EstimateThe independent variable is x.

ANOVA42.571221.28650.160.0011.6974.42444.2696

RegressionResidualTotal

Sum ofSquaresdfMean SquareFSig.

The independent variable is x. 从上表可以得到回归方程为:72ˆ5.8430.0874.4710yxx 由x的系数检验P值大于0.05,得到x的系数未通过显著性检验。 由x2的系数检验P值小于0.05,得到x2的系数通过了显著性检验。 (2)指数曲线

从上表可以得到回归方程为:0.0002tˆ

4.003ye

由参数检验P值≈0<0.05,得到回归方程的参数都非常显著。

Coefficients-.001.001-.449-.891.4234.47E-007.0001.4172.812.0485.8431.3244.414.012

xx ** 2(Constant)

BStd. ErrorUnstandardizedCoefficientsBetaStandardizedCoefficientstSig.

Model Summary.970.941.929.085RR SquareAdjustedR SquareStd. Error ofthe EstimateThe independent variable is x.

ANOVA.5731.57379.538.000.0365.007.6096

RegressionResidualTotal

Sum ofSquaresdfMean SquareFSig.

The independent variable is x.Coefficients

.000.000.9708.918.0004.003.34811.514.000x

(Constant)

BStd. ErrorUnstandardizedCoefficientsBetaStandardizedCoefficientstSig.

The dependent variable is ln(y). 从R2值,σ的估计值和模型检验统计量F值、t值及拟合图综合考虑,指数拟合效果更好一些。 8.3 已知变量x与y的样本数据如表8.16,画出散点图,试用αeβ/x来拟合回归模型,假设: (1) 乘性误差项,模型形式为y=αeβ/xeε

(2) 加性误差项,模型形式为y=αeβ/x+ε。 表8.16 序号 x y 序号 x y 序号 x y 1 4.20 0.086 6 3.20 0.150 11 2.20 0.350 2 4.06 0.090 7 3.00 0.170 12 2.00 0.440 3 3.80 0.100 8 2.80 0.190 13 1.80 0.620 4 3.60 0.120 9 2.60 0.220 14 1.60 0.940 5 3.40 0.130 10 2.40 0.240 15 1.40 1.620 解: 散点图:

(1) 乘性误差项,模型形式为y=αeβ/xeε 线性化:lny=lnα+β/x +ε 令y1=lny, a=lnα,x1=1/x . 做y1与x1的线性回归,SPSS输出结果如下: Model Summaryb

.999a.997.997.04783Model1RR SquareAdjustedR SquareStd. Error ofthe EstimatePredictors: (Constant), x1a. Dependent Variable: y1b. 从以上结果可以得到回归方程为:y1=-3.856+6.08x1 F检验和t检验的P值≈0<0.05,得到回归方程及其参数都非常显著。 回代为原方程为:y=0.021e6.08/x (2)加性误差项,模型形式为y=αeβ/x+ε 不能线性化,直接非线性拟合。给初值α=0.021,β=6.08(线性化结果),NLS结果如下:

ANOVAb10.930110.9304778.305.000a.03013.00210.96014

RegressionResidualTotalModel1

Sum ofSquaresdfMean SquareFSig.

Predictors: (Constant), x1a. Dependent Variable: y1b.

Coefficientsa

-3.856.037-103.830.0006.080.088.99969.125.000(Constant)x1Model1

BStd. ErrorUnstandardizedCoefficientsBetaStandardizedCoefficientstSig.

Dependent Variable: y1a.

Parameter Estimates.021.001.020.0236.061.0445.9656.157

ParameterabEstimateStd. ErrorLower BoundUpper Bound95% Confidence Interval 从以上结果可以得到回归方程为: y=0.021e6.061/x 根据R2≈1,参数的区间估计不包括零点且较短,可知回归方程拟合非常好,且其参数都显著。

8.4 Logistic函数常用于拟合某种消费品的拥有率,表8.17(书上239页,此处略)是北京市每百户家庭平均拥有的照相机数,试针对以下两种情况拟合Logistic回归函数。 (1)已知100u,用线性化方法拟合, (2)u未知,用非线性最小二乘法拟合。 解:(1),100u时,的线性拟合。对0111tybbu函数线性化得到:

ANOVAa4.45822.229.00113.0004.459152.46714

SourceRegressionResidualUncorrected TotalCorrected TotalSum ofSquaresdfMeanSquares

Dependent variable: yR squared = 1 - (Residual Sum of Squares) /(Corrected Sum of Squares) = 1.000.a. 11ln()1.8510.264100y0111ln()lnln100btby,令311ln()100yy,作

3y

关于t的线性回归分析,SPSS输出结果如下:

由表Model Summary得到,0.994R趋于1,回归方程的拟合优度好,由表ANOVA得到回归方程显著,由Coefficients表得到,回归系数都是显著的,得到方程:11ln()1.8510.264100y,进一步计算得到:

00.157b,10.768b(100u) 回代变量得到最终方程形式为: 1ˆ0.010.1570.768ty 最后看拟合效果,通过sequence画图:

Model Summaryb.994a.988.987.16820Model1RR SquareAdjustedR SquareStd. Error ofthe EstimatePredictors: (Constant), ta. Dependent Variable: y3b.

ANOVAb

39.839139.8391408.165.000a.48117.02840.32018

RegressionResidualTotalModel1

Sum ofSquaresdfMean SquareFSig.

Predictors: (Constant), ta. Dependent Variable: y3b.

Coefficientsa

-1.851.080-23.039.000-.264.007-.994-37.526.000(Constant)tModel1

BStd. ErrorUnstandardizedCoefficientsBetaStandardizedCoefficientstSig.

Dependent Variable: y3a.