单因素重复测量方差分析原理和SAS代码实现

发布时间:2024年01月19日

重复测量设计:就是对研究对象某指标进行重复测量,最常见的情况是在前、后分别测量1次,又称为前后测量;当测量次数>=3时,即为重复测量设计。一般采取的统计学方法是重复测量方差分析。

重复测量设计包括单组重复测量设计和多组重复测量设计。

应用条件:要求因变量符合

  1. 正态性
  2. 方差齐性(针对于多组的情况,因为单组只有一个方差);
  3. 另外还需满足球形性,否则需要校正;不满足球形性,可采用以下两种方法:
  1. G-G校正和H-F校正;
  2. 使用多变量方差分析结果

这里建议大家可以直接使用多变量方差分析的结果。

  1. 单组重复测量设计

设计形式如下:没有分组,只是对每个人进行重复测量多次

ID

T0

T1

T2

T3

1

2

3

4

案例:某研究者为研究高强度间歇运动对收缩压变化的影响,招募了5名参与者,要求在参与者每天进行10分钟的高强度间歇运动,分别在基线、第3周、第7周测量血压。

具体如下:

ID

T0

T3

T7

1

130

125

124

2

120

119

117

3

140

130

120

4

110

110

90

5

100

100

99

/*首先写入数据*/

data?dan;

?do?id=1?to?5;

??input?t0 t3?t7@@;

??output;

?end;

?cards;

130 125 124

120 119 117

140 130 120

110 110 90

100 100 99

;

run;

proc?print;quit;

/*单组重复测量方差分析程序*/

proc?glm?data=dan;

?model?t0 t3 t7=/nouni;

/*model语句指定3个因变量,即 t0 t5 t10,

?由于是单组的,因此不存在自变量,自变量也就不指定,

nouni表示不输出单变量分析结果,因为对于重复测量分析,单变量分析意义不大*/

repeated?time 3(0?3?7) contrast(1)/printe?summary;

/*repeated语句指定有3个时间点,(0 3 7)表示写上时间;如果不写,默认的时间间隔是等距的,而此例是非等距的

contrast(1)表示后2个时间点分别与第1时间点比较,

printe输出球性检验,

summary输出各时间点与对照时间点的比较结果*/

run;

1.下面是基本描述

2.下面给出了球形检验和多变量方差分析结果

3.下面给出的是单变量方差分析结果

结果解释:

这里由于满足球形性,因此选择未调整的p=0.0491<0.05,说明血压时间变化有统计学意义。

4.最后给出后两个时间点分别与对照时间点比较的结果。

/*如果我们想进一步了解存在什么样的趋势,*/

/*可以通过在repeated语句中指定polynomial来实现。*/

proc?glm?data=dan;

?model?t0 t3 t7=/nouni;

repeated?time 3(0?3?7) polynomial/printe?summary;

/* polynomial来实现趋势的检验*/

run;

从结果来看,时间的1次项和2次项均无意义。

今天的分享先到这里,下期分享多组重复测量方差分析

文章来源:https://blog.csdn.net/kaiming000/article/details/135686173
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。