非线性回归课件.docxVIP

§3.4 非线性回归与回归诊断 一. 非线性回归 第一类非线性回归： y = a + be x + ε （3-71） y = a + b ln x + ε （3-72） y = a + bx 3 + ε y = a + b1 x + b2 x2 +L+bm xm + ε = a + b1 x1 + b2 x2 + b3 x1 x2 + b4 x12 + b5 x22 + ε 以上模型均假设ε ~ N (0，σ 2 ) 。通常可用适当变换将其化为线性的。 3-6 在维尼纶醛化试验中，固定其它因素 考虑醛浓度与反映时间对醛化度的关系，试验数据如表 3-12。记醛化度为y ，反映时间为x1 ，甲 醛浓度为x2 ，由经验知道，y 与x2 成正比、与x1 成反比，并有 表 3-12 试验数据表 醛 化 时 度 间 3 5 7 12 20 30 甲 醛 32.10 17.8 22.9 25.9 29.9 32.9 35.4 33.00 18.2 22.9 25.1 28.6 31.2 34.1 27.60 16.8 20.0 23.6 28.0 30.0 33.1 3-12 试验数据表 醛 化 时间 度 3 5 7 12 20 30 甲 醛 32.10 17.8 22.9 25.9 29.9 32.9 35.4 33.00 18.2 22.9 25.1 28.6 31.2 34.1 27.60 16.8 20.0 23.6 28.0 30.0 33.1 y = β + β 1 + β x + ε （3-73） 0 1 x1 2 2 试求β 0 、β1 、β 2 的最小二乘估计，并检验回归方程是否有意义（α = 0.05）。 分析：依题意及（3-73）式知，y 与x1 之间的关系不是线性的，但若令t1 = 1 / x1 ，则 y 与t1 之间 的关系是线性的。至于检验回归方程是否有意义 即检验假设H01 ：β 1 = β 2 = 0。 解：（1）参数估计：依题意令t1 = 1 / x1 ， t2 = x2 ，则（3-73）化为下列形式 = β 0 + β1t1 + β 2 t2 + ε 经计算，并将所得数据列入表 3-13, n n lij = ∑ xli xlj ? n x i x j = ∑( xli ? x i )( xlj ? x j ) = l ji l=1 l=1 n n l jy = ∑ yl xlj ? n x j y = ∑( xlj ? x j )( yl ? y ) l=1 l=1 n l yy = ∑ ( yl ? y)2 l=1 S回2 = l1 y b$1 + l2 yb$2 +L+lmyb$m 3-13 数据计算表 序 号 yi t1 t2 yi ? y ti1 ? t1 ti 2 ? t2 1 17.8 1/3 32.10 －8.67 0.1929 1.2 2 18.2 1/3 33.00 －8.27 0.1929 2.1 3 16.8 1/3 27.60 －9.67 0.1929 －3.3 4 22.9 1/5 32.10 －3.57 0.0595 1.2 5 22.9 1/5 33.00 －3.57 0.0595 2.1 6 20.0 1/5 27.60 －6.47 0.0595 －3.3 7 25.9 1/7 32.10 －0.57 0.0024 1.2 8 25.1 1/7 33.00 －1.37 0.0024 2.1 9 23.6 1/7 27.60 －2.87 0.0024 －3.3 10 29.9 1/12 32.10 3.43 －0.0571 1.2 11 28.6 1/12 33.00 2.13 －0.0571 2.1 12 28.0 1/12 27.60 1.53 －0.0571 －3.3 13 32.9 1/20 32.10 6.43 －0.0905 1.2 14 31.2 1/20 33.00 4.73 －0.0905 2.1 15 30.0 1/20 27.60 3.53 －0.0905 －3.3 16 35.4 1/30 32.10 8.93 －0.1071 1.2 17 34.1 1/30 33.00 7.63 －0.1071 2.1 18 33.1 1/30 27.60