《重复测量资料方差分析.ppt》由会员分享,可在线阅读,更多相关《重复测量资料方差分析.ppt(74页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、中国疾病预防控制中心中国疾病预防控制中心重复测量设计方差分析重复测量设计方差分析中国疾病预防控制中心中国疾病预防控制中心方方差差分分析析单因素:完全随机设计单因素:完全随机设计两因素:随机区组设计两因素:随机区组设计多因素:析因、拉丁方、正交设计多因素:析因、拉丁方、正交设计重复测量设计重复测量设计单组单组多组多组中国疾病预防控制中心中国疾病预防控制中心第一节 重复测量资料数据特征第二节 重复测量资料的方差分析第三节 重复测量方差分析案例及SAS实现3中国疾病预防控制中心中国疾病预防控制中心第一节 重复测量资料数据特征中国疾病预防控制中心中国疾病预防控制中心重复测量的定义重复测量的定义 重重复
2、复测测量量(repeated measure)是是指指对对同同一一研研究究对对象象的的某某一一观观察察指指标标在在不不同同场场合合(occasion,如时间点)进行的多次测量。如时间点)进行的多次测量。中国疾病预防控制中心中国疾病预防控制中心(1)随)随时间变时间变化化趋势趋势,如生,如生长发长发育模型育模型(2)时间上的稳定性)时间上的稳定性 分析,如环境监测中对同一地区分析,如环境监测中对同一地区 不同时间点的污染情况进行动态监测不同时间点的污染情况进行动态监测(3)结果变量的可重复性,如实验室中某仪器或实验的)结果变量的可重复性,如实验室中某仪器或实验的 准确性准确性 重复测量资料的研究
3、目的重复测量资料的研究目的中国疾病预防控制中心中国疾病预防控制中心重复测量重复测量设计设计重复测量重复测量设计设计(repeated measurement design)受试者内设计(受试者内设计(within-subject design)是指同一观察对象的某观察指标在不同时间点上进行多次观是指同一观察对象的某观察指标在不同时间点上进行多次观察,医学研究中十分常见。察,医学研究中十分常见。如:儿童生长发育中多次测量如:儿童生长发育中多次测量 临床试验中患者接受治疗后在不同时间测量结果临床试验中患者接受治疗后在不同时间测量结果 流行病学研究中接受干预后多次随访观测结果流行病学研究中接受干预后
4、多次随访观测结果中国疾病预防控制中心中国疾病预防控制中心重复测量资料重复测量资料重复测量资料,是重复测量设计中所获取的多次测量资料。重复测量资料,是重复测量设计中所获取的多次测量资料。不同称谓:不同称谓:纵向数据纵向数据(Longitudinal data),如在生长发育研究中如在生长发育研究中重复测量资料重复测量资料(Repeated measures data),如在临床试验研究中如在临床试验研究中面板数据面板数据(Panel data),如在社会学研究中如在社会学研究中中国疾病预防控制中心中国疾病预防控制中心(1)平衡)平衡资资料:如果每一个体重复料:如果每一个体重复测测量的量的时间时间
5、和次数相同,和次数相同,称称为为平衡平衡资资料。料。(2)非平衡)非平衡资资料:不同个体重复料:不同个体重复测测量的量的时间时间点不同,或重点不同,或重复次数不同,称复次数不同,称为为非平衡非平衡资资料。料。重复测量资料重复测量资料中国疾病预防控制中心中国疾病预防控制中心 同一受试对象的同一观察指标在不同时间点上进行多次测量所得的资料,常用来分析该观察指标在不同时间点上的变化。有时是从同一个体的不同部位(或组织)上重复测量获得的指标的观测值。目的:就是比较不同时间点动态变化趋势的特征目的:就是比较不同时间点动态变化趋势的特征重复重复测量测量资料资料中国疾病预防控制中心中国疾病预防控制中心(1)
6、同一同一观观察察对对象的重复象的重复测测量量值值之之间间是非独立的,后一次是非独立的,后一次测测量的数据可能受前一次量的数据可能受前一次测测量量结结果的影响。果的影响。这这不符合不符合传统传统的的统计统计分析方法中关于独立性的假分析方法中关于独立性的假设设。(2)观观察指察指标标在所在所测测量的量的时间时间范范围围内可能成内可能成趋势趋势性性变变化。化。(3)观观察察值值的的变变异来源异来源较较多,有来自个体内的多,有来自个体内的变变异,个体异,个体间间的的变变异,或更高水平上的异,或更高水平上的变变异;异;变变异可能与异可能与时间时间有关,有关,或其他或其他协变协变量有关。量有关。重复测量资
7、料的特点重复测量资料的特点中国疾病预防控制中心中国疾病预防控制中心重复测量资料的数据特征重复测量资料的数据特征 当对同一受试对象在不同时间重复测量次数p3时,称为重复测量设计或重复测量数据。测测 量量 时时 间间 点点受试者受试者 1 2 p1 y11 y12 y1p2 y21 y22 y2p:n yn1 yn2 y n p中国疾病预防控制中心中国疾病预防控制中心中国疾病预防控制中心中国疾病预防控制中心实例举例1每一根线代表每一根线代表1只兔子只兔子中国疾病预防控制中心中国疾病预防控制中心实例举例2每一根线代表每一根线代表1位病人位病人中国疾病预防控制中心中国疾病预防控制中心典型的重复测量资料
8、典型的重复测量资料研究对象研究对象时间点时间点1时间点时间点2时间点时间点m1y11y12y1m2y21y22y2m3y31y32y3mnyn1yn2ynm中国疾病预防控制中心中国疾病预防控制中心典型的重复测量资料典型的重复测量资料分组因素分组因素研究对象研究对象时间点时间点1时间点时间点2时间点时间点m11y111y112y11m12y121y122y12m1n121y211y212y21m22y221y222y22m2n2p1yp11yp12yp1mp2yp21yp22yp2mpn3中国疾病预防控制中心中国疾病预防控制中心实际中:重复测量资料比独立资料更多见实际中:重复测量资料比独立资料更
9、多见临床研究中,需要观察病人在不同时间的某些生理、生化或病理指标的变化趋势,研究不同时间或疗程的治疗效果。流行病学研究中,观察队列人群在不同时间上的发病情况。研究不同职业、性别人群实施某种控制后,不同时间的多次效果考察。卫生学研究中,纵向观察儿童生长发育规律等,不同地区和环境营养状况。中国疾病预防控制中心中国疾病预防控制中心p简单处理简单处理:对平衡的重复测量资料,分别在各时间点上进对平衡的重复测量资料,分别在各时间点上进行分析行分析 p集合分析集合分析:将各个体的几次不同观察值相加,得到该个体将各个体的几次不同观察值相加,得到该个体的一个综合值,再进行比较分析的一个综合值,再进行比较分析 p
10、两阶段分析两阶段分析:为每个研究对象独立拟合时间函数,再对拟为每个研究对象独立拟合时间函数,再对拟合结果求均数。合结果求均数。重复测量资料的常见误用重复测量资料的常见误用中国疾病预防控制中心中国疾病预防控制中心减少样本含量减少样本含量控制个体变异控制个体变异非实验因素(干扰因素)非实验因素(干扰因素)滞留效应(滞留效应(carry-over effectcarry-over effect)潜隐效应(潜隐效应(latent effectlatent effect)学习效应(学习效应(learning effectlearning effect)重复测量设计优缺点重复测量设计优缺点中国疾病预防控制
11、中心中国疾病预防控制中心21重复测量设计与配对设计的区别?重复测量设计与配对设计的区别?中国疾病预防控制中心中国疾病预防控制中心1.随机区组设计要求每个区组内实验单位彼此独立区组号区组号 A营养素营养素 B营养素营养素 C营养素营养素 1 50.10 58.20 64.50 2 47.80 48.50 62.40 3 53.10 53.80 58.60 4 63.50 64.20 72.50 5 71.20 68.40 79.30 6 41.40 45.70 38.40 7 61.90 53.00 51.20 8 42.20 39.80 46.20 表表 A、B、C 3种营养素喂养种营养素喂养
12、 小白鼠所增体重(克)小白鼠所增体重(克)处理因素只能在区组内随机分配每个实验单位接受处理是不同的见右表:重复测量设计与随机区组设计的区别?重复测量设计与随机区组设计的区别?中国疾病预防控制中心中国疾病预防控制中心2.重复测量设计区组内实验单位彼此不独立中国疾病预防控制中心中国疾病预防控制中心同一受试的血样重复测量结果是高度相关提示:分析存在一定的复杂性。提示:分析存在一定的复杂性。中国疾病预防控制中心中国疾病预防控制中心小结 重复测量资料p数据:对同一受试对象的某个观察指标进行连续观测所得数据:对同一受试对象的某个观察指标进行连续观测所得p至少有两个因素:处理因素、时间因素;至少有两个因素:
13、处理因素、时间因素;p重复测量的试验结果重复测量的试验结果按时间顺序按时间顺序固定排列,不能像随机区固定排列,不能像随机区组设计的处理那样随机排列。组设计的处理那样随机排列。观观察指察指标标在所在所测测量的量的时间时间范范围围内可能成内可能成趋势趋势性性变变化;化;p同一个体不同时间测量值之间同一个体不同时间测量值之间高度相关;高度相关;p观观察察值值的的变变异来源异来源较较多,有来自个体内的多,有来自个体内的变变异,个体异,个体间间的的变变异,或更高水平上的异,或更高水平上的变变异;异;变变异可能与异可能与时间时间有关,或其他有关,或其他协变协变量有关。量有关。中国疾病预防控制中心中国疾病预
14、防控制中心第二节 重复测量资料方差分析中国疾病预防控制中心中国疾病预防控制中心想一想?想一想?同一观察单位具有多个观察值,而这些观察值来自同一受试对象的不同时点(部位等),这类数据间往往有相关性存在,违背了方差分析要求数据满足独立性基本条件。中国疾病预防控制中心中国疾病预防控制中心将总变异总变异分解为:个体间(个体间(between subjects)变异)变异 个体内个体内(within subject)变异变异,其中个体内变异是与重复因素有关的变量。重复测量资料的方差分析总思想重复测量资料的方差分析总思想中国疾病预防控制中心中国疾病预防控制中心平均值之间的多重比较 先采用配对t检验方法,计
15、算需比较的两两均数的t统计量,然后将这些样本统计量t值与Bonferroni临界t值进行比较。确定P值是否大于中国疾病预防控制中心中国疾病预防控制中心重复测量资料重复测量资料多变量重复测量方差分析单变量重复测量方差分析重复测量资料分类repeated measurement data中国疾病预防控制中心中国疾病预防控制中心多变量重复测量方差分析单变量重复测量方差分析重复测量资料分类重复测量资料分类repeated measurement data两组单组重复测量资料重复测量资料中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析模型重复测量方差分析模型 式中,为总均值,g表示第g个处理水平
16、的效应,j表示第j个测量时间点的效应,表示第g组第i个受试者的效应,该效应为随机效应,表示处理因素第g个水平在第j个测量时间点上的效应,也就是处理因素与时间点的交互效应,为随机误差。中国疾病预防控制中心中国疾病预防控制中心单变量重复测量方差分析单变量重复测量方差分析1.单组重复测量 指同一组内(或接受同一种处理)的多个受试者,在多个时间点上的反应变量所作的测量,又称为单变量重复测量。测测 量量 时时 间间 点点受试者受试者 1 2 p1 y11 y12 y1p2 y21 y22 y2p:n yn1 yn2 y n p中国疾病预防控制中心中国疾病预防控制中心2.多组重复测量(多组并不等于多因素)
17、指将受试者按处理的不同水平分为几个组,对这些组内的每一受试者,都在不同时间点对他们的反应变量进行测量。中国疾病预防控制中心中国疾病预防控制中心p重复测量资料方差分析的条件:1 1.正态性正态性 处理因素的各处理水平的样本个体之间是相互处理因素的各处理水平的样本个体之间是相互独立独立的的随机随机样本,其样本,其总体均数服从总体均数服从正态正态分布;分布;2.2.方差齐性方差齐性 相互比较的各处理水平的总体方差相等,即具有方差齐同相互比较的各处理水平的总体方差相等,即具有方差齐同 3.3.各时间点组成的各时间点组成的协方差阵协方差阵(covariance matrix)(covariance ma
18、trix)具有球形性(sphericity)特征。Box(1954)指出,若球形性质得不到满足,则方差分析的)指出,若球形性质得不到满足,则方差分析的F值是有偏的,这值是有偏的,这会造成过多的拒绝本来是真的无效假设会造成过多的拒绝本来是真的无效假设(增加增加型错误型错误)。重复测量资料方差分析的前提条件重复测量资料方差分析的前提条件中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析的球性假定重复测量方差分析的球性假定p方差是指在某一时点上测定值方差是指在某一时点上测定值变变异性的大小,而异性的大小,而协协方差是指在两个不同方差是指在两个不同时时点上点上测测定定值值相互相互变变异性异性的大
19、小。的大小。p如果在某个如果在某个时时点上的取点上的取值值不影响其他不影响其他时时点上的取点上的取值值,则协则协方差方差为为0,反之,反之,则则不不为为0。p由方差由方差协协方差构成的矩方差构成的矩阵阵称称协协方差方差阵阵。中国疾病预防控制中心中国疾病预防控制中心p设设k、l为为两个两个测测定定时时点,点,代表代表协协方差方差阵阵中的元素。当中的元素。当kl时为时为方差,方差,kl时为协时为协方差。共有方差。共有a个个测测定定时时点,将点,将这这a个方差和(个方差和(a-1)/2个个协协方差方差排成排成协协方差方差阵阵V为为:协方差阵的球形假定是指该矩阵主对角线元素(方差)协方差阵的球形假定是
20、指该矩阵主对角线元素(方差)相等、非主对角线元素(协方差)为零。相等、非主对角线元素(协方差)为零。重复测量方差分析的球性假定重复测量方差分析的球性假定中国疾病预防控制中心中国疾病预防控制中心所有两两时间点变量间差值对应的方差相等所有两两时间点变量间差值对应的方差相等对于yi与yj两时间点变量间差值对应的方差可采用协方差矩阵计算为:重复测量方差分析的球性假定重复测量方差分析的球性假定中国疾病预防控制中心中国疾病预防控制中心协方差阵协方差阵 A1A2A3A4A11051015A25201520A310153025A415202540s1-22=10+20-2(5)=20s1-32=10+30-2
21、(10)=20s1-42=10+40-2(15)=20s2-32=20+30-2(15)=20s2-42=20+40-2(20)=20s3-42=30+40-2(25)=20所有两两差值对应的方差相等,说明满足球形假定。重复测量方差分析的球性假定重复测量方差分析的球性假定中国疾病预防控制中心中国疾病预防控制中心p球性假定常用球性假定常用Mauchly法法进进行行检验检验pMauchly法法检验检验的的P值值若大于研究者所若大于研究者所选择选择的的显显著性水准著性水准时时,说说明明协协方差方差阵阵的球形性的球形性质质得到得到满满足。否足。否则则,必,必须对须对与与时间时间有关的有关的F统计统计量
22、的分子、分母自由度量的分子、分母自由度进进行行调调整,以整,以便减少犯便减少犯I类错误类错误的概率。的概率。p注意:注意:该该法法在小样本时很难检出差异在小样本时很难检出差异。重复测量方差分析的球性假定重复测量方差分析的球性假定中国疾病预防控制中心中国疾病预防控制中心对角线元素(方差)相等非主对角线元素(协方差)为零00中国疾病预防控制中心中国疾病预防控制中心球形对称的检验球形对称的检验用Mauchly法检验协方差阵是否为球形 H0:资料符合球形要求,H1:资料不满足球形要求检验的P值若大于研究者所选择的显著性水准时,说明协方差阵的球形性质得到满足。中国疾病预防控制中心中国疾病预防控制中心球形
23、条件不满足怎么办球形条件不满足怎么办常有两种方法可供选择:1.对重复测量ANOVA检验结果中与时间有关的F值的自由度进行调整(调小)2.采用MANOVA(多变量方差分析方法)(1)Geenhouse-Geisser调整系数(G-G)(2)Huynh-Feldt调整系数(H-F)分子自由度分子自由度分母自由度分母自由度中国疾病预防控制中心中国疾病预防控制中心自由度调整方法1中国疾病预防控制中心中国疾病预防控制中心自由度调整方法2中国疾病预防控制中心中国疾病预防控制中心调整规则 调整后的F临界值较原先大,提高了拒绝H0的门槛。减少了犯I类错误的概率。中国疾病预防控制中心中国疾病预防控制中心多变量方
24、差分析统计量目前多数统计软件会给出四个统计量目前多数统计软件会给出四个统计量:pWilks lambdapPillais trace、pHotelling-Lawley tracepRoys largest root。中国疾病预防控制中心中国疾病预防控制中心多变量方差分析统计量pWilkslambda指指标标的含的含义义是最明确的,其是最明确的,其计计算公式算公式为为|H|/(|H|+|E|)。反映。反映了了组组内方差占内方差占总总方差(方差(组间组间方差方差+组组内内方差)的比例方差)的比例。该该值值越小,越小,组间组间方差越大,意味着方差越大,意味着组间组间差异越差异越大大。p其它三个指其
25、它三个指标标的含的含义义并不是很并不是很直直观观:Hotelling-Lawleytrace是是HE-1矩矩阵阵的的对对角角线线方差之和。方差之和。Pillaistrace是是H(H+E)-1矩矩阵阵的的对对角角线线方差之和方差之和。Royslargestroot是是HE-1矩矩阵阵的最大的最大特征特征值值。中国疾病预防控制中心中国疾病预防控制中心第三节 重复测量方差分析案例及SAS实现中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析的SAS语句proc glm;class 分组变量;model 分析变量=分组变量;lsmeans 分组变量;means 分组变量;repeated 重复
26、测量变量名 /;中国疾病预防控制中心中国疾病预防控制中心 repeated 重复测量变量名重复测量变量名 /;【repeated】后必】后必须须指定重复指定重复测测量量变变量名,量名,测测量水平可量水平可选选,转换项转换项也是可也是可选选的。如:的。如:repeatedtime4;表示指定重复表示指定重复测测量量变变量名量名为为time,测测量量水平(即重复次数)水平(即重复次数)为为4次,中次,中间间均以空格隔开。其中均以空格隔开。其中time由由读读者自行命名,也可者自行命名,也可为为其它名称,如其它名称,如repeatedt4,repeatedmeasure4等均可。等均可。重复测量方差
27、分析的SAS语句中国疾病预防控制中心中国疾病预防控制中心repeated语语句后的句后的主要用于各次主要用于各次测测量之量之间间的比的比较较,常用的有两,常用的有两种种转换转换方式:方式:contrast和和polynomial。contrast表示指定一个参照水平,其它均与参照水平比表示指定一个参照水平,其它均与参照水平比较较。如如repeatedtime4contrast(1);表示共有表示共有4次重复次重复测测量,以第量,以第1次次测测量量为为参参照,第照,第2、3、4次次测测量分量分别别与第与第1次次测测量水平量水平进进行比行比较较。polynomial表示正交多表示正交多项项式式对对
28、比,即正交比,即正交转换转换后判断曲后判断曲线线的的n阶阶水平水平是否有是否有统计统计学意学意义义。如:。如:repeatedtime4polynomial;表示表示4次重复次重复测测量,量,分分别对别对一次、二次和三次曲一次、二次和三次曲线线分析是否有分析是否有统计统计学意学意义义,也就是判断,也就是判断测测量量值值随随时间变时间变化大致呈何种化大致呈何种趋势趋势。重复测量方差分析的SAS语句中国疾病预防控制中心中国疾病预防控制中心printe该该选选项项可可输输出出一一些些相相关关矩矩阵阵(不不是是我我们们所所关关心心的),更重要的是给出球性检验的),更重要的是给出球性检验结果。结果。su
29、mmary 输输出出重重复复测测量量各各水水平平间间的的比比较较结结果果。如如果果没没有有该该 选选 项项,即即 使使 前前 面面 指指 定定 了了 contrast或或polynomial,也不会显示比较结果,也不会显示比较结果重复测量方差分析的SAS语句中国疾病预防控制中心中国疾病预防控制中心例某研究例某研究评评价一医价一医疗仪疗仪器在治器在治疗颈疗颈椎病中的椎病中的疗疗效,效,以以VASVAS评评分作分作为评为评价指价指标标,共治,共治疗疗2 2个个疗疗程。从程。从该该研研究中随机究中随机选择选择2020例研究例研究对对象,表象,表1 1列出了列出了2020例研究例研究对对象的治象的治疗
30、疗前、治前、治疗疗1 1个个疗疗程、治程、治疗疗2 2个个疗疗程的程的VASVAS评评分分结结果,果,试试分析分析VASVAS评评分的分的变变化化趋势趋势。重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心ID疗疗前前1个个疗疗程后程后2个个疗疗程后程后ID疗疗前前1个个疗疗程后程后2个个疗疗程后程后1786117742443127663554137654665146545441157666665165457666176538877186659444197661087320875重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心dataexample1;inputvas
31、1-vas3;cards;procglm;modelvas1-vas3=/nouni;repeatedtime3contrast(1)/printesummary;run;重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心p从本例分析从本例分析结结果可以看出,三次果可以看出,三
32、次VAS评评分之分之间间存在一定的存在一定的相关性,不相关性,不应应采用采用单变单变量方差分析。重复量方差分析。重复测测量方差分析量方差分析结结果提示,三次果提示,三次评评分之分之间间存在一定的存在一定的趋势趋势。p如果我如果我们们想想进进一步了解存在什么一步了解存在什么样样的的趋势趋势,可以通,可以通过过在在repeated语语句中指定句中指定polynomial来来实现实现。procglm;modelvas1-vas3=/nouni;repeatedtime3polynomial/printesummary;run;重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心重复测量方差
33、分析案例1中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心图图1观测时间观测时间VAS评评分分疗前第1疗程后第2疗程后01234567重复测量方差分析案例1中国疾病预防控制中心中国疾病预防控制中心例例为为探索磁共振探索磁共振扩扩散加散加权权成像在胃成像在胃肠肠癌肝癌肝转转移化移化疗疗疗疗效中的价效中的价值值,某研究者,某研究者观观察了察了20例化例化疗疗反反应应良良好和化好和化疗疗反反应应不良的胃不良的胃肠肠癌患者分癌患者分别别在化在化疗疗前、化前、化疗疗后后3天、天、7天、天、42天的成像,并天的成像,并计计算表算表观扩观扩散系数散系数A
34、DC值值。表中。表中group=0代表反代表反应应良好良好组组,group=1代代表反表反应应不良不良组组。试试比比较较两两组组人群的人群的ADC值值差异是否差异是否有有统计统计学意学意义义。重复测量方差分析案例2中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例2中国疾病预防控制中心中国疾病预防控制中心dataexample2;inputgadc1-adc4;cards;01.86 1.821.831.8601.94 1.771.881.7311.13 1.211.241.2510.91 1.051.191.16;procglm;classg;modeladc1-adc4=g/no
35、uni;repeatedtime4(03742)polynomial/printesummary;run;注意,当注意,当repeated语句指定语句指定polynomial变换时,如果变换时,如果观测时间间隔不等,必须在观测时间间隔不等,必须在括号中加入具体的时间点。括号中加入具体的时间点。因为正交多项式系数在等间因为正交多项式系数在等间隔和不等间隔时是截然不同隔和不等间隔时是截然不同的。如果不指定,则默认为的。如果不指定,则默认为等间隔,计算的正交多项式等间隔,计算的正交多项式结果会有偏差。结果会有偏差。重复测量方差分析案例2中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例2中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例2中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例2中国疾病预防控制中心中国疾病预防控制中心重复测量方差分析案例2中国疾病预防控制中心中国疾病预防控制中心图图2观测时间观测时间点点ADC值值d0d3d7d420.00.20.40.60.81.01.21.41.6反应良好反应不良重复测量方差分析案例2中国疾病预防控制中心中国疾病预防控制中心
限制150内