功效与样本量的计算

今日头条

1151人已加入

描述

我经常被问到这样的问题:“做假设检验时,需要的样本量是多少”, “我的实验究竟需要多少样本才有意义呢?”,这类问题可以通过功效与样本量计算来解决。

什么是功效?

所有检验都不可能尽善尽美,总存在这样的可能,即当原假设 (H0) 实际为真时您会根据检验结果否定原假设(I 型错误),或当原假设实际为假时您根据检验结果未能否定 H0(II 型错误)。这是因为为了估计总体均值,您必须使用随机样本,而这些随机样本都是随机抽取的。

Minitab

第II类错误(被坑)的概率为β,而1- β我们称之为功效。II类错误是如果有个人有罪,而错误地判他无罪,所以“功效”就可以理解为,如果他实际上确实有罪,而且法官正确地判断他有罪的概率。

假设检验的功效受到检验的样本数量、差值、数据变异性和显著性水平的影响。

Minitab

案例分析

一种镇定剂中的活性成分的目标含量为365mg/ml,每批镇静剂在使用前都要进行检测,以发现明显偏离目标的情况。

数据收集:从6个不同批次中随机抽取6个样品,测量其中的有效成分含量。

Minitab

描述性统计分析

MinitabMinitabMinitab

注意:不能直接从描述性统计量中得到偏离目标含量365mg/ml的结论,即使样本均值确实不是365mg/ml.

单样本t检验

MinitabMinitabMinitab

单样本t检验的结果显示P>0.05,故不拒绝原假设,换句话说,没有足够的证据拒绝有效成分含量均值为365mg/ml。

功效与样本量的计算

在得到的结果是没有发现统计学差异时,报告一下你的检验功效很有必要。因为如果你的功效很低,则说明很可能是样本量不足导致无法检验出有显著差异的结果。

MinitabMinitab

在这里,我不得不对上面的对话框做一些解释(我亲爱的客户朋友经常对此有疑问)。

差值

差值是在检测总体参数的假设值与实际值时您想得到的最小差值,看不懂,没关系,我们举个例子。在上面的对话框中我输入的差值为-2.5和2.5,为什么呢?因为在这里我认为,如果有效成分含量偏离目标值2.5mg/ml时,我就认为是有显著偏离,当然这里的偏离可以偏大(2.5),也可以偏小(-2.5)。

标准差

在功效和样本数量分析中使用的标准差取决于您是否已收集了数据。

1) 尚未收集数据:请使用总体标准差的估计值。您的估计值应以相关研究、设计规范、试验研究、学科知识或相似信息为基础。

2) 如果您已经收集并分析了数据,则使用下表所示的分析结果所提供的值。

您执行的分析的类型 用于标准差的值
单样本 t 或单样本 Z 样本的标准差
双样本 t 合并标准差
配对 t 配对差的标准差
单因子方差分析、2 水平因子设计或 Plackett-Burman 设计 均方误 (MSE) 平方根的平均值

Minitab结果如下;

Minitab

结果显示,当样本量只有6个时,这时候的检验功效只有0.54,或者说犯第二类错误的概率高达0.46。功效太低,当前没有显著差异的结论很可能是样本量不足导致。

如果对于此分析,我想功效不低于0.9,那么这时候我该抽取多少样本量呢?

MinitabMinitab

Minitab结果显示,要想功效不低于0.9,那么这时候我们需要再抽6个样本,一共12个样本才行。

再来一次单样本t检验

现在,又再抽了6个样本,一共12个样本,再来做一次假设检验。

MinitabMinitabMinitab

这时候我们发现,当样本量增加到12时,P<0.05,故拒绝原假设,得到有显著偏离目标值的结论。

小结

检验所需样本量的多少取决于你愿意承担多大犯错误的风险,所以,当你以后分析纠结样本量够不够时,请尝试先做个功效和样本量的计算。

如果检验的功效低,则您可能无法检测到效应并错误地得出不存在任何效应的结论。另外,也要注意如果检验的功效太高,则很小的效应或不需要关注的效应看上去都可能会很显著。

思考

如何提高假设检验的功效?

审核编辑:符乾江

打开APP阅读更多精彩内容
声明:本文内容及配图由入驻作者撰写或者入驻合作网站授权转载。文章观点仅代表作者本人,不代表电子发烧友网立场。文章及其配图仅供工程师学习之用,如有内容侵权或者其他违规问题,请联系本站处理。 举报投诉

全部0条评论

快来发表一下你的评论吧 !

×
20
完善资料,
赚取积分