应用Excel建立二次饱及D-最优设计回归模型探究.docVIP

应用Excel建立二次饱及D-最优设计回归模型探究 　摘要总结了运用Excel建立二次饱和D-最优设计回归模型的过程,包括结构矩阵、转置矩阵、系数矩阵、逆矩阵、回归系数等参数的测算,指出该方法具有设计度高、计算简便等优点,可在农业生产的试验设计中推广应用。 　　关键词二次饱和D-最优试验设计;回归模型;Excel 　　中图分类号S11文献标识码A文章编号 1007-5739(2010)22-0020-02 　　 　　EstablishingSecondSaturatingD-OptimalDesignRegressionModelbyExcel 　　GAO Hong-ze 　　(Puyang Institute of Agricultural Sciences in Henan,Puyang Henan 457000) 　　AbstractThe process of establishing second saturating D-optimal design regression model by Excel was summarized,including the estimate of the parameter as structure matrix,transpose matrix,coefficient matrix,inverse matrix and regression coefficient and so on.This method had the advantage of high accuracy on design and easy to calculate,which could be spread and applied in the experimental design of agricultural production. 　　Key wordssecond saturating D-optimal design;regression model;Excel 　　 　　二次饱和D-最优设计是农业生产中常用的一种试验设计方法,可减少试验次数,提高试验精度,该方法的参数估计精确度高,但是其统计建模因运用矩阵计算而十分繁琐[1-3]。通过实践,利用Excel可轻松实现二次饱和D-最优试验设计[4]。 论文代写 　　1设计过程 　　1.1列出结构矩阵 　　打开Excel 2003,列出该文要分析的试验计划结构矩阵和试验结果(表1),把二次项放在一次项和交互项的中间,即可得到下面的结构矩阵表,用“X”表示: 　　1.2求转置矩阵 　　先求出结构矩阵X的转置矩阵,用X′表示。选中X结构矩阵表数据,单击“编辑”菜单上“复制”,然后单击“编辑”菜单上“选择性粘贴”,选中“转置”后点击“确定”,即可得到转置矩阵X′: 　　1.3利用MMULT函数求系数矩阵(信息矩阵)A 　　系数矩阵的计算公式为:A=X′X。MMULT函数是求两数组矩阵乘积的函数。在Excel表中,首先选中A矩阵要放置的区域,再单击“插入”,然后单击下拉菜单中的“函数”,出现插入函数对话框,在“插入函数”中选中MMULT函数,在函数参数Array1中选中X′矩阵数据组数据,在Array2中选中X矩阵数据组数据,按F2,再按 Ctrl+Shift+Enter,即可得到A矩阵数据: 　　1.4利用MINVERSE函数求A的称相关矩阵(逆矩阵)C 　　称相关矩阵的计算公式为:C=A-1。在Excel表中,首先选中C矩阵要放置的区域,单击“插入”,然后单击下拉菜单中的“函数”,出现插入函数对话框,在“插入函数”中选中MINVERSE函数,在函数参数Array中选中A矩阵数据组数据,按F2,再按Ctrl+Shift+Enter,即可得到C矩阵数据: 　1.5计算回归系数b 　　1.5.1用SUMPRODUCT函数计算常数项矩阵元素B。各列B值=Σ(该列每一编码值×对应的试验指标值)。 　　在Excel中的二因素二次饱和D-最优设计结构矩阵表和试验结果上进行计算。现以求B0进行示范,点中结构矩阵表B0值所放置的位置,再单击 “插入”,然后单击下拉菜单中的“函数”,出现插入函数对话框,在“插入函数”中选中SUMPRODUCT函数,在函数参数Array1中选中x0对应列的编码值,在Array2中选y列对应的数据组数据,然后点击“确定”即可得到B0的值。用同样的方法可求出其他B值。选中求得的B值,进行复制、转置粘贴,即可得到常数项矩阵B: 　　B=B0B1B2B11B22B12=3 253675.299 5361.090 72 222.4732 160.828-164.623 　　1.5.2利用MM