蒲公英 - 制药技术的传播者 GMP理论的实践者

搜索
查看: 5099|回复: 5
收起左侧

[统计应用] 均值的多重比较

[复制链接]
大师
发表于 2018-10-17 21:41:01 | 显示全部楼层 |阅读模式

欢迎您注册蒲公英

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
近期工作中,碰到更换一家供应商,做3批产品,然后做加速,加速结束,需要和原供应商批做对比分析,以判定新供应商与原供应商所生产批次在含量项上,是否相同。拿到4批数据,怎么分析,是两两进行t检验?似乎不太对,于是查,终于找到,应当用“均值的多重比较中的dunnett(邓尼特)法”来做,更适合。现将从minitab帮助中的介绍摘录如下,分享给大家。

通过均值的多重比较可以检查哪些均值不同,并估计它们的差异有多大。在拟合线性模型后,如果您使用统计 > 方差分析表 > 一般线性模型 > 比较,则可以获得均值的多重比较。
在使用多重比较时,必须进行以下选择:
·    配对比较或与对照组进行比较
·    要比较的均值
·    比较方法
·    如何显示结果
配对比较或与对照组进行比较在您没有对照水平,并且想要比较所有均值组合时,请选择配对。选择对照组以将水平均值与对照组均值进行比较。此方法合适后,使用配对比较会缺乏效率,因为配对置信区间更宽,并且假设检验对于给定置信水平而言功效更低。
要比较的均值使用多重比较时,选择要比较的均值是很重要的考虑事项;如果选择不合适,会导致置信水平偏离您的预期。您可能需要考虑以下问题:·    是应该只对具有显著 F 检验的项比较均值还是应该对差异看似很大的
均值集比较均值?·    比较均值应达到设计中的多深程度:仅限于每个因子内,或第一级交互作用项的每个组合内,还是要跨更高级别交互作用项的多个组合?最好是在收集数据之前先确定要比较的均值。如果只比较差异看似很大的均值(这称为数据探测法),则会增大看似不存在任何差异但结果却表明确实存在差异的可能性 [15],[28]。同样,如果将应用多重比较的条件确定为获得显著 F 检验,则多重比较的误差率会比无条件应用多重比较的误差率更高[15],[23]。 但在实际应用中,许多人通常还是将 F 检验作为选择要比较的均值的依据。方差分析 F 检验和多重比较并不是完全独立的评估。例如,如果 F 检验的 p 值为 0.9,则通过多重比较可能无法发现均值之间在统计上有显著的差异。
应将均值比较进行到设计内的多深程度?这需要进行权衡:如果在所有两因子组合水平上比较均值,并且高阶表现得更显著,则比较的均值可能是各种效应混合的结果;如果在太深的水平比较均值,则会由于样本数量太小且比
较次数太多而使比较变得无效。实际上,您可能会决定在您认为有意义的交互作用项的因子水平组合上比较均值。Minitab 会将可以比较均值的项限定于固定项或固定项之间的交互作用项。嵌套被视为是一种交互作用。
如果有两个因子分别为 A 和 B,并指定 A B,则 Minitab 将在每个因子内进行多重比较。如果您指定交互作用 (A * B),Minitab 将显示因子 A 和 B 的所有水平组合的多重比较。
多重比较方法
根据以下条件选择比较过程:您要比较的组均值、您要指定的置信水平类型和您要让结果的保守程度。此上下文中的“保守”指示真正的置信水平可能大于显示的置信水平。
除了 Fisher 方法,多重比较方法也内置了防误报功能。通过多重比较的防误报功能,置信区间的宽度比没有设置保护时的置信区间更宽。现将各种多重比较方法的一些特性总结如下:


比较方法      属性           您指定的置信水平
Tukey    所有配对比较,不保守         同时
Fisher   由于多重比较,没有防误判功能     单个
Dunnett     与对照组进行比较,不保守            同时
Bonferroni      最保守的方法                    同时
Sidak       保守,但保守性略低于Bonferroni      同时

如何显示结果
Minitab 以分组信息表、置信区间和假设检验形式显示多重比较的结果。 分组信息表包含用于分组因子水平的字母列。共享一个字母的水平不存在显著差异。相反,如果它们不共享一个字母,均值则存在显著差异。
查看置信区间时,除了评估统计上的显著性之外,还可以评估均值之间差异的实际显著性。通常,在置信区间不包含零时,您会否定均值之间无差异的原假设。
除了 Fisher 方法,Minitab 会计算假设检验统计量的调整 p 值。在假设集合内对某特定假设的调整 p 值是否定此特定假设所在的最小全族 a- 水平。

