4.6 析因设计一元定量资料方差分析

4.6.1 问题与数据

【例4-6】某研究者欲研究IL-11药对5.5Gy照射小鼠骨髓造血细胞周期（G0/G1期）的影响，选取45只小鼠并将其完全随机地均分成3组，每组15只，分别在5.5Gy剂量照射前给IL-11药、照射后给IL-11药和照射对照（即不给IL-11药）。每一组中的15只小鼠随机地等分成3组，分别在照后6h、12h、24h三个时间点上处死，测量其骨髓造血细胞周期（G0/G1期）。实验结果如表4-6所示，请进行相应的统计分析。

表4-6 IL-11药物对5.5Gy照射小鼠骨髓造血细胞周期（G0/G1期）的影响

4.6.2 对数据结构的分析

实验中涉及两个实验因素：“给药时机”和“处死时间”，受试对象为小鼠。从表面上看，此设计像具有一个重复测量的两因素设计，这样看是把“处死时间”当成了重复测量因素。其实，在不同时间点上获得的数据没有关联，不是来自于同一只小鼠，因为每只小鼠身上只能测到一个数据。因此，此资料不是具有一个重复测量的两因素设计定量资料。由于实验条件为两因素各水平的全面组合，研究者在每种实验条件下均进行了5次独立重复实验，在没有专业依据认为这两个实验因素对观测结果的影响重要程度不同的情况下，此资料应为两因素析因设计定量资料。

4.6.3 分析目的与统计分析方法的选择

对本资料，应采用两因素析因设计定量资料方差分析来处理。这样既可分析两实验因素的主效应，也可分析二者之间的交互效应。

4.6.4 SAS程序

SAS程序名为SASTJFX4_6.SAS。

程序第1步为建立数据集，treat代表“给药时机”，rep代表“受试对象个体”，time代表“处死时间”，outcome代表“骨髓造血细胞周期（G0/G1期）（单位）”。第2步为进行两因素析因设计定量资料方差分析，其中“treat|time”等价于“treat time treat*time”，TDIFF和PDIFF则分别用来输出两两比较的t值和相应的P值。

4.6.5 主要分析结果及解释

这是输出结果的第1部分。由方差分析结果可以看出，两个实验因素“给药时机”、“处死时间”及其交互作用均有统计学意义，即不同给药时机、不同处死时间上骨髓造血细胞周期的值是不完全相同的。

Least Squares Means

这是输出结果的第2部分。本部分首先给出了各实验组在两两比较中的编号，见“LSMEAN Number”列。随后给出了各实验组两两比较的结果，行与列的交叉处即是这两组比较的结果，上行为t值，下行为对应的P值，形式上看两两比较用的是t检验，实际上是方差分析计算出来的。另外，由于篇幅排版限制，对输出结果进行了四舍五入，所以本部分结果比SAS实际输出结果保留的小数点后位数要少一些。由两两比较结果可知：除实验组1与实验组7、实验组2与实验组8、实验组3与实验组6、实验组3与实验组9、实验组5与实验组8、实验组6与实验组9之间观测指标的差异无统计学意义外，其他各组间观测指标的差异均有统计学意义。