三元回归.doc

应用FX-4000P计算器进行三元回归分析的方法姚永兰1，门海元2(1.原平市绿委，山西 原平 034100；2.山西省林业调查规划院，山西 太原 0.0012) 摘要:文中根据FX-4000P计算器的编程原理，编制了林业上常用三元回归方程的通解程序。输出方程参数，相关系数及方差分析表，并通过实例操作系统地叙述了计算的全过程。关鍵词: 计算器 三元回归分析1.概述三元回归分析方法是林业生产和科研中分析处理变量之间相互关系的数学方法，在研究变量之间相互关系时，常按照既定的目标选择变量，收集数据。建立合乎实际的回归方程。并对回归方程的显著性进行检验。而用回归分析的方法找到的变量之间的关系式，并不是确定性的，而是近似性的。通常称为经验公式，找出的这个函数称为回归函数。因此，研究回归关系的主要之点，是要确定回归函数。在林业生产中常会遇到这样的关系，林木的材积与胸径、树高和形数有关；造林的成活率与水分、温度和光照有关；林木的结实量与光照、气候和降水量有关。研究这一类型的问题，就可以使用三元回归的方式来解决。本文求解的三元回归方程有：1. Y=B0+B1X1+B2X2+B3X32.D=B0+B1 (1/A)+B2(1/A)2+B3(1/A)3或者PV= B0+B1 (1/D)+B2(1/D)2+B3(1/D)33.Y=B0+B1X+B2X2+B3X32.基本原理设自变量为X1，X2，X3，因变量为Y，样本数为N。L11=L22=L33=L12=L13=L23=L1Y=L2Y=L3Y=LYY=用最小二乘法整理得正规方程组：L11B1+L12B2+L13B3=L1Y (1)L21B1+L22B2+L23B3=L2Y (2)L31B1+L32B2+L33B3=L3Y (3)用高斯消元法解方程组得参数：B1、B2 、B3。B0=-B11- B22- B33复相关系数为：R=统计量(F)： F=UK(Q(N-K-1)方差： SU2=UK SQ2=Q(N-K-1)式中：U回归平方和；Q-残差平方和；K-自变量个数；N-为样本单元数；SU2-回归方差；SQ2-残差方差。方程2、3都可转化为二元线性回归 Y=B0+B1X1+B2X2+B3X3的一般形式来计算。2.计算程序MCl:“A”:？A:“N”：？N：Lbl 1: MN= GOTO 2:“Y”：？W：A=2= GOTO 4: A=3= GOTO 5: “X1”: ？X: “X2”: ？Y：“X3”: ？Z：Lbl 6:E+WE：F+XF：G+YG：H+ZH：W2+TT：I+X2I：J+Y2J：K+Z2K：XY+LL：X+Z+OO：YZ+PP：XW+QQ：YW+RR：ZW+SS：T-N-1E2D：M+1M：GOTO 1: Lbl 4: “X”: ？X:X-1X：X2Y：YXZ：GOTO 6: Lbl 5: “X”: ？X: X2Y：YXZ：GOTO 6: Lbl 2:I-N-1F2I：J-N-1G2J：K-N-1H2K：L-N-1FGL：O-N-1FHO：P-N-1GHP：Q-N-1EFQ：R-N-1EGR：S-N-1EHS：-LIB：LB+JA：OB+PM：QB+RU：-OIOT：-MAC：K+T+MCW：S+Q(-OI)+UCX：XWY：(U-(MXW)AZ：(Q-LZ-OY) IV：EN-1-FVN-1-GZN-1-HYN-1EVZYVQ+ZR+YSF：(FD)FD-FDF3(D-F)(N-4)PF3P3.应用实例3.1 某自变量X1、X2和X3与因变量Y的观测值见表1，样本数7，试建立三元回归模型。并对模型进行检验。表1 样本资料观测值编号YX1X2X311637680.74226.140742.12321.442632.1426.536723524.343691.67计算操作过程鍵盘操作显示说 明prg N EXEA?启动程序，询问选择第几个回归方程？1 EXEN?输入A，询问样本数？5 EXEY?输入N，询问因变量？16 EXEX1输入Y，询问第一个自变量？37 EXEX2输入X1，询问第二个自变量？68 EXEX3输入X2，询问第三个自变量？0.74 EXEY?输入X3，继续询问下一个因变量？1.67 EXE-47.2458输入最后一个X2，显示回归参数B。 EXE0.62353显示回归参数B1。 EXE0.54886显示回归参数B2。 EXE3.85907显示回归参数B3。 EXE0.99719显示复相关系数R。 EXE74.5914显示回归平方和。 EXE0.42056显示残差平方和。 EXE75.012显示总离差平方和。 EXE24.8638显示回归方差。 EXE0.4206显示残差方差。 EXE59.1208显示统计量。计算结果是：I.回归模型：Y=-47.2458+0.62353X1+0.54886X2+3.85907X3。II.复相关系数：R=0.99719。III.方差分析见表2。表2 方差分析表误差来源离差平方和方差F回归74.591424.863859.1208残差 0.420560.4206总的75.0123.2 某林分材积生长率和平均年龄的调查结果见表3。试用三次多项式D=B0+B1 (1/A)+B2(1/A)2+B3(1/A)3拟合数学模型，并进行检验。表3 林分调查原始记录材积生长率38.118.69.83.31.51.1林分年龄81318232830计算操作过程鍵盘操作显示说 明prg N EXEA?启动程序，询问选择第几个回归方程？2 EXEN?输入A，询问样本数？6 EXEY?输入N，询问因变量？38.1 EXEX?输入Y，询问自变量？8 EXEY?输入X，继续询问下一个因变量？1.1 EXEX?输入y，询问自变量？30 B EXE-3.9131输入x，显示回归参数B0。 EXE-53.826显示回归参数B1。 EXE6809.31显示回归参数B2。 EXE-29525.1显示回归参数B3。 EXE0.99927显示复相关系数R。 EXE1032.82显示回归平方和。 EXE1.5098显示残差平方和。 EXE1034.33显示总离差平方和。 EXE344.27显示回归方差。 EXE0.7549显示残差方差。 EXE456.058显示统计量。运行结果是：I.回归模型：PV=-3.9131-53.826/A+6809.31/A2-29525.1/A3。II.复相关系数：R=0.99927。III.方差分析见表4。表4 方差分析表误差来源离差平方和方差F回归1032.82344.27456.058残差 1.50980.7549总的1034.333.3 某林分平均胸径与平均年龄的实测值见表5，试用模型：Y=B0+B1X+B1X2+B2X3拟合回归方程，并进行检验。表5 林分调查原始记录平均胸径15.114.813.811.78.42.1平均年龄30252015105计算操作过程鍵盘操作显示说 明prg N EXEA?启动程序，询问选择第几个回归方程？3 EXEN?输入A，询问样本数？6 EXEY?输入N，询问因变量？15.1 EXEX?输入Y，询问自变量？30 EXEY?输入X，继续询问下一个因变量？2.1 EXEX?输入y，询问自变量？5 EXE-6.933输入x，显示回归参数B0。 EXE2.1658显示回归参数B1。 EXE-0.07389显示回归参数B2。 EXE-0.00087显示回归参数B3。 EXE0.9997显示复相关系数R。 EXE125.46显示回归平方和。 EXE0.08587显示残差平方和。 EXE125.55显示总离差平方和。 EXE41.821显示回归方差。 EXE0.04294显示残差方差。 EXE974.017显示统计量。运行结果是：I.回归模型：-6.933+2.1658 X0.07389X20.00087 X3。II.复相关系数：R=0.9997。III.方差分析见表6。表6 方差分析表误差来源离差平方和方差F回归