第十六章析因设计和交叉设计资料分析Stata实现

第十六章析因设计和交叉设计资料分析的Stata实现本章应用的 Stata 命令为：方差分析anova varname term / term /.重组交叉设计数据pkshape id sequence periodl period2分析交叉设计数据pkcross outcome例16-3 治疗白细胞降低病人16例，随机分为4组，给予相应治疗措施：基础 药物（对照组）、基础药物+甲药、基础药物+乙药、基础药物+甲药+乙药，3个 月后观察白细胞增加数（109/L）,以评价甲药和乙药以及它们联合用药的作用， 得到结果如表16-3。表16-3治疗白细胞降低四种不同疗法3个月后白细胞增加数第一组（g=l）a = 0, b = 0（基础治疗）第二组（g=2）a = 1, b = 0 （基础治疗+甲药）第三组（g3）a = 0, b = 1 （基础治疗+乙药）第四组（g=4）a = 1, b = 1 （基础治疗+甲药+乙药）1.001.481.422.32Xij0.601.421.282.380.881.110.902.120.701.300.982.21Stata 数据如下：abgx00110010.60010.880010.71021.481021.421021.111021.30131.420131.280130.90130.981142.321142.381142.121142.21Stata 命令为： anova x a b a*b 结果为：Number of obs=16R-squared=0.9218Root MSE=.181968Adj R-squared=0.9022Source |1Partial SSdfMSFProb > F1Model |4.6808249831.5602749947.120.0000a |2.7060249912.7060249981.720.0000b |1.6383999611.6383999649.480.0000a*b |.3364000321.33640003210.160.0078Residual |1.39735000812.0331125011Total |5.0781749915.338544999计算残差predict e,residual对残差进行正态性检验（«二0.05 ）swilk eShapiro-Wilk W test for normal dataVariable | ObsWVz Probz+e |160.948201.0490.096 0.46179P=0.46179>a，按正态分布假定进行统计分析。由于交付作用项a*b的P=0.0078，所以可以认为存在交互作用。对例16-3作简单效应检验（Bonferroni校正）和方差齐性检验（a = 0.10）。oneway x g , b t结果为：|Summary of xg |Mean Std. Dev.Freq.+1 |.795.1791647242 |1.3275.1631716843 |1.145.2456962744 |2.2575.115578324+Total |1.38125 581846216Analysis of VarianceSourceSSdf MSF Prob > FBet ween groups4.6808249831.5602749947.120.0000Within groups.39735000812.033112501Total5.0781749915.338544999Bartlett's test for equal variances: chi2(3) =1.4655 Probchi2 = 0.690Comparison of x by g(Bonferroni)Row Mean-|Col Mean |123+2 |II.53250.0083 I.35-.1825II0.1121.0004 I1.4625.931.1125I0.0000.0000.000方差齐性检验的P=0.690>0.10，可以按方差齐性的假定进行统计分析。两两比较组别1组别2组别1组别2Bonferroni 校正P值不用A药且不用B药(g=1)用A药但不用B药(g=2)-0.53250.008不用A药但用B药(g=3)用A药且用B药(g=4) 1.1125<0.001不用A药且不用B药(g=1)用B药但不用A药(g=3)0.350.112不用B药但用A药(g=2)用A药且用B药(g=4) 0.93<0.001因此可以认为A药和B药均有效，并且有协同交互作用。例16-4 欲研究酰托普利片的降压作用，对照药为安慰剂。选择28例高血压 患者，随机分成两组，每组14例，一组先用安慰剂，经过一段清洗期后再用酰托 普利片（记为AB组），另一组先用酰托普利片，经过一段清洗期后再用安慰剂（记 为BA组）。效应指标为用药前后收缩压的降低值，得到结果如表16-10。表16-10 酰托普利片的降压作用研究（2X2交叉设计）受试者编号试验次序时期1时期2受试者编号试验次序时期1时期21AB142815BA28122AB192216BA22183AB102717BA27184AB202918BA29165AB61519BA23196AB52020BA22117AB142821BA28128AB152522BA1999AB82323BA231310AB172924BA291911AB52725BA27012AB211026BA281213AB251927BA241714AB122228BA170相应的假设检验如下：“ 0.051、顺序间效应差异H : AB和BA两种顺序药物的降压的总体均数相等。 01H : AB和BA两种顺序的降压的总体均数不相等。 112、药物间效应差异H : A和B两种药物的降压的总体均数相等。 02H : A和B两种药物的降压的总体均数不相等。123、时期间效应差异H ：不同时期药物的降压的总体均数相等。 03H ：不同时期药物的降压的总体均数不相等。134、个体间的效应差异H ：不同个体的降压的总体均数相等。04H ：不同个体的降压的总体均数不相等14Stata 数据：idseqp1p21AB14282AB19223AB10274AB20295AB6156AB5207AB14288AB15259AB82310AB172911AB52712AB211013AB251914AB122215BA281216BA221817BA271818BA291619BA231920BA221121BA281222BA19923BA231324BA291925BA27026BA281227BA241728BA170Stata 命令和结果为：pkshape id seq pl p2.tab treat ,su(outcome)|Summary of outcometreat |Mean Std. Dev.Freq.+1 |13.1071436.1904304282 |123.9285714.822236928+Total |18.5178577.748292656.pkcross outcome sequence variable = sequence period variable = period treatment variable = treat carryover variable = carry id variable = idAnalysis of variance (ANOVA) for a 2x2 crossover studySource of Variation| Partial SS丄dfMSFProb > FIntersubjects|Sequence effect|0.8710.870.030.8741Residuals|1888.612634.181.190.3326Intrasubjects+|Treatment effect|1639.4511639.4556.940.0000Period effect|24.45124.450.850.3653Residuals|+748. 612628. 79Total |3301.9855Omnibus measure of separability of treatment and carryover =29.2893%方差分析结果表明：1. 不同顺序间效应的差异无统计学意义(F=0.03,P=0.8741)，表明无延滞效 应。2. 不同时期间效应的差异无统计学意义(F=0.85,P=0.3653)，表明处理的效 应不随时期而变化。3. 不同处理间的差异有统计学意义(F=56.94, PvO.OOOl)，说明A药物与B 药物降压是有差别的。4.个体(顺序)的差异无统计学意义(F=1.19,P=0.3326),说明个体之间没有差异。