实验四方差教案分析.pptVIP

实验四 方差分析实验 流行病与卫生统计学教研室 【例】 为研究煤矿粉尘作业环境对尘肺的影响，将24只Wistar 大鼠随机分到甲、乙、丙三个组，每组8只，分别在地面办公楼、煤炭仓库和矿井下染尘，12周后测量大鼠全肺湿重（g）。 案例分析 甲组 4.2 3.3 3.7 4.3 4.1 3.3 3.5 4.1 乙组 4.5 4.4 3.5 4.2 4.6 4.2 4.4 4.1 丙组 5.6 3.6 4.5 5.1 4.9 4.7 4.8 4.4 表 三组大鼠的全肺湿重（g） 请思考以下问题 该实验属何种设计方案？处理因素是什么？有几个水平？观察指标是什么? 能否采用 t 检验比较不同作业环境中的大鼠全肺湿重是否有差异？ 对于小样本多组均数的比较不能采用t检验进行两两比较 原因： 割裂整体设计，只见树木，不见森林 增大一型错误的概率 可以采用方差分析的方法进行分析 一、方差分析的应用条件 方差分析的应用条件 独立性：各样本是相互独立随机的样本 正态性：各样本都来自正态总体 方差齐性：各样本的总体方差相等 二、方差分析思想 方差分析思想 根据资料设计的类型及研究目的，可将总变异分解为两个或多个部分，每个部分的变异可由某因素的作用来解释。通过比较可能由某因素所至的变异与随机误差，即可了解该因素对测定结果有无影响。 变异分解(以完全随机设计为例) 组间变异 总变异 组内变异 SS总=SS组间+SS组内 ?总=?组间＋?组内 总变异——SS总（离均差平方和 ） ?总=N-1 组间变异——SS组间 ?组间=k-1 MS组间 =SS组间/(k-1) 组内变异——SS组内 ?组内=N-k MS组内=SS组内/(N-k) 如果三组样本均数来自同一总体(H0:?1=?2=?3)，那么从理论上说组间变异应该等于组内变异，因为两者均只反映随机误差(包括个体差异)，这时若计算组间均方与组内均方的比值： F= MS组间 /MS组内 则F值在理论上应等于1，但由于抽样误差的影响，F通常接近1，而并不正好等于1。相反，若三种疗法效果不同，则组间变异就会增大，F值则明显大于1，要大到什么程度才有统计学意义呢？可通过查附表3 方差分析用F界值表得到P值，将其与事先规定的?值比较后作出判断。 三、不同类型的方差分析 完全随机设计资料的方差分析 随机区组设计资料的方差分析 1、完全随机设计资料的方差分析 2、随机区组设计资料的方差分析 四、多个样本均数的两两比较 当方差分析的结果为拒绝H0，接受H1时，只说明g个总体均数不全相等。若想进一步了解哪些两个总体均数不等，需进行多个样本均数间的两两比较或称多重比较（multiple comparison）。 为什么不一开始就进行均数的两两比较？ SNK-q检验：样本均数两两之间的全面比较，探索性研究。 Dunnett-t检验：它适用于k–1个实验组与对照组均数的比较，验证性研究。 Bonfferoni t检验 ：两两比较方法中最为保守的，不适合分组很多、比较次数太多的资料。 五、方差齐性检验 Bartlett ?2检验 Levene检验 定量资料假设检验方法的选择 P0.05 P0.05 P0.05 P0.05 No Yes 多个样本均数比较 两样本均数比较 方差齐否？ 正态否？ F检验、t检验 秩和检验 变量变换 不拒绝Ho 两两比较 案例1 课本案例7-1、7-2 作业： P148页： 分析计算题： 第1题 第4题 * * 表8-1 三种不同喂养方式下大白鼠体重喂养前后差值（g） 正常钙（0.5%） 中剂量钙（1.0%） 高剂量钙（1.5%） 合计 332.96 253.21 232.55 297.64 235.87 217.71 312.57 269.30 216.15 295.47 258.90 220.72 284.25 254.39 219.46 307.97 200.87 247.47 292.12 227.79 280.75 244.61 237.05 196.01 261.46 216.85 208.24 286.46 238.03 198.41 322.49 238.19 240.35 282.42 243.49 219.56 n 12 12 12 36（N） 293.37 239.49 224.78 252.55() 606.15 350.51 540.31 1364.52（） 表8-1 完全随机设计资料方差分析表 变异来源 SS ( MS F 总变异 N-1 组间变异 k-1 SS组间/(k-