医学统计学与软件方差分析.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