4.4 拉丁方设计一元定量资料方差分析

4.4.1 问题与数据

【例4-4】某研究者欲比较不同浓度的氯化钠溶液（1%、2%、4%、8%、16%，分别以数字1、2、3、4、5代替）静脉注射对家兔血压的影响，选取5只体重2～3kg的家兔进行实验，每只家兔分别注射不同浓度的氯化钠溶液各1次，注射量以2ml/kg计算，每次注射10s，观察并记录每次注射前后血压的最大变化值。两次注射之间间隔一段时间，使家兔血压恢复至初始水平，实验结果如表4-4所示。已知因素间交互作用可以忽略不计，请判定资料设计类型，并进行合适的统计分析。

表4-4 家兔静脉注射5种不同浓度NaCl溶液后血压升高值

4.4.2 对数据结构的分析

在本资料中，受试对象为5只家兔，实验因素为NaCl溶液浓度，行区组因素为兔号，列区组因素为注射次序，观测指标为血压升高值，故此资料为拉丁方设计一元定量资料。

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

此资料为拉丁方设计一元定量资料，应选用其相应的方差分析进行处理。

4.4.4 SAS程序

SAS程序名为SASTJFX4_4.SAS。

程序中第1步是建立数据集，rabbit代表“家兔编号”，order代表“注射次序”，dose代表“不同浓度的氯化钠溶液”，y代表“血压升高值”。第2步是调用ANOVA过程进行方差分析，并对5种浓度的氯化钠溶液的升血压效果进行两两比较。

4.4.5 主要分析结果及解释

这是输出结果的第1部分。可以看出，不同家兔和不同注射顺序对血压升高值影响之间的差异均无统计学意义，而不同浓度的氯化钠溶液对血压升高值影响之间的差异则有统计学意义（F=478.26，P<0.0001）。

Student-Newman-Keuls Test for y

这是结果输出的第2部分。由第1部分结果可知，不同浓度的氯化钠溶液对血压升高值的影响不全相同，这里给出了5种浓度的氯化钠溶液对血压升高值影响两两比较的结果。由“SNK Grouping”列可以看出，5种浓度的氯化钠溶液对血压升高值影响之间的差异均有统计学意义。由均值大小可知各浓度氯化钠溶液对血压升高值影响的程度由大到小分别为16%>8%>4%>2%>1%。

其实，从第1部分方差分析的结果可以看出，不同家兔和不同注射顺序对血压升高的影响之间的差异没有统计学意义，可在SAS程序中的class和model语句中删除rabbit和order，这样可以增大误差项的自由度，使结果更加稳定可靠。重新运行程序，得检验结果如下：

结果与前面相同，专业结论同上。