《第十六章析因设计和交叉设计资料分析Stata实现》由会员分享，可在线阅读，更多相关《第十六章析因设计和交叉设计资料分析Stata实现（6页珍藏版）》请在装配图网上搜索。
1、第十六章析因设计和交叉设计资料分析的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
2、）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.21St
3、ata 命令为： 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.03
4、31125011Total |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.46179a，按正态分布假定进行统计分析。由于交付作用项a*b的P=0.0078，所以可以认为存在交互作用。对例16-3作简单效应检验（Bonferroni校正）和方差齐性检验（a = 0.10）。oneway x g , b t结果为：|Sum
5、mary 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.338544999Bartletts test for e
6、qual 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.6900.10，可以按方差齐性的假定进行统计分析。两两比较组别1组别2组别1组别2Bonferroni 校正P值不用A药且不用B药(g=1)用A药但不用B药(g=2)-0.53250.008不用A药但用B药(g=3)用A药且
7、用B药(g=4) 1.11250.001不用A药且不用B药(g=1)用B药但不用A药(g=3)0.350.112不用B药但用A药(g=2)用A药且用B药(g=4) 0.93 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),说明个体之间没有差异。
第十六章析因设计和交叉设计资料分析Stata实现
