应用回归分析第九章部分答案word精品文档7页

第9章 非线性回归

9.1 在非线性回归线性化时，对因变量作变换应注意什么问题？

答：在对非线性回归模型线性化时，对因变量作变换时不仅要注意回归函数的形式， 还要注意误差项的形式。如：

(1) 乘性误差项，模型形式为

e y AK L αβε=， (2) 加性误差项，模型形式为y AK L αβε=+

对乘法误差项模型（1）可通过两边取对数转化成线性模型，（2）不能线性化。

一般总是假定非线性模型误差项的形式就是能够使回归模型线性化的形式，为了方便通常省去误差项，仅考虑回归函数的形式。

9.2为了研究生产率与废料率之间的关系，记录了如表9.14所示的数据，请画出散点图，根据散点图的趋势拟合适当的回归模型。

表9.14 生产率x （单位/周） 1000 2000 3000 3500 4000 4500 5000 废品率y （%） 5.2 6.5 6.8 8.1 10.2 10.3 13.0 解：先画出散点图如下图：

从散点图大致可以判断出x 和y 之间呈抛物线或指数曲线，由此采用二次方程式和指数函数进行曲线回归。

（1）二次曲线

SPSS 输出结果如下：

从上表可以得到回归方程为：72ˆ 5.8430.087 4.4710y

x x -=-+⨯ 由x 的系数检验P 值大于0.05，得到x 的系数未通过显著性检验。 由x 2的系数检验P 值小于0.05，得到x 2的系数通过了显著性检验。

（2）指数曲线

从上表可以得到回归方程为：0.0002t ˆ 4.003y

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

从R2值，σ的估计值和模型检验统计量F值、t值及拟合图综合考虑，指数拟合效果更好一些。

9.3 已知变量x与y的样本数据如表9.15，画出散点图，试用αeβ/x来拟合回归模型，假设：

(1)乘性误差项，模型形式为y=αeβ/x eε

(2)加性误差项，模型形式为y=αeβ/x+ε。

表9.15

序号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β/x eε

线性化：lny=lnα+β/x +ε令y1=lny, a=lnα,x1=1/x .

做y1与x1的线性回归，SPSS输出结果如下：

从以上结果可以得到回归方程为：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结果如下：

从以上结果可以得到回归方程为：y=0.021e6.061/x

根据R2≈1，参数的区间估计不包括零点且较短，可知回归方程拟合

非常好，且其参数都显著。

9.4 Logistic 回归函数常用于拟合某种消费品的拥有率，表8.17（书上239页，此处略）是北京市每百户家庭平均拥有的照相机数，试针对以下两种情况拟合Logistic 回归函数。

（1）已知100u =，用线性化方法拟合，

（2）u 未知，用非线性最小二乘法拟合。根据经济学的意义知道，u 是拥有率的上限，初值可取100；b0>0,0

解：（1），100u =时，的线性拟合。对011

1t y b b u =+函数线性化得到：

11ln() 1.8510.264100y -=--0111ln()ln ln 100b t b y -=+，令311ln()100y y =-，作3y 关于t 的线性回归分析，SPSS 输出结果如下：

由表Model Summary 得到，0.994R =趋于1，回归方程的拟合优度好，由表ANOVA 得到回归方程显著，由Coefficients 表得到，回归系数都是显著的，得到方程：11ln() 1.8510.264100

y -=--，进一步计算得到：00.157b =，10.768b =（100u =）

回代变量得到最终方程形式为： 1ˆ0.010.1570.768

t y =+⨯ 最后看拟合效果，通过sequence 画图：

由图可知回归效果比较令人满意。

（2）非线性最小二乘拟合,取初值100u =,00.157b =,10.768b =： 一共循环迭代8次，得到回归分析结果为：

0.995R =>0.994，得到回归效果比线性拟合要好，且：91.062u =，00.211b =，10.727b =，

回归方程为：1

1

0.211*0.727 91.062t

y=

+

。

最后看拟合效果，由sequence画图：

得到回归效果很好，而且较优于线性回归。

9.5表9.17（书上233页，此处略）数据中GDP 和投资额K 都是用定基居民消费价格指数（CPI ）缩减后的，以1978年的价格指数为100。

（1） 用线性化乘性误差项模型拟合C-D 生产函数；

（2） 用非线性最小二乘拟合加性误差项模型的C-D 生产函数；

（3） 对线性化检验自相关，如果存在自相关则用自回归方法改进；

（4） 对线性化检验多重共线性，如果存在多重共线性则用岭回归方法改进； 解：（1）对乘法误差项模型可通过两边取对数转化成线性模型。

ln y =ln A + α ln K + β ln L

令y ′=ln y ,β

0=ln A ,x 1=ln K ,x 2=ln L ,

y ′=β0+ α x 1+ β

x 2+ ε SPSS 输出结果如下：

模型综述表 从模型综述表中可以看到，调整后的为0.993，说明C-D 生产函数拟合效果很好，也说明GDP 的增长是一个指数模型。

方差分析表

从方差分析表中可以看到，F 值很大，P 值为零，说明模型通过了检验，这与上述分析结果一致。

系数表

根据系数表显示，回归方程为：

尽管模型通过了检验，但是也可以看到，常数项没有通过检验，但在这个模型里，当lnK 和lnL 都为零时，lnY 为-1.785，即当K 和L 都为1时，GDP 为0.168，也就是说当投入资本和劳动力都为1个单位时，GDP 将增加0.168个单位，这种解释在我们的承受范围内，可以认为模型可以用。

最终方程结果为：

y=0.618K 0.801 L 0.404

（2） 用非线性最小二乘法拟合加性误差项模型的C-D 生产函数；

上述假设误差是乘性的，现假设误差是加性的情况下使用非线性最小二乘法估计。初值采用（1）中参数的结果，SPSS 输出结果如下：

参数估计表

SPSS 经过多步迭代，最终得到的稳定参数值为P=0.407，a=0.868，b=0.270

y=0.407K 0.868 L 0.270