第22卷第2期 四川农业大学学报 Vol.22 No.2 2004年06月 Journal of Sichuan Agricultural University Jum.2004 Computer Programs of the Method of Combined Analysis for the Data of Crop Variety Regional Trials in Several Years LIU Yong-jian,MING Dao-xu.HUANG YU-bi,FAN Zhi-xia (College of Agriculture,Sichun Agricultural University.Yaan 625014,Sichuan,China Abstract:Computer programs were compiled to analyze the data of crop variety regional trials in several years by means of Visual BASIC 6.0 and the SAS System for Windows version 6.12.By debugging the two programs,it indicates that analysis of variance can be done,the least squares means (LSM)of varieties,sites and their combinations can be estimated and the significance of difference between LSMs can be tested.The results analyzed by the two programs are equal correspondingly Key words:crop variety regional trials:computer program 作物品种区域试验多年资料联合分析法的计算机程序 刘永建，明道绪，黄玉碧，范志霞 四川农业大学农学院，四川雅安625014 摘要：分别用Visual BASIC6.0和The SAS System for Windows V6.12编制了作物品种区域试验多年资料联合分析 方法的计算机程序。程序调试结果表明，两种系统开发的程序不仅可以进行方差分析，而且可以计算品种、试点和 品种试点组合的最小二均数，并可对其差异显著性进行检验，且计算结果完全相同。 关健词：作物品种区域试验：计算机程序 中图分类号：S11*4:S322.1 文就标识码，A 文章编号：1000-2650(2004)02-0095-02 作物品种区域试验由于参试品种的参试年份、 于评价参试品种的生产能力及其适应区域。因此， 参试地点可能不完全相同，其至还可能出现缺区的 可以将方差分析数学模型 情况，造成试验资料的不平衡性。对于这类资料，明 =红+a:+月+Ya+(a3)g+(ay)a+ 道绪等根据最小二乘原理提出了相应的分析方 (1) 法。但是，至今为止还没有该方法的实际应用 (BY)+(a3)+t+e脚 简化为 为了便于运用该方法分析作物品种区域试验多年资 料，对参试品种作出科学评价，本文拟利用Vsul y=μ+a+月+Ys+（a3)g+（)+ BASIC6.The SAS System for Windows V6. t+e测 (2) 12编制该方法的计算机程序。 在(1)和(2)中，y为第i年份第j地区第k品种第 l区组的观测值：以为总体平均数：a为第i年份效 1数学模型 应：B,为第j地区效应：Y为第k品种效应：(a3) 作物品种区域试验由于参试年份、参试地区可 为第i年份第j地区的互作效应：(ay),为第i年份 能不同，有时可能出现缺区，所得到的资料可能是非 第k品种的互作效应：()，：为第j地区第k品种的 平衡资料。利用明道绪等提出的联合分析方法 互作效应：()为第i年份第j地区第k品种的 分析这类资料时，必须满足观测值个数大于待估计 互作效应：为第i年份第j地区内第1区组效应 的独立参数个数，而且作物品种区域试验的目的在 e抛为随机误差、相互独立且都服从N(0,g2):为 ·袋雅器自袋科辛特金项用心307
!"#$%&’()("*(+#,"-&.’/’&."0"-!"#123’043+56,2, -"(&.’7+&+"-!("$8+(2’&69’*2"3+5:(2+5,23;’&：!"$($*"32$(3)"06E:7EC4?$A+%%)0% ,.$,>"=("%(/$$0H-I3 */0A$,$3,$?EJ.$($3+#,3/0/#4@$?A4,.$,>"=("%(/3B6E:7编制了作物品种区域试验多年资料联合分析 方法的计算机程序。程序调试结果表明，两种系统开发的程序不仅可以进行方差分析，而且可以计算品种、试点和 品种试点组合的最小二乘均数，并可对其差异显著性进行检验，且计算结果完全相同。 关键词：作物品种区域试验；计算机程序 中图分类号：-::K;；-L77E: 文献标识码：’ 文章编号：:999M7689（799;）97M99N8M97 作物品种区域试验由于参试品种的参试年份、 参试地点可能不完全相同，甚至还可能出现缺区的 情况，造成试验资料的不平衡性。对于这类资料，明 道绪等根据最小二乘原理提出了相应的分析方 法［:］。但是，至今为止还没有该方法的实际应用。 为了便于运用该方法分析作物品种区域试验多年资 料，对参试品种作出科学评价，本文拟利用 B)3+/# C’-D!6E9 ［;，8］和J.$-’--43,$3B6E :7 ［L］编制该方法的计算机程序。 : 数学模型 作物品种区域试验由于参试年份、参试地区可 能不同，有时可能出现缺区，所得到的资料可能是非 平衡资料。利用明道绪等［:］提出的联合分析方法 分析这类资料时，必须满足观测值个数大于待估计 的独立参数个数，而且作物品种区域试验的目的在 于评价参试品种的生产能力及其适应区域。因此， 可以将方差分析数学模型 !"#$%O!K""K##K$$K（"#）"#K（"$）"$K （#$）#$K（"#$）"#$K%"#%K&"#$% （:） 简化为 !"#$%O!K""K##K$$K（"#）"#K（#$）#$K %%K&"#$% （7） 在（:）和（7）中，!"#$%为第"年份第#地区第$品种第 %区组的观测值；! 为总体平均数；""为第"年份效 应；##为第#地区效应；$$ 为第$ 品种效应；（"#）"# 为第"年份第#地区的互作效应；（"$）"$为第"年份 第$品种的互作效应；（#$）#$为第#地区第$品种的 互作效应；（"#$）"#$为第"年份第#地区第$品种的 互作效应；%"#%为第"年份第#地区内第%区组效应； &"#$%为随机误差、相互独立且都服从 &（9，’7）；%%为 第77卷 第7期 799;年96月 四川农业大学学报 P"+(0/#"&-)*.+/0’%()*+#,+(/#10)2$(3),4 B"#E77Q"E7 P+0E799; !收稿日期：799;M98M78 基金项万方数据 目：国家自然科学基金项目（LN7R9LSR）
第2期 刘永建（等）：作物品种区域试验多年资料联合分析法的计算机程序 127 第区组效应。 表，检验试点间、品种间和试点品种组合间的差异显 著性。 2 Visual BASIC6.0程序编制 显著性检验模块：可对各试点最小二乘均数间 2.1程序编辑 各品种最小二乘均数间和试点品种组合最小二乘均 Visual BASIC是一种可视化的、面向对象和采 数间进行差异显著性检验。 用事件驱动方式的结构化高级程序设计语言。V- 限于篇幅，计草机原程序代码略 al basic6.0是它的一个版本，可以用来建立32 2.2结果输出 位的应用程序。 利用该程序分析作物品种区域试验多年资料， 用Visual BASIC6.0编制作物品种区域试验多 输出的结果有：①方差分析表（见表1）：②品种、试 年资料联合分析方法的计算机程序，主要由数据输 点以及试点品种组合的最小二乘均数：③品种、试点 入模块、正规方程组简缩模块、最小二乘均数估计模 以及试点品种组合最小二乘均数间显著性检验的 块、方差分析模块和显著性检验模块5部分构成（见 值。 图1)。 表1作物品种区域试验多年资料的方差分析表 Table 1 The table of anova for the data of crop variet 效据输入 regional traits in several vears 建立正规方程YX可 变异来源 正规方程的阀绵，得儿码 试当A 求解简缩正规方程。得的估 求的x 品种V 品种组合 乘法，可进子 小二结计 分析 年份X试点YearX Are路 试点×品种Area×Variety 总的Total 3 The SAS System for Windows V6.2程序编制 图1作物品种区域试验多年资料联合分析法的Vil BASIC流程图 3.1程序编辑 The SAs Sytem for Windows V6 12 SAS 数据输入模块：在该模块中，须首先输入年份 STAT中提供了GLM(General Linear Model)过程 其推算参数的理论基出是最小二乘法，最适用于不 数、试点数、品种数、小区数，再输入两个数据文件名 平衡的试验设计。因此，用GLM过程来解决作物 varicty0l.dat和varicty02.dat,其中varicty0l.dh 品种区域试验多年资料联合分析法的计算应该是可 由小区产量构成，无观察值的小区输入“0”：vai 行的，我们编写了分析程序（见图2）。 02dt由“0”"和“1”构成，其中“0”表示相应小风 无观察值，“1”表示相应小区有观察值。在此基础 上，构建正规方程组的结构矩阵X和小区产量列向 anety.dat 量Y。 variety block yied@@ 正规方程组简缩模块：由于正规方程组是不 area variety blok 秩的，本身没有解。但是，根据作物品种区域试验资 od yield=year arca varicty block yearrrvri 料的数学模型和有关限制条件可以对正规方程组进 35H 行行列变换，进而可得简缩方程组。 iety area run: 最小二乘均数估计模块：解简缩方程组，并利用 限制条件，可得试点、品种和试点品种组合的最小 2作物品种区域试多年资料联合分析法的S4S程 乘均数。 for the data gional traits in several v 方差办：在该模块中，可以得到方差分析
第!区组效应。 ! "#$%&’()*+,-./程序编制 !.0 程序编辑 "#$%&’()*+,是一种可视化的、面向对象和采 用事件驱动方式的结构化高级程序设计语言。"#1 $%&’()*+,-./是它的一个版本，可以用来建立2! 位的应用程序。 用"#$%&’()*+,-./编制作物品种区域试验多 年资料联合分析方法的计算机程序，主要由数据输 入模块、正规方程组简缩模块、最小二乘均数估计模 块、方差分析模块和显著性检验模块3部分构成（见 图0）。 图0 作物品种区域试验多年资料联合分析法的"#$%&’()*+,流程图 4#50 678"#$%&’()*+,9’:;:9>78?8>7:@:978@&>&:9C=85#:B&’>=&#>$#B$8E8=&’C8&=$ 数据输入模块：在该模块中，须首先输入年份 数、试点数、品种数、小区数，再输入两个数据文件名 E&=#8>C/0.@&>和 E&=#8>C/!.@&>，其中 E&=#8>C/0.@&> 由小区产量构成，无观察值的小区输入“/”；E&=#1 8>C/!.@&>由“/”和“0”构成，其中“/”表示相应小区 无观察值，“0”表示相应小区有观察值。在此基础 上，构建正规方程组的结构矩阵" 和小区产量列向 量#。 正规方程组简缩模块：由于正规方程组是不满 秩的，本身没有解。但是，根据作物品种区域试验资 料的数学模型和有关限制条件可以对正规方程组进 行行列变换，进而可得简缩方程组。 最小二乘均数估计模块：解简缩方程组，并利用 限制条件，可得试点、品种和试点品种组合的最小二 乘均数。 方差分析模块：在该模块中，可以得到方差分析 表，检验试点间、品种间和试点品种组合间的差异显 著性。 显著性检验模块：可对各试点最小二乘均数间、 各品种最小二乘均数间和试点品种组合最小二乘均 数间进行差异显著性检验。 限于篇幅，计算机源程序代码略。 !.! 结果输出 利用该程序分析作物品种区域试验多年资料， 输出的结果有：!方差分析表（见表0）；"品种、试 点以及试点品种组合的最小二乘均数；#品种、试点 以及试点品种组合最小二乘均数间显著性检验的$ 值。 表0 作物品种区域试验多年资料的方差分析表 6&A’80 678>&A’8:9)FG")9:=>78@&>&:9C =85#:B&’>=&#>$#B$8E8=&’C8&=$ 变异来源 *:%=C 年份K试点 J8&=K)=8& 试点K品种 )=8&K"&=#8>C 误差 L==:= 总的 6:>&’ 2 678*)**C$>8?9:=M#B@:;$"-.0!程序编制 2.0 程序编辑 678*)**C>8?9:=M#B@:;$"-.0!在 *)*／ *6)6中提供了 NOP（N8B8=&’O#B8&=P:@8’）过程， 其推算参数的理论基础是最小二乘法，最适用于不 平衡的试验设计。因此，用 NOP 过程来解决作物 品种区域试验多年资料联合分析法的计算应该是可 行的，我们编写了分析程序（见图!）。 >#>’8‘678,:?D%>8=Q=:5=&?:9678P8>7:@:9,:?A#B8@)B&’C$#$ 9:=678R&>&:9,=:D"&=#8>CS85#:B&’6=#&’$#B*8E8=&’J8&=$’； @&>&"&=#8>C； #B9#’8‘"&=#8>C.@&>’； #BD%>C8&=&=8&E&=#8>CA’:CA’:I； ?:@8’C#8’@TC8&=&=8&E&=#8>CA’:C／$$0 $$!$$2$$U； ’$?8&B$E&=#8>C&=8&E&=#8>C!&=8&／>5#99； =%B； 图! 作物品种区域试验多年资料联合分析法的*)*程序 4#5! 678*)*D=:5=&?:9>78?8>7:@:978@&>&:9C=85#:B&’>=&#>$#B$8E8=&’C8&=$ 第!期 刘永建（等）：作物品种区域试验多年资料联合分析法的计算机程序 0!V 万方数据
第2期 刘水建（等）：作物品种区域试验多年资料联合分析法的计算机程序 3.2结果输出 程序都能正常运行。 利用该程序分析作物品种区域试验多年资料 分析结果中，都得到了方差分析表和品种、试点 输出的结果有：①方差分析表（见表2）：②品种、试 以及试点品种组合的最小二乘均数及其显著性检验 点以及试点品种组合的最小二乘均数及其显著性检 的t值 验的t值。 我们对分析结果进行了比较，两个程序分析结 表2作物品种区域试验多年资料的方差分析表 果完全相同。 Table 2 The a of cro 参文款 [1明道绪，黄 模型Modd Ming D X.Huang Y B.Wang C. 区组ock F份Ycw Vanety YearX Area Ming DX.Liu Y J.Hunng Y B.et a.The method of least [4]潭强.BASIC语言 ”结构化程序设计教程M门.北京：中 4程序调试 [5刘 不SC程序设计教程（第二版汇M.北京：清 分别用编制的Visual BASIC程序和SAS程序 大学出版社，2003】 对2001年和2002年国家玉米区域试验部分资料进 (本文审稿：潘光堂 行分析调试（具体分析结果，我们将另文介绍），两个 万方数掘
!"# 结果输出 利用该程序分析作物品种区域试验多年资料， 输出的结果有：!方差分析表（见表#）；"品种、试 点以及试点品种组合的最小二乘均数及其显著性检 验的!值。 表# 作物品种区域试验多年资料的方差分析表 $%&’(# $)(*%&’(+,-./0-,+1*)(2%*%+,31+4 5%16(*71(86+9%’*1%6*:69:(5(1%’7(%1: 变异来源 ;+<%1(:,+1%9%’7S6982%*%T6*)=6::6985%’<(:69*T+U,%3*+1(V4(1U 6=(9*T6*)1%92+=6S(2&’+3B2(:689%926*:%44’63%*6+9［K］"&’!( &)%*+*,-’(.-+-’(，JLLM，#!（G）：G!#@G!L" ［!］ 彭昭英"世界统计与分析全才;-;系统应用开发指南［?］" 北京：希望电子出版社，#III" ［G］ 谭浩强"A-;HF语言———结构化程序设计教程［?］"北京：中 国科学技术出版社，JLLI" ［W］ 刘炳文"06:<%’A-;HF程序设计教程（第二版）［?］"北京：清 华大学出版社，#II!" （本文审稿：潘光堂） 第#期 刘永建（等）：作物品种区域试验多年资料联合分析法的计算机程序 万方数据
作物品种区域试验多年资料联合分析法的计算机程序 日万方藏据文秋传： 刊名： 四川农业大学学报P回 英文刊名： JOURNAL OF SICHUAN AGRICULTURAL UNIVERSITY 年，卷（期） 2004,22(2) 引用次数： 1次 文献6 1.明道绪。黄玉碧作物品种区城试验多年资料的联合分析方法[期刊论文]-作物学报1997() 乙，明道结，刘水建有缺失数据两因素随机区组试验资料的最小二柔分桥法及其应用[期刊论文]-作物学报1997（④） 3.彭照英世界统计与分析全才5S系统应用开发指南2000 4.谭浴强BASICE语言一结构化程序设计教程1990 5.刘炳文Visua1BAS1C程序设计教程2003 相似文献) 1.期刊论文刘永建，明道，黄玉碧，范志霞作物品种区域试验多年资料联合分析法的应用一四川农业大学学报 2004.22(3) 1.刘水建.明道绪。黄玉碧.范志度作物晶种区城试验多年资料联合分析法的应用[期刊论文]-四川农业大学学报 2004(03) 本文链接：http:/a.g.vanfangdata.com.cn/Periodical.scmydxxb200402008.aspx 下载时间：2010年3月5日
作物品种区域试验多年资料联合分析法的计算机程序 作者： 刘永建， 明道绪， 黄玉碧， 范志霞 作者单位： 四川农业大学,农学院,四川,雅安,625014 刊名： 四川农业大学学报 英文刊名： JOURNAL OF SICHUAN AGRICULTURAL UNIVERSITY 年，卷(期)： 2004，22(2) 引用次数： 1次 参考文献(5条) 1.明道绪.黄玉碧 作物品种区域试验多年资料的联合分析方法[期刊论文]-作物学报 1997(1) 2.明道绪.刘永建 有缺失数据两因素随机区组试验资料的最小二乘分析法及其应用[期刊论文]-作物学报 1997(4) 3.彭昭英 世界统计与分析全才SAS系统应用开发指南 2000 4.谭浩强 BASIC语言—结构化程序设计教程 1990 5.刘炳文 Visual BASIC程序设计教程 2003 相似文献(1条) 1.期刊论文 刘永建.明道绪.黄玉碧.范志霞 作物品种区域试验多年资料联合分析法的应用 -四川农业大学学报 2004,22(3) 用作物品种区域试验多年资料联合分析方法的计算机程序对2001～2002年度国家玉米品种区域试验的部分资料进行分析.结果表明,该程序不仅可以 进行方差分析,而且可以计算品种、试点和品种试点组合的最小二乘均数,并可对其差异显著性进行检验. 引证文献(1条) 1.刘永建.明道绪.黄玉碧.范志霞 作物品种区域试验多年资料联合分析法的应用[期刊论文]-四川农业大学学报 2004(03) 本文链接：http://d.g.wanfangdata.com.cn/Periodical_scnydxxb200402008.aspx 下载时间：2010年3月5日
