SPSS学习笔记之——重复测量的多因素方差分析

1、SPSS学习笔记之一一重复测量的多因素方差分析（2012-08-14 22:30:39）转载标签：spss方差分析 统计 杂 分类：SPSS谈1、概述重复测量数据的方差分析是对同一因变量进行重复测量的 一种试验设计技术。在给予一种或多种处理后，分别在不同 的时间点上通过重复测量同一个受试对象获得的指标的观 察值，或者是通过重复测量同一个个体的不同部位（或组织）获得的指标的观察值。重复测量数据在科学研究中十分常见。分析前要对重复测量数据之间是否存在相关性进行球形检验。如果该检验结果为 P 0.05，则说明重复测量数据之间 不存在相关性，测量数据符合Huynh-Feldt条件，可以用单因素方差分析

2、的方法来处理；如果检验结果P 0.05，则说明重复测量数据之间是存在相关性的，所以不能用单因素方 差分析的方法处理数据。在科研实际中的重复测量设计资料 后者较多，应该使用重复测量设计的方差分析模型。球形条件不满足时常有两种方法可供选择：（1）采用MANOVA（多变量方差分析方法）；（2）对重复测量ANOV检验结果中与时间有关的F值的自由度进行调整。2、问题新生儿胎粪吸入综合征（MAS是由于胎儿在子宫内或着生 产时吸入了混有胎粪的羊水，从而导致呼吸道和肺泡发生机 械性阻塞，并伴有肺泡表面活性物质失活，而且肺组织也会 发生化学性炎症，胎儿出生后出现的以呼吸窘迫为主，同时 伴有其他脏器受损现象的一组

3、综合征11。血管内皮生长因子（vascular endothelial growth factor，VEGF）是一种有丝分裂原，它特异作用于血管内皮细胞时，能够调节血管内 皮细胞的增殖和迁移，从而使血管通透性增加。而本实验旨 在通过观察分析给予外源性肺表面活性物质治疗前后胎粪 吸入综合征患儿血清中 VEGF勺含量变化，评价药物治疗的效 果。将收治的诊断胎粪吸入综合症的新生儿共42名。将患儿随机分为肺表面活性物质治疗组 （PS组）和常规治疗组（对照组）， 每组各21例。PS组和对照组两组所有患儿均给予除用药外的 其他相应的对症治疗。PS组患儿给予牛肺表面活性剂 PS70mg/kg治疗。采集PS组

