生物统计学第九章Word文档格式.docx

1、 class temp time; model effect=temp time; means temp time/duncan alpha=0.05; The SAS System Analysis of Variance Procedure Class Level Information Class Levels Values TEMP 3 1 2 3 TIME 3 1 2 3 Number of observations in data set = 9Dependent Variable: EFFECT Sum of MeanSource DF Squares Square F Valu

2、e Pr FModel 4 87.0707778 21.7676944 12.56 0.0155Error 4 6.9321778 1.7330444Corrected Total 8 94.0029556 R-Square C.V. Root MSE EFFECT Mean 0.926256 4.741881 1.31645 27.7622Source DF Anova SS Mean Square F Value Pr TEMP 2 78.7202889 39.3601444 22.71 0.0066TIME 2 8.3504889 4.1752444 2.41 0.2058 Duncan

3、s Multiple Range Test for variable: NOTE: This test controls the type I comparisonwise error rate, not the experimentwise error rate Alpha= 0.05 df= 4 MSE= 1.733044 Number of Means 2 3 Critical Range 2.984 3.050 Means with the same letter are not significantly different. Duncan Grouping Mean N TEMP

4、A 31.477 3 3 B 27.570 3 2 C 24.240 3 1 Duncan Grouping Mean N TIME A 28.447 3 3 A A 28.440 3 2 A 26.400 3 1从方差分析结果可以得知，温度是极显著的影响因素，时间是不显著因素。在Duncan检验中，温度的三个水平之间差异是显著的。时间的三个水平间差异不显著。本实验是二因素固定模型设计，如果设置重复，会得到两个因素之间的交互作用（如果存在的话），其结果能更好地说明问题。以上方差分析的结果可以归纳成下表：变差来源平方和自由度均方FP温度（temp）78.720 288 9239.360 144

5、422.710.006 6时间（time）8.350 488 94.175 244 42.410.205 8误差6.932 177 841.733 044 4总和94.002 955 689.2 研究浙江蜡梅大苗移栽技术，处理方式包括移栽后的不同覆盖方式和做床方法，统计每100株移栽苗的成活率，结果见下表53：做床方法精细作床仅挖穴覆盖方法遮 阴9385未遮阴9081根据以往经验在覆盖方法与作床方法之间不存在交互作用，对上述结果做方差分析。请注意，这里的结果是百分数。本例需对数据做反正弦变换，程序和结果如下：data plum; do cover=1 to 2; do seedbed=1 to

6、 2; input y ; surrate=arsin（sqrt（y/100）*180/3.14159265;93 8590 81 class cover seedbed; model surrate=cover seedbed; COVER 2 1 2 SEEDBED 2 1 2 Number of observations in data set = 4 SURRATEModel 2 64.5953445 32.2976722 90390.15 0.0024Error 1 0.0003573 0.0003573Corrected Total 3 64.5957018 R-Square C.

7、V. Root MSE SURRATE Mean 0.999994 0.027238 0.01890 69.3987COVER 1 9.4515512 9.4515512 26451.66 0.0039SEEDBED 1 55.1437933 55.1437933 99999.99 0.0001从结果可以看出，覆盖方式和做床方式都是极显著因素。以上结果可以归纳成下表：覆盖方法（cover）9.451 551 2126 451.660.003 9做床方法（seedbed）55.143 793 399 999.990.000 1误 差0.000 357 3总 和64.595 701 839.3 为

8、了研究不同NaCl质量浓度对小麦愈伤组织生长的影响。 配制质量浓度分别为0、0.1、0.3和0.5的NaCl MS培养基，接种15天后，测定每块愈伤组织平均增重百分率，结果见下表54：材料名称NaCl质量浓度 / H870634/%G8901/%极早熟/%中国春/%103.8063.7067.3267.100.199.3156.2752.2452.300.352.2645.0124.1734.300.518.3815.3720.4013.22对上述结果进行方差分析。作者已经给出四种实验材料都是盐敏感型小麦，但是不同的小麦品种必定对盐的抗性不同。也就是说，品种与盐浓度之间存在交互作用，更完善的实

9、验应当怎样设计？这是一个固定模型设计，程序不再给出，结果如下。 CONCEN 4 1 2 3 4 MATERIAL 4 1 2 3 4 Number of observations in data set = 16 INCREASEModel 6 10411.9081 1735.3180 14.18 0.0004Error 9 1101.6046 122.4005Corrected Total 15 11513.5126 R-Square C.V. Root MSE INCREASE Mean 0.904321 22.54545 11.0635 49.0719CONCEN 3 8374.158

