非线回归.pptVIP

§8.3 非线性模型 §8.3 非线性模型 用非线性最小二乘法求得的三个参数估计值为 k=150.0，a=0.012 43，b=0.892 7 其中k=150.0为回归模型估计的国内航线里程增长上限。 图8.4是用Excel绘制的国内航线里程趋势预测图，其中粗实线是观测值，虚细线是预测值。 图8.4 龚珀兹曲线拟合国内航线里程趋势图 §8.3 非线性模型 上面的过程中回归迭代给出了一些警告错误，这是由于回归迭代过程中的参数取值超出了允许范围造成的，可以通过对参数取值范围增加一些限制解决，用非线性回归对话框中的Constraints按钮实现。 另外如果参数的初值给得不够准确回归迭代不收敛时，对参数增加一些限制可能就收敛了。例如本例中把k的初值改为130回归迭代不能收敛，但是增加一个k116的限制回归迭代就收敛了。 §8.3 非线性模型 龚珀兹模型和几种常见的非线性回归模型可以用三和值法求解，见参考文献[15]第13章。 在正态误差假定下，非线性回归的最小二乘估计与极大似然估计是相同的，而极大似然估计具有好的大样本性质，例如渐近无偏性、渐近正态性、一致性等。因而非线性最小二乘估计值比三和值更精确，可以把三和值法的参数估计值作为求解非线性最小二乘的初值。 §8.3 非线性模型 【例8.6】 下表8.9是我国从1950——2005年历年大陆总人口数，试用威布尔（Weibull）曲线拟合数据并做预测。威布尔曲线如下： 其中参数k是变量发展的上限，参数a 0, 0 b 1, c0。 §8.3 非线性模型 表8.9 我国历年大陆总人口数 单位：亿人 年份 t y 年份 t y 1950 1 5.5196 1978 29 9.6259 1951 2 5.63 1979 30 9.7542 1952 3 5.7482 1980 31 9.8705 1953 4 5.8796 1981 32 10.0072 1975 26 9.242 2003 54 12.9227 1976 27 9.3717 2004 55 12.9988 1977 28 9.4974 2005 56 13.0756 §8.3 非线性模型 根据人口学的专业预测，我国人口上限为16亿人，因此取k的初值=16，取b的初值=0.5，取c的初值=1。 对以上初值把t=1时（即1950年）=5.5196代入，得， 用以上初值做非线性最小二乘，得下面的输出结果8.7。从中看到，人口上限为k=15.76亿人，这与人口学预测的人口上限16亿人完全一致。图8.5是用Excel绘制的人口趋势预测图，其中粗实线是观测值，虚细线是预测值。 §8.3 非线性模型 §8.3 非线性模型 §8.3 非线性模型 图8.5 威布尔模型预测我国人口趋势图 §8.3 非线性模型 【例8.6】 柯布—道格拉斯生产函数研究。在计量经济学中有一种熟知的C-D（Cobb—Douglas）生产函数 其中，y为产出,K(资本)、L(劳力)为两个投入要素，A0为效率系数、a和b为K和L的产出弹性，A、a、b 均为待估参数。 §8.3 非线性模型 a是产出对资本投入的弹性系数，度量在劳动投入保持不变时资本投入增加1%时产出增加的百分比。 b是产出对劳动投入的弹性系数，度量在资本投入保持不变时劳动投入增加1%时产出增加的百分比。 两个弹性系数之和 a+b 表示规模报酬（returns to scale）。 a+b =1表示规模报酬不变，即1倍的投入带来1倍的产出； a+b ＜1表示规模报酬递减，即1倍的投入带来少于1倍的产出； a+b ＞1表示规模报酬递增，即1倍的投入带来大于1倍的产出。 §8.3 非线性模型 §8.3 非线性模型 年份 t GDP K L lnGDP lnK lnL 1978 1 3624.1 1377.9 40152 8.195361 7.228316 10.60043 1979 2 4038.2 1474.2 41024 8.303554 7.295871 10.62191 1980 3 4517.8 1590.0 42361 8.415780 7.371489 10.65398 1981 4 4862.4 1581.0 43725 8.489287 7.365813 10.68568 1982 5 5294.7 1760.2 45295 8.574462 7.473183 10.72095 1983 6 593