回复

使用道具 举报

大师
 楼主| 发表于 2018-10-17 21:42:16 | 显示全部楼层
多重比较的 Dunnett 法
在方差分析中用于为每个因子水平的均值与控制组均值之间的差异创建置信区间。如果区间包含 0,则两个比较的均值之间不存在显著差异。可以为所有比较指定整体误差率,Dunnett 法可以针对每个单个比较相应地确定置信水平。

例如,您正在研究三种减肥药丸,从而确定它们与无效对照剂之间是否有显著差异。在双盲试验中,50 个人服用药丸 A,50 个人服用药丸 B,50 个人服用药丸 C,还有 50 个人服用无效对照剂。服用无效对照剂的组被指定为控制组。您记录了每组的平均减肥情况,并使用 Dunnett 法执行方差分析,从而确定三种药丸中的任何一种产生的减肥效果与无效对照剂是否有显著差异。Dunnett 法产生三个置信区间:一个反映组 A 与服用无效对照剂的组之间的平均减肥效果的差异,一个反映组 B 与服用无效对照剂的组之间的平均减肥效果的差异,还有一个反映组 C 与服用无效对照剂的组之间的平均减肥效果的差异。可以将所有三个比较的整体误差率都设置为 0.10,这样一来,所有同时比较的置信水平都为 90%。

药丸 A 与无效对照剂之间差异的置信区间包含 0;因此,您断定组 A 与服用无效对照剂的组的减肥效果之间没有差异。药丸 B 与无效对照剂之间差异的置信区间仅包含负数;因此,您断定组 B 中的研究对象与服用无效对照剂的组中的研究对象相比,减肥效果较差。换句话说,药丸 B 与减肥效果有一定冲突。最后,药丸 C 与无效对照剂之间差异的置信区间仅包含正数;因此,您断定药丸 C 与无效对照剂相比,产生的减肥效果更明显。作为此项研究的结果,您推荐药丸 C。
回复

使用道具 举报

大师
 楼主| 发表于 2018-10-17 21:44:52 | 显示全部楼层
例子: 2.png    1.png
回复

使用道具 举报

药士
发表于 2018-10-17 23:24:40 | 显示全部楼层
学习一下统计分析
回复

使用道具 举报

药徒
发表于 2018-10-18 11:24:39 | 显示全部楼层
P值都是大于0.05,说明没有显著的差异,说明两家供应商的产品差不多,可以接受新供应商

回复

使用道具 举报

药生
发表于 2018-10-20 09:40:53 | 显示全部楼层
比较4个均值,需要单独进行6[4(4-1)/2)]次t检验。因此,犯类型I错误概率=1-(1-0.05)^6=0.26, 大于0.05。
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

×发帖声明
1、本站为技术交流论坛,发帖的内容具有互动属性。您在本站发布的内容:
①在无人回复的情况下,可以通过自助删帖功能随时删除(自助删帖功能关闭期间,可以联系管理员微信:8542508 处理。)
②在有人回复和讨论的情况下,主题帖和回复内容已构成一个不可分割的整体,您将不能直接删除该帖。
2、禁止发布任何涉政、涉黄赌毒及其他违反国家相关法律、法规、及本站版规的内容,详情请参阅《蒲公英论坛总版规》。
3、您在本站发表、转载的任何作品仅代表您个人观点,不代表本站观点。不要盗用有版权要求的作品,转贴请注明来源,否则文责自负。
4、请认真阅读上述条款,您发帖即代表接受上述条款。

QQ|手机版|蒲公英|ouryao|蒲公英 ( 京ICP备14042168号-1 )  增值电信业务经营许可证编号:京B2-20243455  互联网药品信息服务资格证书编号:(京)-非经营性-2024-0033

GMT+8, 2025-4-4 03:51

Powered by Discuz! X3.4运维单位:苏州豚鼠科技有限公司

Copyright © 2001-2020, Tencent Cloud.

声明:蒲公英网站所涉及的原创文章、文字内容、视频图片及首发资料,版权归作者及蒲公英网站所有,转载要在显著位置标明来源“蒲公英”;禁止任何形式的商业用途。违反上述声明的,本站及作者将追究法律责任。
快速回复 返回顶部 返回列表