主要内容

sampsizepwr

样品大小和测试能力

描述

sampsizepwr考虑到另一个两个值,计算假设测试的样本大小,电源或备选参数值。例如,考虑到替代假设的参数值,可以计算获得特定功率所需的示例大小。

nout= sampsizepwr (testtypep0p1返回样本大小,nout的规定类型的双面试验所需的testtype当具有0.90时的功率(当拒绝零假设时拒绝零假设的概率)0.90时,当显着性水平(当拒绝零假设时拒绝零假设的概率)是0.05。p0指定null假设下的参数值。p1指定在备择假设下测试的单个参数的值或值的数组。

例子

nout= sampsizepwr (testtypep0p1PWR.返回样本大小,nout,这对应于指定的电源,PWR.,和备择假设下的参数值,p1

例子

pwrout= sampsizepwr (testtypep0p1[],n返回为样本大小时所获得的功率n当true参数值为时p1

例子

P1OUT.= sampsizepwr (testtypep0[],PWR.n通过指定的示例大小返回可检测的参数值,n,以及规定的权力,PWR.

例子

___= sampsizepwr (testtypep0p1PWR.n名称,值使用一个或多个名称-值对参数返回前面的任何参数。例如,您可以更改测试的显著性级别,或指定右尾或左尾测试。名称-值对可以以任何顺序出现,但必须从第6个参数位置开始。

例子

全部折叠

公司运行制造过程,填充带100毫升液体的空瓶子。为了监控质量,公司随机选择几瓶,测量内部的液体体积。

确定公司必须使用的样本大小t-test检测100 mL和102 mL之间的差异,功率为0.80。假设标准偏差为5ml。

nout = sampsizepwr(“t”,[100 5],102,0.80)
nout = 52

为了检测100毫升和102毫升的平均容积差,必须对52瓶进行0.80倍的测试。

生成一个功率曲线来可视化样本大小如何影响测试的功率。

nn = 1:10 0;pwrout = sampsizepwr (“t”, 102年100年[5],[],nn);图;情节(nn、pwrout“b -”,nout,0.8,“罗”) 标题(“权力与样本量”)包含(“样本”) ylabel ('力量'

图中包含一个坐标轴。标题为Power versus Sample Size的轴包含2个类型为line的对象。

一位员工想在她办公室附近买一栋房子。她决定不考虑任何平均早上通勤时间超过20分钟的房子。右边检验的零假设是H0: μ. = 20,备择假设为HA: μ. > 20。所选显著性水平为0.05。

要确定平均通勤时间,员工每天早上每天早上在高峰时段到她的办公室,员工每天早上一周内的时间,所以她的样品大小为5.她认为标准差, σ. ,等于5。

该员工认为25分钟的真实平均通勤时间与她的目标20分钟限制相差太大,所以她想要在真实平均通勤时间为25分钟时检测出显著差异。找出得出平均通勤时间不大于20分钟的错误结论的概率。

计算测试的功率,然后从1中减去电源以获得 β

power = sampsizepwr(“t”,[20 5],25,[],5,'尾巴'“对”);= 1 -乘方
β= 0.4203

β Value表示该员工错误推断早晨通勤时间不大于20分钟的概率为0.4203。

员工认为这个风险太高了,她希望得出错误结论的概率不超过0.01。计算员工必须进行的测试次数,以获得0.99的幂次。

nout = sampsizepwr(“t”,[20 5],25,0.99,[],'尾巴'“对”
nout = 18

结果表明,她必须从一个候选房子进行18次试驾才能达到这个功率级别。

该员工决定她只有时间进行10次试驾。她还认为得出错误结论的概率为0.05。计算在平均通勤时间中产生可检测差异的最小真参数值。

p1out = sampsizepwr(“t”,[20 5],[],0.95,10,'尾巴'“对”
p1out = 25.6532.

考虑到员工的目标功率水平和样本大小,她的测试检测到与至少25.6532分钟的平均通勤时间有显著差异。

计算样本容量,n,需要区分p= 0.30来自p= 0.36,使用功率为0.8的二项式测试。

napprox = sampsizepwr (“p”, 0.30, 0.36, 0.8)
警告:值N>200是近似值。将功率绘制成N的函数可以显示具有所需功率的较低的N个值。
napprox = 485

结果表明,0.8的功率需要样本大小为485.然而,该结果是近似的。

制作一个剧情,以查看是否更小n值提供了所需的0.8的幂。

nn = 1:50 0;pwrout = sampsizepwr (“p”, 0.3, 0.36, [], nn);nexact = min (nn (pwrout > = 0.8))
nexact = 462.
pwrout图绘制(神经网络,“b -”, (napprox nexact], pwrout ([napprox nexact]),“罗”) 网格

图中包含一个坐标轴。轴包含2个类型的型号。

结果表明,样本量为462的这个测试也提供了0.8的幂次。

一位农民希望测试两种不同类型肥料对豆类作物产量的影响。他目前使用肥料A,但相信肥料B可能会提高作物产量。由于肥料B比肥料A昂贵,而农民希望限制他在该实验中用肥料B治疗的计划数量。

农民在每个处理组中使用2:1的植物比例。他用A肥料测试了10株作物,用b肥料测试了5株作物。使用A肥料的平均产量为每株1.4 kg,标准差为0.2。施用B肥的平均产量为每株1.7公斤。检验的显著性水平为0.05。

计算测试的功率。

pwr = sampsizepwr(《终结者2》,[1.4 0.2],1.7,[],5,“比”2)
PWR = 0.7165.

农民希望将测试的力量增加到0.90。计算他必须用每种肥料治疗多少植物。

n = sampsizepwr(《终结者2》,[1.4 0.2],1.7,0.9,[])
n = 11

为了将测试的力量增加到0.90,农民必须用每种类型的肥料测试11种植物。

农民希望减少必须使用B肥料处理的植株数量,但将试验功率保持在0.90,并保持初始的各处理组植株比例为2:1

使用每个处理组中植物的比例为2:1,计算农民必须测试多少植物才能获得0.90的幂。使用前一次测试中获得的平均值和标准差值。

[n1out, n2out] = sampsizepwr (《终结者2》,[1.4,0.2],1.7,0.9,[],“比”2)
n1out = 8.
n2out = 16.

为了获得0.90的幂,农民必须用肥料a处理16株作物,用肥料B处理8株作物。

输入参数

全部折叠

测试类型,指定为以下之一。

  • “z”- - - - - -z-对已知标准差的正态分布数据进行检验。

  • “t”- - - - - -t- 对于具有未知标准偏差的正常分布数据。

  • 《终结者2》——两个示例集中t-对标准偏差未知且方差相等的正态分布数据的检验。

  • “var”- Chi-Square对正常分布数据的差异测试。

  • “p”-测试p二项分布的参数(成功概率)。的“p”测试是一种离散测试,对于这种测试,增加样本量并不总是增加功率。为n大于200的值,可能存在小于返回的值n还产生指定功率的值。

null假设下的参数值,指定为标量值或标量值的双元素数组。

  • 如果testtype“z”“t”,然后p0是一个两个元素数组[mu0,sigma0]分别在零假设下的平均值和标准偏差。

  • 如果testtype《终结者2》,然后p0是一个两个元素数组[mu0,sigma0]在零假设和可选假设下,第一个样本的均值和标准差。

  • 如果testtype“var”,然后p0是零假设下的方差。

  • 如果testtype“p”,然后p0的价值p在零假设下。

数据类型:|双倍的

可选假设下的参数值,指定为标量值或标量值数组。

  • 如果testtype“z”“t”,然后p1是替代假设下的平均值的价值。

  • 如果testtype《终结者2》,然后p1是替代假设下第二个样品的平均值的值。

  • 如果testtype“var”,然后p1为备择假设下的方差。

  • 如果testtype“p”,然后p1的价值p在备择假设下。

如果您指定p1作为一个数组,那么sampsizepwr返回一个数组noutpwrout这和p1

要返回备选参数值,P1OUT.,指定p1使用空括号([]),如语法描述

数据类型:|双倍的

测试的功率,指定为范围(0,1)中的标量值或范围(0,1)中的标量值数组。检验的力量是在特定显著性水平下,当备择假设为真时,拒绝原假设的概率。

如果您指定PWR.作为一个数组,那么sampsizepwr返回一个数组noutP1OUT.这和PWR.

要返回一个功率值,pwrout,指定PWR.使用空括号([]),如语法描述

数据类型:|双倍的

样本大小,指定为正整数值或正整数值数组。

如果testtype《终结者2》,然后sampsizepwr假设两个样本尺寸是相等的。对于不等的样本尺寸,指定n作为两个样本中较小的一个,并使用“比”名称-值对参数指示样本大小比例。例如,如果较小的样本量为5,较大的样本量为10,则指定nAs 5, and the“比”名称-值对为2。

如果您指定n作为一个数组,那么sampsizepwr返回一个数组pwroutP1OUT.这和n

数据类型:|双倍的

名称-值对的观点

指定可选的逗号分隔的对名称,值论点。姓名参数名和价值为对应值。姓名必须出现在引号内。可以以任意顺序指定多个名称和值对参数name1,value1,...,namen,valuen

例子:“阿尔法”,0.01,“尾巴”,“对”指定具有0.01显著性水平的右尾检验。

测试的意义价值,指定为逗号分隔对组成“α”和范围(0,1)的标量值。

例子:“阿尔法”,0.01

数据类型:|双倍的

两个样本的样本尺寸比t-test,指定为逗号分隔的对“比”一个大于等于1的标量。的价值比率等于n2 / n1, 在哪里n2是否样本量较大,并且n1是较小的样本量。

为了归还能量,pwrout,或可选参数值,P1OUT.,指定两个样本尺寸的较小n和使用“比”表示样本量比例。

例子:'比例',2

测试类型,指定为逗号分隔的对,由'尾巴'以及以下其中之一:

  • “两个”-选择项不等于的双边检验p0

  • “对”- 单面测试比替代方案大于p0

  • “左”-单侧检测小于p0

例子:'尾巴','对'

输出参数

全部折叠

样本大小,作为正整数值或正整数值数组返回。

如果testtypet2,你用“比”名称值对参数指定两个不等样本大小的比率,然后nout返回两个样本大小中较小的一个。

或者,要返回两个样本大小,请将此参数指定为[n1out, n2out].在这种情况下,sampsizepwr返回较小的样本大小为n1out,以及更大的样本大小n2out

如果您指定PWR.p1作为一个数组,那么sampsizepwr返回一个数组nout这和PWR.p1

测试获得的功率,作为范围(0,1)中的标量值或作为范围(0,1)中的标量值数组返回。

如果您指定np1作为一个数组,那么sampsizepwr返回一个数组pwrout这和np1

替代假设的参数值,返回为标量值或作为标量值数组。

当计算P1OUT.“p”测试,如果在给定的空假设和显著性级别下没有可拒绝的选项,函数将显示一条警告消息并返回

介绍了R2006b