U(1)=3.612 U(2)=11.33 U(3)=17.66 U(4)=1.211
偏相关系数 ρ(i): ρ1,234=-0.5171 ρ2,134=-0.7307 ρ3,124= 0.8005 ρ4,123= 0.3302
各方程项对回归的贡献(按偏回归平方和降序排列): U(3)=17.66, U(3)/U=67.02% U(2)=11.33, U(2)/U=43.01% U(1)=3.612, U(1)/U=13.71% U(4)=1.211, U(4)/U=4.596%
残差分析:
残 差 分 析 表 № 观 测 值 1 1.000 2 0.000 3 0.000 4 3.000 5 2.000 6 0.000 7 4.000 8 5.000 9 0.000 10 0.000 11 0.000 12 0.000
----- 计算第3组回归方程 -----
尚未选入的因变量是y3
第33步, 引入方程项: y3 已引入因变量y的序号: 3 已引入自变量X的序号:
第34步, 自变量X引入或剔除判别: 各项的判别值(升序排列):
回 归 值 0.3945 0.6236 2.290 1.302 2.252 -8.343e-2 3.872 4.162 -0.1488 0.3906 0.1278 -0.1829 观测值-回归值 0.6055 -0.6236 -2.290 1.698 -0.2520 8.343e-2 0.1280 0.8380 0.1488 -0.3906 -0.1278 0.1829 (回归值-观测值)/观测值×100(%) -60.55 - - -56.60 12.60 - -3.200 -16.76 - - - - Vx(4)= 2.044e-2 Vx(6)= 3.077e-2 Vx(7)= 6.082e-2 Vx(2)= 6.558e-2 Vx(3)= 8.076e-2 Vx(5)= 9.328e-2 Vx(1)= 0.1769
未引入项中, 第1项[X(1)]Vx值(≥0)的绝对值最大, 引入检验值Fax(1)=2.149, 引入临界值Fx=2.000, Fax(1)>Fx, 可以引入第1项。
第35步, 引入方程项: X(1) 已引入因变量y的序号: 3 已引入自变量X的序号: 1
第36步, 自变量X引入或剔除判别: 各项的判别值(升序排列): Vx(1)=-0.2149 Vx(7)= 9.151e-4 Vx(2)= 3.530e-3 Vx(5)= 6.790e-3 Vx(6)= 3.149e-2 Vx(3)= 0.2267 Vx(4)= 0.2269
已引入项中, 第1项[X(1)]Vx值(<0)的绝对值最小, 剔除检验值Fex(1)=2.149, 剔除临界值Fx=2.000,
Fex(1)>Fx, 不能剔除第1项, 检查是否可以引入其他自变量。 未引入项中, 第4项[X(4)]Vx值(≥0)的绝对值最大, 引入检验值Fax(4)=2.641, 引入临界值Fx=2.000, Fax(4)>Fx, 可以引入第4项。
第37步, 引入方程项: X(4) 已引入因变量y的序号: 3 已引入自变量X的序号: 1,4
第38步, 自变量X引入或剔除判别: 各项的判别值(升序排列): Vx(1)=-0.5393 Vx(4)=-0.2935 Vx(5)= 2.782e-3 Vx(7)= 6.635e-3 Vx(6)= 1.355e-2 Vx(2)= 0.2129 Vx(3)= 0.2336
已引入项中, 第4项[X(4)]Vx值(<0)的绝对值最小, 剔除检验值Fex(4)=2.641, 剔除临界值Fx=2.000,
Fex(4)>Fx, 不能剔除第4项, 检查是否可以引入其他自变量。 未引入项中, 第3项[X(3)]Vx值(≥0)的绝对值最大, 引入检验值Fax(3)=2.439, 引入临界值Fx=2.000, Fax(3)>Fx, 可以引入第3项。
第39步, 引入方程项: X(3) 已引入因变量y的序号: 3 已引入自变量X的序号: 1,3,4
第40步, 自变量X引入或剔除判别: 各项的判别值(升序排列): Vx(1)=-0.8738 Vx(4)=-0.3051 Vx(3)=-0.3049 Vx(7)= 1.016e-2 Vx(6)= 3.539e-2 Vx(5)= 8.790e-2 Vx(2)= 0.2678
已引入项中, 第3项[X(3)]Vx值(<0)的绝对值最小, 剔除检验值Fex(3)=2.439, 剔除临界值Fx=2.000,
Fex(3)>Fx, 不能剔除第3项, 检查是否可以引入其他自变量。 未引入项中, 第2项[X(2)]Vx值(≥0)的绝对值最大, 引入检验值Fax(2)=2.560, 引入临界值Fx=2.000, Fax(2)>Fx, 可以引入第2项。
第41步, 引入方程项: X(2) 已引入因变量y的序号: 3
已引入自变量X的序号: 1,2,3,4
第42步, 自变量X引入或剔除判别: 各项的判别值(升序排列): Vx(4)=-0.7808 Vx(3)=-0.4028 Vx(2)=-0.3657 Vx(1)=-0.3559 Vx(7)= 1.163e-2 Vx(6)= 1.640e-2 Vx(5)= 9.282e-2
已引入项中, 第1项[X(1)]Vx值(<0)的绝对值最小, 剔除检验值Fex(1)=2.491, 剔除临界值Fx=2.000,
Fex(1)>Fx, 不能剔除第1项, 检查是否可以引入其他自变量。 未引入项中, 第5项[X(5)]Vx值(≥0)的绝对值最大,
引入检验值Fax(5)=0.6139, 引入临界值Fx=2.000,
Fax(5)≤Fx, 不能引入第5项, 检查是否可以引入其他因变量。
第43步, 因变量y引入或剔除判别: 各项的判别值(升序排列): Vy(3)=-1.800
已引入项中, 第3项[y3]Vy值(<0)的绝对值最小, 剔除检验值Fey(3)=3.151, 剔除临界值Fy=2.500, Fey(3)>Fy, 不能剔除第3项, 变量筛选结束。
变量筛选结果:
因变量: 总 数=1, 引 入 数=1
自变量: 检验项数=7, 预期引入项数=4, 实际引入项数=4, 实际引入项数=预期引入项数
第3组回归方程结果:
回归方程:
y3 = b(0) + b(1)*X(1) + b(2)*X(2) + b(3)*X(3) + b(4)*X(4)
回归系数 b(i): b(0)=-38.98 b(1)= 0.9282 b(2)=-0.2961 b(3)=-2.865 b(4)= 26.02
标准回归系数 B(i): B(1)= 0.5097 B(2)=-0.7273 B(3)=-0.4006 B(4)= 0.8986
复相关系数 γ=0.8018 决定系数 γ^2=0.6429 调整的决定系数 γ^2a=0.5090
变量分析:
变 量 分 析 表 变异来源 回 归 剩 余 总 和 平 方 和 U=3877 Q=2153 L=6030 自 由 度 K=4 N-1-K=7 N-1=11 均 方 U/K=969.2 Q/(N-1-K)=307.6 均 方 比 F=3.151 样本容量N=12, 显著性水平α=0.05, 检验值Ft=3.151, 临界值F(0.05,4,7)=4.120
剩余标准差 s=17.54
回归系数检验值: t检验值(df=7): t(1)= 1.578 t(2)=-1.600 t(3)=-1.679 t(4)= 2.338
F检验值(df1=1, df2=7): F(1)= 2.491 F(2)= 2.560 F(3)= 2.819 F(4)= 5.466
偏回归平方和 U(i): U(1)=766.3 U(2)=787.5 U(3)=867.3 U(4)=1681
偏相关系数 ρ(i): ρ1,234= 0.5123 ρ2,134=-0.5175 ρ3,124=-0.5358 ρ4,123= 0.6622
各方程项对回归的贡献(按偏回归平方和降序排列): U(4)=1681, U(4)/U=43.37% U(3)=867.3, U(3)/U=22.37% U(2)=787.5, U(2)/U=20.31% U(1)=766.3, U(1)/U=19.77%
残差分析:
残 差 分 析 表 № 观 测 值 1 15.39 2 10.84 3 13.82 4 34.57 5 7.410 6 9.120 7 8.640 8 81.64 回 归 值 20.37 14.10 40.57 33.80 4.595 -1.222 25.12 54.93 观测值-回归值 -4.980 -3.260 -26.75 0.7700 2.815 10.34 -16.48 26.71 (回归值-观测值)/观测值×100(%) 32.36 30.07 193.6 -2.227 -37.99 -113.4 190.7 -32.72 9 23.95 10 33.60 11 8.930 12 58.97
23.38 17.04 14.04 60.14 0.5700 16.56 -5.110 -1.170 -2.380 -49.29 57.22 1.984 得到3组回归方程。
显著性水平α=0.05, 自变量引入、剔除的临界值43步计算后得到三组回归方程: 第一组: 第二组: 第三组:
由计算结果看出, 森林覆盖率人口密度密度
、人均收入增长率
及年采伐面积
受相同自变量影响, 主要影响因素为山地比例的影响; 针叶林比例
及火灾频数
、
。 , , , ,
, 因变量引入、剔除的临界值
, 经
及月平均最低气温主要受公路
的
及月平均最低气温影响; 造林面积主要受公路密度、人均收入增长率及山地比例
影响。