医学统计学与SAS软件方差分析.ppt

医学统计学与SAS软件 佳木斯大学公共卫生学院 流行病与卫生统计教研室 方差分析 在SAS软件中，通常进行方差分析采用 ANOVA和GLM过程，前者运算速度较快， 但功能较为有限；后者运算速度较慢，但 功能强大。二者语法应用上没有大的差异 。 data tmp; do a=1 to 3; do i=1 to 12; input x; output; end; end; datalines; 332.96 297.64 312.57 295.47 274.25 307.97 292.12 244.61 261.46 286.46 322.49 282.42 263.21 235.87 269.30 258.90 254.39 200.87 227.79 237.05 216.85 238.03 238.19 243.49 232.55 217.71 216.15 220.72 219.46 247.47 280.75 196.01 208.24 198.41 240.35 219.56 ; run; proc anova data=tmp; class a; model x=a ; means a /hovtest snk; quit; proc glm data=tmp; class a; model x=a ; means a /hovtest snk lsd dunnett(3); quit; LSD两两比较结果 Dunnett比较结果，与第3组比较 随机区组设计的方差分析 data ex4_4; input x a b ; cards; 0.82 1 1 0.65 2 1 0.51 3 1 0.73 1 2 0.54 2 2 0.23 3 2 0.43 1 3 0.34 2 3 0.28 3 3 0.41 1 4 0.21 2 4 0.31 3 4 0.68 1 5 0.43 2 5 0.24 3 5 ; proc anova; class a b; model x=a b; means a/snk; quit; 析因设计的方差分析 data xiyin; do a=1 to 2; do b=1 to 2; do n=1 to 4; input x; output; end; end; end; datalines; 0.163 0.199 0.184 0.198 0.127 0.168 0.152 0.150 0.124 0.151 0.127 0.101 0.101 0.192 0.079 0.086 ; run; proc glm data=xiyin; class a b; model x = a b a*b; means a /snk; quit; 练习题： 1、某医生为了解一种降血脂新药的临床疗效，按统一纳入标准选择120名高血 脂患者，采用完全随机设计方法分为4组，6周治疗后结果如下，问4组低密度 脂蛋白含量的总体均数有无差别？ 分组测量值 安慰剂4.38 6.26 5.84 3.77 3.98 5.37 3.9 7.19 2.4 4.04 3.63 3.68 4.23 5.86 7.28 5.83 4.27 6.21 5.88 3.75 5.28 4.11 3.22 3.46 6.57 3.48 2.67 5.57 4.86 3.49 2.4g药物 6.02 7.46 6.11 7.38 9.36 7.88 5.64 5.54 4.19 8.59 7.17 7.96 5.97 6.55 6.39 5.1 8.71 4.45 5.59 6.72 6.99 4.75 8.85 7.5 7.03 8.67 7.09 8.91 6.51 2.7 4.8g药物 9.33 6.62 7.73 4.6 7.85 7.23 6.75 4.5 8.99 5.67 7.02 7.19 6.54 9.52 8.41 9.14 7.04 3.51 10.07 6.95 5.61 7.63 9.2 4.87 5.71 6.53 7.74 8.38 8.26 6.12 7.2g药物 5.4 7.78 9.38 9.57 7.68 10.07 10.98 9.43 8.91 9.1 8.81 13.19 10.42 8.97 8.37 8.09 10.07 8.08 11.81 9.5 11.93 8.48 10.57 9.38 11.57 12.9 7.97 8.18 10.87 6.42 2、某研究者采用随机区组设计进行试验，比较三种抗癌药 物对小白鼠肉瘤的抑瘤效果，先将15只染有肉瘤的小白鼠按 体重大小配成5个区组，每个区组内3只小白鼠随机接受三种 抗癌药物，以肉瘤的重量为指标，试验结果见表。问三种不 同药物的抑瘤效果有无差别？ 区组A药B药C药 10.820.651.98 20.730.541.50 30.430.341.05 40.410.210.93 50.680.431.35 将20只家兔随机分成4组，每组5只，进行神经损伤后的 缝合实验，处理由两个因素组合而成，A因素为缝合方法 ，有两个水平， 一个水平为外膜缝合，记作a1，另一个 水平为束膜缝合，记作a2；B因素为缝合后的时间，有两 个水平，一水平为缝合后1个月，记作b1，另一水平为缝 合后2个月，