不等距重复测量设计方差分析用sas和spss实现的对比

《不等距重复测量设计方差分析用sas和spss实现的对比》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、维普资讯http://www.cqvip.com·352·中国卫生统计2007年8月第24卷第4期不等距重复测量设计方差分析用SAS和SPSS实现的对比*第四军医大学卫生统计学教研室(71o032)李运明曹文君陈长生△【提要】目的探讨不等距重复测量设计方差分析的条件及结果解释。方法用SAS和SPSS软件实现不等距重复测量设计方差分析。结果给出了方差分析条件、SAS和SPSS程序及主要结果解释。结论只有满足方差分析条件、正确合理解释结果，才能保证不等距重复测量设计方差分析的有效性。【关键词】重复测量统计软件方差分析重复测量数据(repeatedmeasuresdata)是

2、医学领域矩阵A十分常见的一种数据形式，是对同一研究对象的某个测量指标或多个测量指标，在不同时间点的多次重复测量结果。该类数据具有系统结构的特点，即数据间存在一定的自相关性和随机误差至少分布于两个层次水平上，I1(71常常不满足常规统计分析方法所要求的独立性假定，使lL11得其分析方法有别于一般的统计分析方法u，2J。重复测矩阵B量设计资料的方差分析是实际工作中和教科书中最常用的方法，但其应用及其结果解释仍存在着不同程度的％一些问题，尤其对于测量时间间隔不相等的资料，即不等距重复测量设计资料【3】。本文拟结合生物医学实例，着重对著名统计软件SASPROCGLM和SPSSG

3、LM进行不等距重复测量设计方差分析的前提条件、过程和结果解释等方面作一比较。2．SAS和SPSS程序前提条件及程序SAS和SPSS软件中可处理重复测量设计的过程有：SASPROCGLM，SASPROCMIXED，SPSSGLM1．重复测量设计方差分析前提条件和SPSSMIXED。本文采用实际工作中最常用的重复测量设计方差分析有：一元方差分析和多元SASPROCGLM和SPSSGLM进行不等距重复测量方差分析。重复测量设计一元方差分析，除假定组间设计方差分析对比，其程序见表1。标准正交对比协方差阵(covariancematrix)齐同(齐性条件)外，还需假定共同的标准正

4、交对比协方差阵为实例分析j，，即球对称条件【引。对球对称条件，通常可假定原一项药物代谢动力学研究，目的是对比某种药物始资料的协方差矩阵满足复合对称性【5(compound的不同剂型在体内的代谢速度。剂型分胶囊型和片剂symmetry)，即各测量时间点上的测量值的方差相同型。将16名受试对象随机分为两组，每组8名。一组(用2表示)，任何两测量时间点上测量值的协方差也给予胶囊，另一组给予片剂，分别在服药后1小时、2相同(用(71表示)，如矩阵A的形式。重复测量设计小时、4小时、6小时及8小时测定血中的药物浓度。多元方差分析与传统的多元方差分析对资料的假定相测定结果见表2。S

5、AS和SPSS主要结果对比及解释同[63，即：①反应变量服从多变量正态分布；②组间方如下：差一协方差矩阵齐同，是否满足此假定可用Box’SM1．Mauchly’S球对称检验齐性检验。重复测量设计多元方差分析通常不对协方SAS给出了两种变换后的球对称检验结果(表差矩阵的结构做出具体假定，其结构可为一般形式，如3)：①是对原始数据按照REPEATED语句中设定的矩阵B的形式。对比矩阵进行变换后的检验结果；②是按照标准正交对比变换后的检验结果。本文中REPEATED语句中*基金项目：国家自然科学基金资助项目(39900126)陕西省自然科学设定的对比矩阵就是POLYNOMIA

6、L【3】，所以①与②基金资助项目(2003Fl1)△通讯作者：陈长生，Email：chencs(~fllllnt1．edu．cn的检验结果一致。SPSS仅给出了进行标准正交对比维普资讯http://www.cqvip.comChineseJournalofHealthStatistics。Aug2007．Vo1．24．No．4．‘353’变换后的检验结果，同时给出了校正系数￡的三种估协方差阵不满足球对称性，因此一元方差分析时需对计值(表4)[6j。本例P<0．01，共同的标准正交对比组内效应的自由度用校正系数￡进行校正。表1SAS／SPSS软件处理不等距重复测量设计资料

7、的程序SAS程序SPSS程序DATAexample：GLMDOtreat=1TOm；／*m为处理分组数*／t：imel⋯timekBYtreatDOid=1TOn；／*n为每组的观测单位数*／／rWSFA(OR=timekPolynomial(tl⋯tk)INPuTtimel⋯timek@@；／*k为重复测量次数*／／METHOD=SrYPE(3)OUTPUT；END；END；CARDS；／rpLOT=PROFILE(time*treat)数据行／PRINT=ETASQHOMoGENEITY；RUN；／CRITERIA=ALPHA(．05