10、67 2791.38622 22.81 0.0002MATERIAL 3 2037.74942 679.24981 5.55 0.0196方差分析结果指出，盐浓度是极显著的影响因素，不同基因型的愈伤组织是显著影响因素。根据以往的研究工作的经验，盐浓度与基因型之间可能存在交互作用，最理想的设计应当设置重复，从总平方和中分离出交互作用平方和，问题可以说明得更确切。以上结果可以归纳成下表。浓度间8 374.158 72 791.386 222.810.000 2品系间2 037.749 4679.249 85.550.019 61 101.604 69122.400 511 513.512 79.4

11、 为了研究植物的光合作用，设计了一个实验。将烟草的两个变种种植在田间，利用CO2-depletion技术检测单位叶面积捕获CO2的比率。实验共涉及两个因素：一个是变种，选用了两个变种；另一个是抽样时期，在整个生长季共进行40次田间抽样。这是一个无重复两因素实验设计，方差分析表如下55：抽样时期3.356390.060 433.25*变 种0.015 78.66*0.071 10.001 822.44379 注：*P TREAT 2 26.5151696 13.2575848 1.57 0.2233TYPE 1 3.2191227 3.2191227 0.38 0.5409TREAT*TYPE

12、2 40.3929921 20.1964961 2.40 0.1075Source DF Type III SS Mean Square F Value Pr TREAT 2 10.5603197 5.2801598 0.63 0.5408TYPE 1 0.0089915 0.0089915 0.00 0.9741在方差分析表中我们选用I型可估函数，从F的显著性概率可以得出，不论是养生处理、小鼠类型还是两者的交互作用都是不显著因素。上述结果可以归纳成下表：处理间26.515 169 613.257 584 81.570.223 3类型间3.219 122 70.380.540 9处理类型40.

13、392 992 1 20.196 496 12.400.107 5261.061 904 8318.421 351 8331.189 189 2369.6 野生型C57BL/6及STAT-1-/-型小鼠胰岛，在移入自发糖尿病的NOD#小鼠中之后的存活天数见下表56：养生处理未处理CsA*0 0 2 5 511 11 12 1313 15 17 0 5 8 12 12 158 8 8 10 10 11 185 10 11 11 12 16 206 10 10 1310 125 13 1410 11 11 1212 13# NOD：nonobese diabetic（非肥胖糖尿病）。 * IL-1

14、ra： * CsA：cyclosporine A（环孢菌素 A）。本题与第5题的程序基本一样，下面只给出计算的结果。Model 7 81.7459473 11.6779925 0.55 0.7939Error 39 833.4880952 21.3714896Corrected Total 46 915.2340426 0.089317 45.64659 4.62293 10.1277TREAT 3 67.9880108 22.6626703 1.06 0.3770TYPE 1 8.0410256 8.0410256 0.38 0.5432TREAT*TYPE 3 5.7169109 1.90

15、56370 0.09 0.9656TREAT 3 60.4851648 20.1617216 0.94 0.4290TYPE 1 8.0762698 8.0762698 0.38 0.5423本题的两个主效应和它们的交互作用都是不显著因素。67.988 010 822.662 670 31.060.337 08.041 025 60.543 25.716 910 9 1.905 637 00.090.965 6833.488 095 221.371 489 6915.234 042 6469.7 一项音乐心理学研究，实验是这样设计的：为了避免熟悉的音乐环境，实验安排在两种非典型的音乐练习和演出

16、环境中进行。一种环境是在剧场底层敞开的大厅中（环境A），另一种是在办公室中（环境B）。要求实验参与者学习并回忆所学习的练习曲。学习和回忆包括在相同环境中（AA，BB）和不同环境中（AB，BA），评判学习和回忆的得分，从而判断得分与环境之间的关系57。该实验是一个典型的两因素交叉分组实验设计，方差分析表如下：学习环境180.2671.4410.275回忆环境640.2675.1200.064学习环境1 008.6008.0650.030750.3336125.0562 579.467问：（1）本实验共有几次重复？（2）本实验属于哪一种模型？（3）本实验的两个因素中哪些因素是显著因素？在本实验中显著因素的意义是什么？你可以得到什么结论？（1） 因为本实验共有4种条件组合，df误差（组合1重复数1）（组合2重复数1）（ 组合3重复数1）（ 组合4重复数1）重复数46。因此，重复数6410。 （2） 属固定模型。因为：根据作者所用的检验统计量，由推断，环境的水平是人为选定的。 （3） 只有交互作用是显著的。说明音乐的学习是与环境的两种特定水平有关的。结论：音乐的学习属本例的两种特定环境依赖型记忆。作者的结论是：音乐的学习属环境依赖型记忆。这样的叙述不够严格，只有随机模型才能够得到这样的结论，固定模型只能说“音乐的学习属本例的两种特定环境依赖型记忆”。9.8 与上一实验类似，这次