样本量估算
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
1.单因素二水平设计定量资料的非劣效性检验时样本量的估算
1.1计算公式:
非劣效性检验应当采用单侧的检验水准α,假定允许的第二类错误概率不超过β,则非劣效性检验每组需要的样本含量为:
22211)/()(2θδβα-+=--L S u u n (1-1)[1]
2221)
/()(2δβαe s z z n n ⨯+==(1-2)[2] 1.2式中各参数代表的意义,n 为每组样本含量,α-1u 、β-1u 为单侧标准正态离差界值,S 为估计的共同标准差,L δ为非劣界值,且L δ<0,θ为试验组与对照组总体均值差值的估计值。
说明:单因素二水平设计定量资料的非劣效性检验时样本量的估算公式与上式完全类似,只需将非劣界值L δ(L δ<0)替换成优效界值u δ(u δ>0)即可。
1.3例题:某利尿新药拟进行Ⅱ期临床试验,与阳性药按1:1的比例安排例数,考察24h 新药利尿量(ml )是否不差于阳性药。根据以往的疗效和统计学的一般要求,取α=0.05,β=0.20,非劣效界值L δ=﹣60ml ,已知两组共同标准差S =180ml ,假定新药与阳性对照药总体利尿量的差值θ=﹣20ml ,问每组需要多少病例?
将05.01-u =1.645,20.01-u =0.845,s=180,L δ=﹣60,θ=﹣20代入公式,得:
22211)/()(2θδβα-+=--L S u u n =2
(1.645+0.845)2×1802/(﹣60﹣(﹣20))2≈251.1,取n=252,即每组需要252例。
2.单因素二水平设计定性资料的非劣效性检验时样本含量的估算
2.1计算公式:
非劣效性检验应当采用单侧检验,检验水准为α,假定允许的第二类错误概率不超过β,试验组与对照组总体率的差值为C T ππθ-=
(T π、C π未知时可用样本频率估计),两组的平均有效率为2/)(C T πππ+=,非劣界值为u δ<0,则在两组样本含量相等的情况下,非劣效性检验每组需要的样本含量为:
2211)/()1()(2θδππβα--+=--L u u n (2-1)[1]
2合合221/)-1()(2δβαp p z z n n +==(2-2)[2]
说明:单因素二水平设计定性资料的优效性检验时样本含量的估计公式与式(2-1)完全类似,只需将非劣界值L δ(L δ<0)替换成优效界值u δ(u δ>0)即可。
2.2例题:某新药拟进行Ⅱ期临床试验,与阳性药按1:1的比例安排例数,考察新药临床治愈率是否不差于阳性药。根据以往的疗效和统计学的一般要求,取α=0.05,β=0.20,非劣界值δ=﹣0.15,平均有效率P=0.80,并假定两组总体有效率相等,问每组需
要多少病例?
代入公式(2-1)或者(2-2)得n=2(1.645+0.845)2×0.8×(1-0.8)/(﹣0.15)2≈88,即每组需要88例。
3.单因素二水平设计定量资料的等效性检验时样本含量的估算
3.1计算公式
对于等效性检验,按照双侧的检验水准α(若等效性检验采用双单侧检验,则没个单侧检验的检验水准为α/2),允许犯第二类错误的概率不超过β,则结果变量为定量变量时等效性检验每组需要的样本含量的计算公式(3-1)[1]为:
1)1(])[()1(])[(2/1222/122-⎪⎭
⎪⎬⎫⎪⎩⎪⎨⎧-++-+⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧-+--=--αασδμμφσδμμφu k kn u k kn power C T C T C T C T 这是一个超越方程。可通过迭代法求解出满足要求的试验组样本含量T n 。式中T μ、C μ分别为试验组与对照组的总体均值,未知时可用样本均值代替;δ是等效界值,且δ>0;2/1α-u 为标准正态离差界值;T n 为试验组的样本含量;k 为对照组与试验组样本含量的比值,即T C n n /=k ,或者C n =k T n ;2C σ是试验组与对照组的合并方差,未知时可用样本的合并方差2C S 代替。power 不小于1-β。
3.2例题:某利尿新药拟进行Ⅱ期临床试验,与阳性药按1:1的比例安排例数,考察24h 新药利尿量(ml )是否相当于阳性药。根据以往的疗效和统计学的一般要求,取α=0.05,β=0.20,非劣效界值 =﹣60ml ,已知两组共同标准差S=180ml ,假定新药与阳性对照药总体利尿量的差值θ=﹣20ml ,问每组需要多少病例?
k=1,T n =100,代入公式power=Φ(1.183)+Φ(-0.389)-1≈0.8816+0.3483-1=0.2299<0.80,所以T n =100过小,需继续增大样本含量进行迭代尝试……当T n =318时得power=0.8014>0.80.故试验组需要319例,对照组也需要319例。
4.单因素二水平设计定量资料的等效性检验时样本含量的估算
4.1计算公式(4-1)[1]:
1)1()1(])[()1()1(])[(2/122/12-⎪⎭
⎪⎬⎫⎪⎩⎪⎨⎧--++-+⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧--+--=--ααππδππφππδππφu k kn u k kn power T C T T C T 公式中的各参数:双侧检验水准α,允许犯第二类错误的概率不超过β,试验组与对照组总体率的差值为θ=πT -πC (πT 、πC 未知时可用样本频率估计),两组的平均有效率为π=(πT +πC )/2,等效界值为δ>0,对照组与试验组样本含量的比值为k ,即T C n n /=k ,或者C n =k T n 。power 不小于1-β。
例题:某新药拟进行Ⅱ期临床试验,与阳性药按1:1的比例安排例数,考察新药临床治愈率是否相当于阳性药。根据以往的疗效和统计学的一般要求,取α=0.05,β=0.20,非劣界值δ=﹣0.15,平均有效率P=0.80,并假定两组总体有效率相等,问每组需要多