随机区组设计方差分析概述

随机区组设计又称为配伍设计,该方法属于两因素方差分析(Two-Way ANOVA),用于多个样本均数间的比较,比如动物按体重、窝别等性质配伍,然后随机地分配到各个处理组中,即保证每一个区组内的观察对象的特征尽可能相近。同一受试对象在不同时间点上观察,或同一样品分成多份,每一份给予不同处理的比较也可用随机区组设计进行分析。

随机区组设计分组原则:在某些研究中,先将受试对象按可能影响试验结果的属性分组(非随机组),分组的原则是将属性相同或相近的受试对象分在同一组内,如将病人按年龄/性别/职业或病情分组,或者将动物按性别/体重分组,然后采取随机化的方法对每个组内的受试对象分配各种处理。如此以来,可使得区组内的观察单位同质性好,使各比较组的可比性强,使组间均衡性好,处理因素的效应更容易检测处理。

随机区组设计方差分析用于分析两个或两个以上因素是否对不同水平下样本的均值产生显著的影响;检验多个因素取值水平的不同组合之间,因变量的均值是否存在显著性差异。其既可以分析单个因素的作用(主效应),也可以分析因素之间的交互作用(交互效应),还可以进行协方差分析,以及各因素变量与协变量之间的交互作用。

若有两个因素A与B,因素A与B间不存在交互作用,那么可以对因素A和B各自进行独立分析,在后续分析中去除不显著的因素。如果方差分析结果显示因素A和B间存在交互作用,则需对数据进行进一步分析,具体包括:

SPSS实现随机区组设计方差分析

示例:研究3种不同的避孕药A/B/C在体内的半衰期,考虑到窝别对结果的影响,采用随机区组设计方案。将同一窝别的3只雌性大白鼠随机分配到A/B/C 3组,测定该药在血液中的半衰期(小时),试分析3种药物的半衰期有无不同?

独立样本检验中的假设检验汇总(连续变量假设检验)(1)

1. 示例分析:

2. 数据录入:

独立样本检验中的假设检验汇总(连续变量假设检验)(2)

3. 建立假设:建立检验假设,确定检验水准 α

4. 随机区组设计方差分析简要

(1) 打开 分析—一般线性模型—单变量

独立样本检验中的假设检验汇总(连续变量假设检验)(3)

(2) 参数说明

独立样本检验中的假设检验汇总(连续变量假设检验)(4)

a. 因变量:试验研究的结果变量,为数值型

b. 固定因子:指因素所有水平在样本中都出现了,一般为分类变量,如血型、不同治疗方案等

c. 随机因子:指因素所有可能的水平在样本中没有都出现,多为连续性变量,如研究不同温度(30°/40°/50°)对某药物的影响。

d. 协变量:对因变量可能有影响,需要在分析时对其作用加以控制的连续性变量。

e. 权重:用于加权的最小平方分析,权重变量可用于样本赋以不同的权重。

f. 模型:

独立样本检验中的假设检验汇总(连续变量假设检验)(5)

独立样本检验中的假设检验汇总(连续变量假设检验)(6)

g. 对比:对比用来检验因子的水平之间的差值

独立样本检验中的假设检验汇总(连续变量假设检验)(7)

h. 图:可绘制一个或多个因素变量作用后的因变量的均值分布图

独立样本检验中的假设检验汇总(连续变量假设检验)(8)

i. 事后比较:可对因素的各个水平进行两两比较,同前面讲述的方差分析中两两比较,具体含义在此不再阐述。拒绝无效假设后,需具体判断哪些均属不同就要做多重比较,只有水平数大于3才可选择

独立样本检验中的假设检验汇总(连续变量假设检验)(9)

j. EM平均值:

独立样本检验中的假设检验汇总(连续变量假设检验)(10)

k. 保存 与 选项

独立样本检验中的假设检验汇总(连续变量假设检验)(11)

5. 数据结果与说明

(1) 输出处理因素(药物)及区组因素的取值和样本量。

独立样本检验中的假设检验汇总(连续变量假设检验)(12)

(2) 方差分析结果:从下表可看出,drug的F=7.204,P=0.016<0.05,拒绝H0,接受H1,认为三种避孕药的半衰期不完全相同;区组因素block的P=0.174>0.05,认为区组因素没有统计学意义;药物因素drug的p=0.016<0.05,认为不同药物之间有统计学意义。

独立样本检验中的假设检验汇总(连续变量假设检验)(13)

(3) 多重比较结果:从LSD中可看出,药物a与药物c之间存在显著性差异。

独立样本检验中的假设检验汇总(连续变量假设检验)(14)

(4)下图显示了不同区组不同药物的半衰期的均值图,可见药物a半衰期最短,药物c半衰期最长。

独立样本检验中的假设检验汇总(连续变量假设检验)(15)

  1. 语法

UNIANOVA x BY drug block /METHOD=SSTYPE(3) /INTERCEPT=INCLUDE /POSTHOC=drug block(TUKEY LSD BONFERRONI) /PLOT=PROFILE(block*drug) TYPE=LINE ERRORBAR=NO MEANREFERENCE=NO YAXIS=AUTO /CRITERIA=ALPHA(.05) /DESIGN=drug block.

,