小样本资料的差异显著性检验课件

书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 1 第六章 小样本资料的差异显著性检验 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 2 本章主要介绍小样本时单个均数 两个均数的假设检验 单个率 两个率间的假设检验 应重点掌握各种情况下的t检验方法 正确区分成组资料和配对资料 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 3 在上一章中 我们系统介绍了抽样分布和统计推断的基本原理和基本方法 即通过随机抽样的方法获得某一特定总体的随机样本 用这一样本进行试验 并对试验后所得资料进行分析 通过统计推断来定性或定量地分析研究总体的特征本章主要介绍总体方差未知 且样本较小时不同资料类型的统计推断 差异显著性检验 假设检验 的具体方法 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 4 第一节单个平均数的假设检验 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 5 单个平均数的假设检验是检验一个样本所属的总体平均数与一个特定 已知 总体平均数间是否存在显著差异的一种统计方法 也可理解为检验一个样本是否来自某一特定 已知 总体的统计分析方法根据统计假设检验的基本原理可知 假设检验的关键是根据统计量的分布计算实得差异 即表面效应 由抽样误差造成的概率测验的统计量分布服从u 分布或t 分布 所以单个平均数的假设检验可分为u test和t test两种 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 6 一 总体方差已知时单个平均数的假设检验 当总体方差已知 不管样本多大 均可用u 分布计算实得差异由抽样误差造成的概率 所以称u testu 检验 u test 的方法和步骤见前一章内容 请逐一回忆一下检验公式和检验步骤 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 7 二 总体方差未知时单个样本平均数的假设检验 一 当总体方差未知 而样本较大时 可以使用u 分布计算实得差异由抽样误差造成的概率因此 其检验还是u test 请回忆一下以上两种u test的共同点和不同点 公式的不同 哪里不同 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 8 二 总体方差未知 且样本较小时的单个样本平均数的假设检验实际上 在很多情况下 总体方差往往是未知的 而由于各种条件的限制 试验的样本又不可能很大 即只能用小样本来作试验 或调查时抽取的样本较小因此 总体方差未知 样本较小是试验中最常见的一种情况 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 9 在第四章讨论t 分布时 我们已经知道 总体方差未知 且样本较小时 可以用代替 其统计量就不再服从标准正态分布 而是服从t 分布 请回忆一下t 分布曲线及其特点 t 分布曲线受自由度制约 不同自由度下的t 分布曲线其形状是不同的 因此不同自由度下算得的t 值落在某一范围内的概率值也随自由度的不同而不同 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 10 下面我们用例子来说明这一类型的t test例 某一地区根据多年资料 得出4龄草鱼的平均体重为 4 3kg 今在该地捕得一批同龄草鱼 n 12 体重为 4 05 4 25 4 10 4 20 4 30 4 20 4 35 4 10 4 25 4 15 4 30 4 00kg问 该批草鱼在体重方面是否符合该地区4龄草鱼的特征 即 该批草鱼饲养水平是否正常 该例样本量不大 n 12 30 因此符合总体方差未知 且是小样本的情况 应使用t test来进行检验 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 11 已知 4 3kg 样本 4 19kg 0 109kg检验步骤如下 设kgkg计算和t值 查附表4 t分布表 得知 自由度为df 12 1 11时的 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 12 本例中所得所得t值的概率因此 应否定无效假设 接受备择假设即 该批草鱼的体重显著不符合该地区4龄草鱼的体重特征由于所得t值较大 还需要和t0 01相比较 得 因此 应在0 01水平上否定无效假设 即该批草鱼的体重极显著不符合该地区4龄草鱼的体重特征 故应寻找原因 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 13 再举一例 药典规定 每100g某药物中应含有60mg的总黄酮 现对某药厂生产的某一批次的这种药进行检测 得如下数据 试问该批次药物的总黄酮含量合格吗 58 359 258 660 660 359 559 158 061 058 9n 10下面我们作统计分析由于是小样本 且总体方差为未知因此应使用t test进行分析 请同学们先自行立题 计算 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 14 首先计算平均数和标准差 得 第一步 设立无效假设设 第二步 计算t值 查t值表 得 即第三步 结论 接受无效假设 即该批次药物的总黄酮含量符合药典规定 试想一下 该题可以用一尾检验吗 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 15 上例也可以设立一个正常范围 检查样本的平均值是否在这一范围内我们有三个公式用于单个样本的统计假设检验 请说出这三个公式各自的使用场合和它们之间的区别 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 16 第二节两个样本平均值相比较的统计假设检验 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 17 很多情况下 我们不只是将样本平均数与总体平均数相比较而是做一个试验 这个试验中设置两个组 一个组作试验 施加某种试验条件 即处理 另一个组作对照 试验完了将这两个组的试验数据进行比较这种试验可以采用两种方法进行 一种方法是这两个组是相互独立的 一组作处理 一组作对照另一种方法是配对试验 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 18 下面我们先讨论两个组是相互独立的情况这种统计假设检验的方法称为成组数据的比较一 成组数据的平均值比较方法介绍 在一个总体中 随机地抽取两个样本 在这两个样本中 被抽取的个体是相互独立的 且基本条件 如品种 日龄 性别 体况等 都应一致 均匀 必要时须作适当的调整 尽可能使两个组在样本量上一致 各组的基本情况一致随机地指定一组作处理 另一组作对照 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 19 一定要注意其中的任意一组作处理 一组作对照 而不能事先规定哪一组做处理 哪一组作对照试验过程中注意记录资料 结束以后整理资料并进行统计分析这样得到的资料称为成组数据这样的数据组间 组内都是独立的成组数据的t test其公式是 其中 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 20 或 也可以这样写 其中 当两样本量相等时 则 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 21 在小样本时 两样本平均数差的标准化是服从t分布的因此 由于我们在之前的无效假设是 备择假设是 因此这一式子可以简化为 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 22 得到成组数据所进行的试验称为完全随机设计法 仅两个组 两个组的样本量可以相等 也可以不等 但应尽量接近这种两个组是完全独立的试验 如 一组用药 一组不用药一组用试验药物 一组用常规药物一组用试验剂量 一组用常规剂量一组是引进品种 品系 一组是本地品种 等等就称为完全随机设计试验这里 前面的一组称为处理 后一组称为对照 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 23 下面我们以实例来说明成组数据的比较A B两厂生产某同类鱼药 现作24小时鱼体内残留量检测 得如下数据 试分析哪一厂家的该类药物的残留量大A厂 49 348 149 851 250 050 749 948 550 451 6B厂 48 349 748 248 847 347 750 449 2显然 这两个厂家的药物是相互独立的 试验所用动物也是独立的 样本较小因此应使用成组数据的t test进行分析 同学们先行分析 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 24 先做预备计算 将两个样本的平均值 方差等计算出来 设立无效假设 计算t值 并作比较 所得t值出现的概率因此 否定无效假设 接受备择假设 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 25 即 A B两厂生产的该类药物的24h血液残留量差异显著 下面应针对这一现象作出专业解释 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 26 在这种检验中 我们求t值时用的是合并均方 合并均方只有在两总体方差 我们一般用样本均方估计总体方差 相同 即两方差差异不显著的情况下才能得到这里我们总假定两个均方差异不显著 但如果两方差差异显著 就不能合并 两均方是否相等 必须用下一章的F test进行检验才能知道当两均方不等 称为方差不齐 时 用以下方法计算和t值 当用的作判断的临界值当时须用Cochran cox法 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 27 首先计算在 水平上显著的临界值式中 若就否定 否则 就接受由于处于之间 因此 只有在实际计算得到的在之间时 才需要计算 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 28 附 两均方是否齐性的判别方法 如果表示两均方齐性 否则就是不齐 例题见P52该例实际是不需要计算的 为什么 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 29 二 配对数据的平均值比较方法介绍 配对数据来自于配对试验 配对试验根据具体试验情况可以有好多种方法 1 将两个品种 性别 日龄 体况等一致的动物 最好是有血缘关系的同胞或半同胞 配成一对 任意一只进入试验组 另一只进入对照组 配成若干对 试验中做好记录 这种记录到的数据就是配对数据注意 对子内的两个个体应尽可能一样 但对子之间应有较大的差异 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 30 2 选取若干个动物 每一个动物在试验前测定一次 试验后测定一次 这样的两次记录就是配对数据 如同一个人吃早饭前后各测一次血糖值 这同一个人的两次血糖值就是一对数据 若干个人就有若干个数据对3 同一个体的不同部位可以配成一对 如一个兔子的左右体侧 就可以配成一对 若干个兔子的不同体表就是若干对4 同一个动物的不同试验时期所施加的不同处理形成一个对子 如一只猪在第一试验期施加A处理 在第二试验期施加B处理 或第一试验期施加处理 第二试验期作对照 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 31 凡此种种 都是配对试验 因此配对试验既有空间上的配对 亦有时间上的配对 配对试验所得到的数据 就称为配对数据配对数据的比较方法不同于成组数据的比较方法其t test的公式为 其中 而是对子数 这里的样本量不是2n 而是n因此自由度不是2 n 1 而是n 1 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 32 在对配对资料进行统计分析时 首先要计算每一对对子内两数据之差 即d值 然后对d值进行分析下面我们以实际例子来说明配对数据的差异性检验研究和的关系 挑选体况基本相似的一对全同胞小鼠 其中任意一只放入 组 另一只放入 组 共挑选了8对 分成了两组一组饲喂正常饲料 另一组饲喂缺乏的饲料 试验结束后检测小鼠肝脏中的含量 得如下数据 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 33 对子号12345678正常35502000300039503800375034503050缺24502400180032003250270025001750差1100 400120075055010509501300第一步 设立无效假设 即饲料中VD缺乏不影响肝脏中VA的储存量 对 饲料中缺乏VE会严重影响肝脏中VA的储存量简写为 设 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 34 第二步 计算 首先计算各组数据的差然后计算 得得 即t值出现的概率否定无效假设 接受备择假设 即饲料中缺乏VE会严重影响小鼠肝脏对VA的储存量 或正常饲料与缺乏VE的饲料两者对肝脏VA的储存量差异极显著 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 35 随机抽取10名健康青年男子 测量一日中不同时间的血压 舒张压 变化 得如下数据 试分析不同时间血压的差异性序号123456789106时808375787189778172759时85938983808685797881差 5101459 38 266上一例是两个试验动物配成一对 这一例是同一个试验材料不同的时间内测定的数据配成一对 这是在水产诊断试验中经常使用的对比方法下面是统计分析结果 同学们可先行分析 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 36 设立无效假设 计算t值 先计算每一个人的差值 写于表格的最下面一行再计算差的平均值和标准误 计算t值 即出现t值的概率因此 否定无效假设 接受备择假设 上午9时血压舒张压极显著高于上午6时 或 上午6时和9时的血压舒张压有极显著的差异 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 37 差异不显著 显著 极显著也可以这样表示 不显著 n s 显著 极显著 如 成组数据差异显著性检验和配对数据差异显著性检验的区别在于 1 在作试验设计和实施试验时 两者的差别就已经确定了 试验设计是配对的 就不能用成组比较法来分析所得资料 反之 试验设计是成组的 就不能用配对比较法来分析所得资料 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 38 2 如果是分析别人的资料 如我们在统计学习中碰到的例题等 一定要根据例题 或习题 的内容来加以判断 独立抽取两个样本 将所抽得的样本随机分为两组 抽取一对相似的试验动物 一个个体在某一时间段作什么 另一时间段作什么 某一试验前测定一次 试验后测定一次 等 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 39 两者的区别千万不能搞错 配对试验如使用成组比较法 易发生 型错误 即不能鉴别应属显著的差异 这是人为地扩大了标准误 存伪 成组试验的资料如果乱配成对子后用配对比较法 易使不显著的差异检验成显著的差异 弃真 成组试验应采用供试单位基本一致的材料 彼此独立 随机抽样 随机分组 随机供试 两组样本尽可能一致配对试验对子内应尽可能一致 对子间应适当扩大距离 以使试验有更大的适应性 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 40 第三节百分资料的u test 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 41 在动物科学研究中 特别是水产临床诊断研究中 有很多资料是属于二项分布的这样的资料一般可用率来表示 如受精率 出苗率 孵化率 死亡率 存活率 治愈率 淘汰率 等计算这一类率时所用的样本一般很大 样本量较小时 所得到的率其实用意义不大当p或1 p不太小 而np或n 1 p 不小于5时 率的分布也服从正态分布 因此 率的检验可采用u test率的检验有单个样本率和两个率的比较两种 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 42 一 单个样本率的假设检验这是检验某一个率所属总体与理论率或期望率是否一致的假设检验 即某一个样本率是否符合一预设的总体率其中 由于一般计算率的样本都比较大 因此可使用u test进行检验 与之相比较的两个临界值就是时取1 96 时取2 58即当时 差异不显著 时 差异显著 时 差异极显著 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 43 例 某鱼药厂称该厂生产的治疗鱼烂鳃病的某种新药其疗效达85 今有一养鱼场购买了这种药 500尾病鱼经治疗后成活了415尾 该厂关于这一新药的宣传是否属实 设立无效假设 接受无效假设 即该厂生产的这种新药其疗效在95 的水平上可以得到承认 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 44 又例 孟德尔的豌豆花试验中 孟德尔曾假设其红花和白花的比例是3 1 今一个试验中 得红花447朵 白花165朵 问 这次试验的结果符合这一假设吗 作无效假设 这一次试验在95 的水平上可被接受 即本次试验的结果是符合孟德尔假设的 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 45 二 两个样本率差异的显著性检验这一类型的假设检验往往是检验两个样本率和所属总体率和是否一致 或者说两个样本率是否来自一个总体设两个样本容量分别是和 两样本中发生某一阳性事件的次数分别是和 则两样本率分别是和这两个样本率所属总体的率分别是和我们通过两样本率的差来推断和是否相等 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 46 这一类型的无效假设是其备择假设是用以检验的公式为 其中 由于一般样本很大 因此可用u test得到的u值与1 96和2 58相比时 差异不显著时 差异显著时 差异极显著 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 47 例 试验某种新药对杀灭鱼体外寄生虫的效果 用常规药物作对比试验 常规药物施与860头虫体 死亡585头 该新药施与920虫体 死亡672头 新药的疗效如何 新药的杀灭率为常规药物的杀灭率为设立假设 设 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 48 否定无效假设 即新药的除虫效果性状好于常规药物 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 49 又例 试验用EM制剂处理鱼苗可否提高鱼苗的成活率 处理组 施用EM 382尾鱼苗成活309尾 对照组 未施用任何药剂 278尾成活204尾 EM有效果吗 处理组成活率 对照组成活率 设 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 50 否定无效假设 接受备择假设 即EM可以显著提高鱼苗的成活率 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 51 三 小样本率假设检验的校正质量性状的率的计算 一般应使用大样本 但有些特殊情况下 样本不够大 只能得到一个小样本的率 如n 25 且np 5时 资料服从二项分布 不能再使用常规的u test 这时应进行校正 近似地进行u test 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 52 a 单个小样本率的校正性检验其公式是 b 两个小样本率的校正性比较其公式是 其中 是较大的样本率 是较小的样本率 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 53 第四节总体平均值的 小样本资料 区间估计 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 54 复习 标准误的作用复习 置信区间 置信度 置信半径 置信上限 置信下限当总体方差为已知时 或写成 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 55 当总体方差未知 但样本较大时 或写成 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 56 当总体方差未知 且样本较小时 总体平均值的置信区间估计应当用以下公式进行 或写成 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 57 第五节百分率的置信区间估计 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 58 一般情况下 率的估计都是大样本 因此 不管总体方差是否已知 我们都可以用下式来估计总体率的置信区间 或写成 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 59 从以上置信区间的设置 我们可以看出 置信区间的宽窄与置信度 样本容量 标准误等有关 因此 为了得到一个较好的置信区间用以估计总体平均值或总体率应当选择一个合理的置信度应当扩大样本容量应当提高试验的技术和精确度以降低试验误差 书山有路勤为径 学海无涯苦作舟专业分享 敬请收藏 60 思考与习题 1 某渔场进了一批用作饲料原料的鱼粉 该批鱼粉标定