《析因设计和交叉设计资料分析Stata实现》由会员分享，可在线阅读，更多相关《析因设计和交叉设计资料分析Stata实现（6页珍藏版）》请在金锄头文库上搜索。
1、第十六章 析因设计和交叉设计资料分析本章应用的STATA命令为：方差分析anova varname term / term / .重组交叉设计数据pkshape id sequence period1 period2分析交叉设计数据pkcross outcome例16-3 治疗白细胞降低病人16例，随机分为4组，给予相应治疗措施：基础药物（对照组）、基础药物+甲药、基础药物+乙药、基础药物+甲药+乙药，3个月后观察白细胞增加数（109/L），以评价甲药和乙药以及它们联合用药的作用，得到结果如表16-3。表16-3 治疗白细胞降低四种不同疗法3个月后白细胞增加数第一组(g=1)（基础治疗）第二组
2、(g=2)(基础治疗+甲药)第三组(g3)(基础治疗+乙药)第四组(g=4)(基础治疗+甲药+乙药)1.000.600.880.701.481.421.111.301.421.280.900.982.322.382.122.21解：STATA数据如下：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 = 16 R-squared = 0.921
3、8 Root MSE = .181968 Adj R-squared = 0.9022 Source | Partial SS df MS F Prob F -+- Model | 4.68082498 3 1.56027499 47.12 0.0000 | a | 2.70602499 1 2.70602499 81.72 0.0000 b | 1.63839996 1 1.63839996 49.48 0.0000 a*b | .336400032 1 .336400032 10.16 0.0078 | Residual | .397350008 12 .033112501 -+- Tot
4、al | 5.07817499 15 .338544999 计算残差predict e,residual对残差进行正态性检验()swilk e Shapiro-Wilk W test for normal data Variable | Obs W V z Probz-+- e | 16 0.94820 1.049 0.096 0.46179P=0.46179a，按正态分布假定进行统计分析。对例16-3作简单效应检验（Bonferroni校正）和方差齐性检验（）。oneway x g , b t结果为： | Summary of x g | Mean Std. Dev. Freq.-+- 1
5、| .795 .17916472 4 2 | 1.3275 .16317168 4 3 | 1.145 .24569627 4 4 | 2.2575 .11557832 4-+- Total | 1.38125 .5818462 16 Analysis of Variance Source SS df MS F Prob F-Between groups 4.68082498 3 1.56027499 47.12 0.0000 Within groups .397350008 12 .033112501- Total 5.07817499 15 .338544999Bartletts test
6、 for equal variances: chi2(3) = 1.4655 Probchi2 = 0.690 Comparison of x by g (Bonferroni)Row Mean-|Col Mean | 1 2 3-+- 2 | .5325 | 0.008 | 3 | .35 -.1825 | 0.112 1.000 | 4 | 1.4625 .93 1.1125 | 0.000 0.000 0.000方差齐性检验的P=0.6900.10，可以按方差齐性的假定进行统计分析。两两比较组别1组别2组别1组别2Bonferroni校正P值不用A药且不用B药(g=1)用A药但不用B药(
7、g=2)0.53250.008不用A药但用B药(g=3)用A药且用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.930.001因此可以认为A药和B药均有效，并且有协同交互作用。例16-4 欲研究酰托普利片的降压作用，对照药为安慰剂。选择28例高血压患者，随机分成两组，每组14例，一组先用安慰剂，经过一段清洗期后再用酰托普利片（记为AB组），另一组先用酰托普利片，经过一段清洗期后再用安慰剂（记为BA组）。效应指标为用药前后收缩压的降低值，得到结果如表16-10。表16-10 酰托普
8、利片的降压作用研究（22交叉设计）受试者编号试验次序时期1时期2受试者编号试验次序时期1时期21AB142815BA28122AB192216BA22183AB102717BA27184AB202918BA29165AB61519BA23196AB52020BA22117AB142821BA28128AB152522BA1999AB82323BA231310AB172924BA291911AB52725BA27012AB211026BA281213AB251927BA241714AB122228BA170相应的假设检验如下：1、顺序间效应差异：AB和BA两种顺序药物的降压的总体均数相等。：AB
9、和BA两种顺序的降压的总体均数不相等。2、药物间效应差异：A和B两种药物的降压的总体均数相等。：A和B两种药物的降压的总体均数不相等。3、时期间效应差异：不同时期药物的降压的总体均数相等。：不同时期药物的降压的总体均数不相等。4、个体间的效应差异：不同个体的降压的总体均数相等。：不同个体的降压的总体均数不相等。解：STATA数据：idseqp1p21AB14282AB19223AB10274AB20295AB6156AB5207AB14288AB15259AB82310AB172911AB52712AB211013AB251914AB122215BA281216BA221817BA271818BA291619BA231
