: 针叶林比例(%); : 造林面积(千亩/年); : 年采伐面积(千亩/年); : 火灾频数(次/年)。
注: 本例子引自 裴鑫德 编著,《多元统计分析及其应用》, 北京农业大学出版社, 1990年
本软件给出的回归分析有关的结果如下(与回归分析无关的内容未列出):
指标1名称: 森林覆盖率 单位: % 指标2名称: 针叶林比例 单位: % 指标3名称: 造林面积 单位: 万亩/年 指标4名称: 年采伐面积 单位: 千亩/年 指标5名称: 火灾频数 单位: 次/年 因素1名称: 山地比例 单位: %
因素2名称: 人口密度 单位: 人/平方公里 因素3名称: 人均收入增长率 单位: 元/年 因素4名称: 公路密度 单位: 100米/公顷 因素5名称: 前汛期降水量 单位: 厘米/年 因素6名称: 后汛期降水量 单位: 厘米/年 因素7名称: 月平均最低温度 单位: ℃
回归分析采用双重筛选逐步回归法, 显著性水平α=0.05 自变量引入、剔除的临界值Fx=2.000 因变量引入、剔除的临界值Fy=2.500
对指标1~5拟建立回归方程分别为:
y1 = b(0) + b(1)*X(1) + b(2)*X(2) + b(3)*X(3) + b(4)*X(4) + b(5)*X(5) + b(6)*X(6) + b(7)*X(7)
y2 = b(0) + b(1)*X(1) + b(2)*X(2) + b(3)*X(3) + b(4)*X(4) + b(5)*X(5) + b(6)*X(6) + b(7)*X(7)
y3 = b(0) + b(1)*X(1) + b(2)*X(2) + b(3)*X(3) + b(4)*X(4) + b(5)*X(5) + b(6)*X(6) + b(7)*X(7)
y4 = b(0) + b(1)*X(1) + b(2)*X(2) + b(3)*X(3) + b(4)*X(4) + b(5)*X(5) + b(6)*X(6) + b(7)*X(7)
y5 = b(0) + b(1)*X(1) + b(2)*X(2) + b(3)*X(3) + b(4)*X(4) + b(5)*X(5) + b(6)*X(6) + b(7)*X(7)
----- 计算第1组回归方程 -----
第1步, 引入方程项: y1 已引入因变量y的序号: 1 已引入自变量X的序号:
第2步, 自变量X引入或剔除判别: 各项的判别值(升序排列): Vx(3)= 4.541e-2 Vx(5)= 0.2868 Vx(7)= 0.4082 Vx(4)= 0.4104 Vx(6)= 0.4731 Vx(2)= 0.5998 Vx(1)= 0.8810
未引入项中, 第1项[X(1)]Vx值(≥0)的绝对值最大, 引入检验值Fax(1)=74.00, 引入临界值Fx=2.000, Fax(1)>Fx, 可以引入第1项。
第3步, 引入方程项: X(1) 已引入因变量y的序号: 1 已引入自变量X的序号: 1
第4步, 自变量X引入或剔除判别: 各项的判别值(升序排列): Vx(1)=-7.400 Vx(5)= 1.359e-3 Vx(6)= 2.254e-2 Vx(3)= 4.720e-2 Vx(2)= 0.2260 Vx(7)= 0.2306 Vx(4)= 0.2372
已引入项中, 第1项[X(1)]Vx值(<0)的绝对值最小, 剔除检验值Fex(1)=74.00, 剔除临界值Fx=2.000,
Fex(1)>Fx, 不能剔除第1项, 检查是否可以引入其他自变量。 未引入项中, 第4项[X(4)]Vx值(≥0)的绝对值最大, 引入检验值Fax(4)=2.798, 引入临界值Fx=2.000, Fax(4)>Fx, 可以引入第4项。
第5步, 引入方程项: X(4) 已引入因变量y的序号: 1 已引入自变量X的序号: 1,4
第6步, 自变量X引入或剔除判别: 各项的判别值(升序排列): Vx(1)=-5.492
Vx(4)=-0.3109 Vx(5)= 1.370e-8 Vx(6)= 6.972e-3 Vx(2)= 4.284e-2 Vx(3)= 9.849e-2 Vx(7)= 0.2536
已引入项中, 第4项[X(4)]Vx值(<0)的绝对值最小, 剔除检验值Fex(4)=2.798, 剔除临界值Fx=2.000,
Fex(4)>Fx, 不能剔除第4项, 检查是否可以引入其他自变量。 未引入项中, 第7项[X(7)]Vx值(≥0)的绝对值最大, 引入检验值Fax(7)=2.718, 引入临界值Fx=2.000, Fax(7)>Fx, 可以引入第7项。
第7步, 引入方程项: X(7) 已引入因变量y的序号: 1 已引入自变量X的序号: 1,4,7
第8步, 自变量X引入或剔除判别: 各项的判别值(升序排列): Vx(1)=-4.767 Vx(4)=-0.3513 Vx(7)=-0.3398 Vx(2)= 5.297e-2 Vx(5)= 6.120e-2 Vx(6)= 0.1234 Vx(3)= 0.1380
已引入项中, 第7项[X(7)]Vx值(<0)的绝对值最小, 剔除检验值Fex(7)=2.718, 剔除临界值Fx=2.000,
Fex(7)>Fx, 不能剔除第7项, 检查是否可以引入其他自变量。 未引入项中, 第3项[X(3)]Vx值(≥0)的绝对值最大, 引入检验值Fax(3)=1.120, 引入临界值Fx=2.000,
Fax(3)≤Fx, 不能引入第3项, 检查是否可以引入其他因变量。
第9步, 因变量y引入或剔除判别: 各项的判别值(升序排列): Vy(1)=-13.75 Vy(5)= 0.2198 Vy(3)= 0.2859 Vy(2)= 0.5592 Vy(4)= 0.5895
已引入项中, 第1项[y1]Vy值(<0)的绝对值最小, 剔除检验值Fey(1)=36.68, 剔除临界值Fy=2.500,
Fey(1)>Fy, 不能剔除第1项, 检查是否可以引入其他因变量。 未引入项中, 第4项[y4]Vy值(≥0)的绝对值最大,
引入检验值Fay(4)=3.350, 引入临界值Fy=2.500, Fay(4)>Fy, 可以引入第4项。
第10步, 引入方程项: y4 已引入因变量y的序号: 1,4 已引入自变量X的序号: 1,4,7
第11步, 自变量X引入或剔除判别: 各项的判别值(升序排列): Vx(1)=-8.013 Vx(4)=-0.8942 Vx(7)=-0.7008 Vx(6)= 0.1388 Vx(2)= 0.3799 Vx(5)= 0.4204 Vx(3)= 0.4595
已引入项中, 第7项[X(7)]Vx值(<0)的绝对值最小, 剔除检验值Fex(7)=2.453, 剔除临界值Fx=2.000,
Fex(7)>Fx, 不能剔除第7项, 检查是否可以引入其他自变量。 未引入项中, 第3项[X(3)]Vx值(≥0)的绝对值最大, 引入检验值Fax(3)=2.551, 引入临界值Fx=2.000, Fax(3)>Fx, 可以引入第3项。
第12步, 引入方程项: X(3) 已引入因变量y的序号: 1,4 已引入自变量X的序号: 1,3,4,7
第13步, 自变量X引入或剔除判别: 各项的判别值(升序排列): Vx(1)=-7.721 Vx(7)=-1.304 Vx(4)=-0.8840 Vx(3)=-0.8502 Vx(6)= 0.1487 Vx(5)= 0.3635 Vx(2)= 0.5014
已引入项中, 第3项[X(3)]Vx值(<0)的绝对值最小, 剔除检验值Fex(3)=2.551, 剔除临界值Fx=2.000,
Fex(3)>Fx, 不能剔除第3项, 检查是否可以引入其他自变量。 未引入项中, 第2项[X(2)]Vx值(≥0)的绝对值最大, 引入检验值Fax(2)=2.514, 引入临界值Fx=2.000, Fax(2)>Fx, 可以引入第2项。
第14步, 引入方程项: X(2)
已引入因变量y的序号: 1,4
已引入自变量X的序号: 1,2,3,4,7
第15步, 自变量X引入或剔除判别: 各项的判别值(升序排列): Vx(1)=-6.973 Vx(7)=-1.798 Vx(3)=-1.301 Vx(2)=-1.006 Vx(4)=-9.920e-2 Vx(6)= 0.1563 Vx(5)= 0.5073
已引入项中, 第4项[X(4)]Vx值(<0)的绝对值最小, 剔除检验值Fex(4)=0.2480, 剔除临界值Fx=2.000, Fex(4)≤Fx, 需要剔除第4项。
第16步, 剔除方程项: X(4) 已引入因变量y的序号: 1,4 已引入自变量X的序号: 1,2,3,7
第17步, 自变量X引入或剔除判别: 各项的判别值(升序排列): Vx(1)=-6.465 Vx(2)=-2.438 Vx(7)=-1.794 Vx(3)=-1.293 Vx(4)= 9.025e-2 Vx(6)= 0.1959 Vx(5)= 0.4328
已引入项中, 第3项[X(3)]Vx值(<0)的绝对值最小, 剔除检验值Fex(3)=3.878, 剔除临界值Fx=2.000,
Fex(3)>Fx, 不能剔除第3项, 检查是否可以引入其他自变量。 未引入项中, 第5项[X(5)]Vx值(≥0)的绝对值最大, 引入检验值Fax(5)=1.907, 引入临界值Fx=2.000,
Fax(5)≤Fx, 不能引入第5项, 检查是否可以引入其他因变量。
第18步, 因变量y引入或剔除判别: 各项的判别值(升序排列): Vy(1)=-16.34 Vy(4)=-6.328 Vy(5)= 0.2662 Vy(2)= 0.4603 Vy(3)= 0.6480
已引入项中, 第4项[y4]Vy值(<0)的绝对值最小,