4、及对照组患儿0小时，治疗后24小时 和72小时静脉血2ml，离心并提取上清液后保存备用并记录 血清中VEGF勺含量变化情况。结果如下：* 3 1血借中YEGF的劇冶哉据gMX1.09 1 7720 094OlSM 1 314 1 5160 987 1 03 0 8360 7460 55(5 1.12 11.211 0.983 0.934 O.S28 1.141594 l.OsMf24 W1 703 L54S 1.23 1 135 1A22 1.123 1.345 1.D9S L235 h33fi 1.3C7 1.241.214 1 234 1 313 1 464 1 二1 1.356 1 29

5、8 1434 0J76raw1.703 L54S 1.23 1 435 1 422 1 823 1 345 1 412 1 435 1.336 1367 1213541J341.413 1 4641,271 1 356119& 1 334 1423绢(M册0,88 L03 1.2261.55S0 S54 1 030.778 L030 856078 0.868 L030,640.E 1,272 1 DOO I 0S0 8621.050 SS42Ab时1.203 L081.19 1J351J22 L123221 0.P87 1 235 1 236 L432 LW1.214 0.934 1.213

6、L164 1.21 0 556 1298 1 134 0 97672刈耐L20? L368 1.423 1.335 1322 LI23 1.4S3 1387 1.235 L236 L432 111214 1 334 L213 1 364 L21 0,956 1 2PS 1.134 0-9763、统计分析建立数据文件变量视图：i名称粪型|小数标签II值ngroup数值80分组1,对照组Lt1数 (N)83(P卜时无t2数值7324小时无t3feH(N)83茂小时无数据视图:菜单选择:首先进入如下对话框，在“被试内因子名称”中输入“time“级别数”输入3，因为每个患者重复测量了 3次。后点击“添

7、加”按钮。此时下方“定义”按钮变为可用，点击进入下列对话框：将“ group ”选入“因子列表”框，t1-t3分别选入“全体内变量（time ） ”框内，如下图所示：点击右上角“模型”按钮，进入以下对话框，选择“设定” 将“ time ”选入“全体内模型”框，“ group ”选入“群体 间模型”框，“构建项”选择“主效应”。下方的平方和选“类型III，这是对于平衡数据。如果两组样本量不等， 则选择“类型IV重复度量：模璽itsflla0全因子()發定! time构建项|time群低间回1主沁群住间環型卫）叵Il/ij. groupgroup主体内辿）：群伟内視型廻）:平方和Q）点击“继续”返

8、回，点击“绘制”按钮。进入下面对话框:将“ time ”选入“水平轴”，group选入“单图”，然后点 击“添加”按钮，下面框中会显示“time*group ”。点击“继续”返回，点击“两两比较”按钮，将group选入右侧“两两比较检验”框中，选中复选框“LSD”。点击“继续”返回，点击“选项”按钮，进入下面对话框: 将time选入“显示均值框”，选中“比较主效应”复选框， 选中下方“描述统计”复选框。下方显著性水平设为0.05。点击“继续”返回，点击“确定” 输出结果。4、结果解读:描迷性纯计量分组均值标堆偏差N0小时对IS趙.97167.35672021治疗爼97924.20085721总

9、计97545.2859504224小时对黑组f.2B49516373321洽疗组1 15076.14153721总计1 22106.16409842吃小时对照組1 41171.13744121治疗组1.25567.1 3868121总计1.33369.15753342这是一个关于各个时间点的两组数据描述性统计。主体内姑風MauchlY 的 W近似初dfStgEpsilon1renhouse- G 歸斫HuynbF 创 dttima14拥gJ.001763 1807M谢的暮JEflE益號卜度量MEASURE栓址債锻盃即恆灌正交暮换田实鼻的谓菱忖方差矩阵与一丨睾住味降应比例.a 6用于進整豐署愷+

10、均植桂的冒由庫.隹”主馬内SM鶴1T古屆中豐十修正肓的检竝.1b 设计:JSE * group 主体内您计bma这是球形检验结果，p=0.0010.05，所以不满足球形分布假设，需要进行多变量方差分析或者自由度调整，SPSS接下来会给出以上两种结果。效应值F假设df误差dfStimeFilQi肋跟踪.57941.29712.00039.000Wilks 的 Lambda.321141 297J2.DOO39 000H血III叩的跟跨2.11841.297J2.00039 000冃眄的量大根2.11841 297J2.00039 000time * groupPillai 1?1,05 inI.

11、943-32.00039 000Wilks 的 Lambda.909V94Ra2 00039 000H血III陶的跟備1001.94BJ2.ODO39 000.1001.94BJ2.00039.000怎精确毛直计量b.设计：戳距+ grojp 主体內设计:time这是进行多变量方差分析的结果，给出了4种统计量，它们的检验结果一致，time的P0.05，说明时间和分 组无交互作用，说明时间因素（即0小时、24小时、72小时） 的作用不随分组（即治疗组和对照组）的不同而不同。JTSMEASURE 1川型平方和df均力Flime采用的再常度2B222 I1.41142 720Greenhouse-G

12、eisset2.B221 5251.94942.720HuynKFeldt2 8221 6141748X2 720TH2.B221 COO282242720lime-group 采用科瑋洋虜J 5920B02414Greenhouse-GeisserJ 591 52S1042 414Huynh-F$ltl|.1591 6U0992 4T4下限J591 000.1592.414俣差(time)采用的球形庫16+J80Greenhouse-Gisser2 64 2B1 054.043Hirynb-Flcft2 64 264 571041Tffi2.64540000D06所谓“主体内”，即是重复测量

13、的各个时间点。上表是用各 个时间点进行分组的方差分析表，给出4种统计量，第一种为满足球星假设的情况，后三种对自由度进行了校正，本题 目中不满足球形分布假设，只能看下面的三种检验方法。结 果解释同上一个表。主体间效应的检验度量 MEASURE 1 转换的变量呼均值:源川型平方和df均方FSig叢距174.5511174.512898.326.000group.2641.2644.384.043饯差2.40940,060这是对分组的方差分析，对变量进行如下的变换：y=(t1+t2+t3)/sqrt(3)。P=0.0430.05，说明有治疗组与对照组之间有统计学差异。IMEASURED岸Jf占第小師

14、均伯Atimeill.w辿肾伉这个图可以直观地看出测量指标随时间的变化趋势。治疗组与对照组两组资料随时间变化的趋势大致相同，治疗组血清 中VEG的含量较对照组呈下降趋势，说明治疗组的效果优于 对照组。我们还可以给出在每个时间点上两个分组之间的比较，需要 用到多变量方差分析：操作步骤如下：跟之前操作类似，不 赘述，看图就行參变量凶指定楼型主因子徃）越设定-bHIIMM 1 hlMMM&hiiri1曲刊我泾屋l(u)-世9叫group构读顶Sffi(P)：L国平方盹）慣S1川 匸4在植型中程含聲结果输出丰佛间效应的橙验源因变量川型平方和df均方FSig.校正欖型口小时.001 J100100793324佃J67fc1.1677 139011256c125613.4