灵敏度与特异度关联条件下Youden指数的样本量估计及meta分析方法

来源 :南方医科大学 | 被引量 : 9次 | 上传用户:zj888666
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
研究背景:诊断试验(diagnostic test)的评价方法始终是生物统计研究领域的热点。Youden指数是诊断试验评价领域应用最多的综合性评价指标之一。然而,最早提出的Youden指数并未考虑灵敏度和特异度关联的真实情况,亦未有配对设计的对应方法。为此,本人和导师等于2015年提出了基于灵敏度与特异度关联的Youden指数统计推断方法(以下简称“改进Youden指数方法”),并提出了配对诊断试验Youden指数的统计推断方法,解决了Youden指数应用中存在的上述问题。新的统计推断方法较原有方法对Ⅰ类错误的控制更好,同时具有更高的检验效能,也更符合应用中的实际。目前,我们所提出的方法对应的样本量估计方法和meta分析方法尚未见报道。这两类方法的产生无疑对诊断试验应用领域具有重要意义,对诊断试验评价的方法学研究亦有重要贡献。研究目的:针对单样本、两独立样本和配对样本三种设计,在已经建立的改进Youden指数方法的基础上,提出相应的样本量估计方法。针对单样本、两独立样本和配对样本三种设计,在考虑是否存在异质性的条件下,建立改进Youden指数方法的meta分析方法。方法:本研究所有工作均在同一个诊断试验的灵敏度与特异度关联假设下进行,灵敏度和特异度均服从二项分布。本研究所有工作,均考虑诊断试验设计中单样本,两独立样本和配对样本三种设计。1.对于样本量估计,在前期工作的基础上,根据近似正态检验样本量估计的基本原理,通过数学推导,得到三种设计下样本量估计方法。通过Monte Carlo模拟,考察所提出样本量估计方法的检验效能,同时探讨配对设计中,因配对产生关联性对样本量估计的影响。2.对于合并效应量估计,在前期工作的基础上,考虑meta分析中纳入研究间存在异质性和不存在异质性两种情况。对于不存在异质性的情况,采用固定效应模型,并考虑以Mantel Haenszel法,最小二乘法,极大似然法等三种方法进行效应量合并;对于存在异质性的情况,采用随机效应模型,并考虑以DerSimonian-Laird法进行效应量合并。通过Monte Carlo模拟,考察不同模型下,不同加权方法所得合并效应量置信区间的覆盖率,研究其统计学性能。结果:1.改进Youden指数方法的样本量估计对单样本Youden指数(J),设灵敏度为Sen,特异度为Spe, Var(*)为方差,诊断试验中病例组样本量为n1,对照组样本量为n0,且有n0/n1=r,总样本量为N。同时设检验水准为α,Ⅱ类错误为β,检验效能为1-β,采用双侧检验(以下样本量估计的设定均仿此)。根据前期研究的结果,已知:Var(J)= Sen2Var(Spe)+Spe2Var(Sen)基于近似正态检验及其样本量、方差及Ⅰ/Ⅱ类错误的关系,得到病例组样本量的估计为:进而,由n0/n1=r可以得到对照组样本量的估计和总样本量的估计。对两独立样本诊断试验,设两个独立的Youden指数之差为DJ。其中,第i个诊断试验的Youden指数为Ji,病例组样本量为ni1,对照组样本量为ni0,总样本量为Ni。设,n10=r1n11,n21=r2n11,n20=r3n11。根据前期研究的结果,已知:基于近似正态检验及其样本量、方差及Ⅰ/Ⅱ类错误的关系,得到其中一个诊断试验病例组样本量的估计为:进而,根据两个诊断试验各病例组与对照组样本量的关系,可以得到各个诊断试验对照组样本量的估计和总样本量的估计。对配对样本诊断试验,设配对的Youden指数之差为DJ,诊断试验中病例组样本量为n1,对照组样本量为n0,且有n0/n1=r,总样本量为N。其中,第i个诊断试验的Youden指数为Ji。根据前期研究的结果,已知:基于近似正态检验及其样本量、方差及Ⅰ/Ⅱ类错误的关系,令:得到病例组样本量的估计为:进而由n0/n1=r可以得到对照组样本量的估计和总样本量的估计。Monte Carlo模拟的结果显示,基于改进Youden指数方法,本研究所提出的单样本、两独立样本和配对样本的样本量估计方法,其检验效能稳定,始终控制在设定的检验效能附近,样本量估计结果可靠,可以满足实际应用的要求。2. Youden指数meta分析方法在不存在异质性时,考虑固定效应模型,采用三种加权方法进行合并效应量估计。(1)基于Mantel Haenszel法(MH法)的合并效应量估计a.对单样本诊断试验,设meta分析纳入的第i个诊断试验的MH法权重为WMHi,其病例组样本量为ni1,对照组样本量为ni0,灵敏度为Seni,特异度为Spei, Youden指数为Ji,其合并效应量为JMH,则其权重:则基于MH法加权的单样本Youden指数的合并效应量为:其合并效应量的方差为:其(1-α)%置信区间为:b.对两独立样本诊断试验,设meta分析纳入的第i个研究的权重为WMHi,其中一个诊断试验的病例组样本量为ii1,对照组样本量为ni0,其灵敏度为Seni1,特异度为Spei1, Youden指数为Ji1,;另一个诊断试验的病例组样本量为mil,对照组样本量为mi0,其灵敏度为Seni2,特异度为Spei2; Youden指数为Ji2,两个Youden指数之差为DJi,其合并效应量为DMH,则其权重:则,基于MH法加权的两样本Youden指数之差的合并效应量为:其合并效应量的方差为:其(1-α)%置信区间为:c.对配对样本诊断试验,设meta分析纳入的第i个研究的权重为WMHi,病例组样本量为ni1,对照组样本量为ni0,其中一个诊断试验的灵敏度为Seni1,特异度为Spei1, Youden指数为Ji1;另一个诊断试验的灵敏度为Sena,特异度为Spei2; Youden指数为Ji2,两个Youden指数之差为DJi,其合并效应量为DMH,则其权重:其基于MH法加权的配对样本Youden指数的合并效应量为:其合并效应量的方差为:其(1-α)%置信区间为:(2)基于最小二乘法(逆方差法)法的合并效应量估计a.对单样本诊断试验,设meta分析纳入的第i个诊断试验的最小二乘法权重为WLSi,其病例组样本量为ni1,灵敏度为Seni,特异度为Spei, Youden指数为Ji,其合并效应量为JLS,则:基于最小二乘法加权的单样本Youden指数的合并效应量为:其合并效应量的方差为:其(1-α)%置信区间为:b.对两独立样本诊断试验,设meta分析纳入的第i个研究的最小二乘法权重为WLSi,其中一个诊断试验的灵敏度为Seni1,特异度为Spei1, Youden指数为Ji1;另一个诊断试验的灵敏度为Seni2,特异度为Spei2;Youden指数为Ji2,两个Youden指数之差为DJi,其合并效应量为DLS,则其权重:基于最小二乘法加权的两独立样本Youden指数差的合并效应量为:其合并效应量的方差为:其(1-α)%置信区间为:c.对配对样本诊断试验,设meta分析纳入的第i个研究的最小二乘法权重为WLSi,设其中一个诊断试验的灵敏度为Seni1,特异度为Spei1, Youden指数为Ji1,;另一个诊断试验的灵敏度为Seni2,特异度为Spei2; Youden指数为Ji2,两个Youden指数之差为DJi,其合并效应量为DLS,则:基于最小二乘法加权的配对样本Youden指数的合并效应量为:其合并效应量的方差为:其(1-α)%置信区间为:(3)基于极大似然法(ML法)的合并效应量估计方法极大似然法所得合并效应量及其方差的估计无显式。a.对单样本诊断试验,在固定效应模型下,设meta分析纳入的第i个诊断试验病(共k个)其灵敏度为Seni,特异度为Spei, Youden指数为Ji,其合并效应量为JML,建立对数似然方程:求偏导数,并令其值等于0后得到方程组:解上述方程组即可得到Youden指数的合并效应量及其方差的估计。根据对数似然函数的Fisher信息矩阵的逆矩阵,可以解出其合并效应量方差的估计值。b.对两独立样本和配对样本Youden指数的合并效应量,设meta分析纳入的第i个研究,其中一个诊断试验的灵敏度为Seni1,特异度为Spei1, Youden指数为Ji1,;另一个诊断试验的灵敏度为Seni2,特异度为Spei2; Youden指数为Ji2,两个Youden指数之差为DJi,其合并效应量为DML,建立对数似然方程:求偏导数后得到:对方程组求解即可得到合并效应量的极大似然估计值。根据对数似然函数的Fisher信息矩阵的逆矩阵,可以解出其合并效应量方差的估计值。虽然两独立样本和配对样本的方程是一致的,但是其Fisher信息矩阵截然不同。两独立样本设计下,微分算子中两个随机变量的下标不同时,其元素值为0,而对于配对样本设计,其Fisher信息矩阵所有元素均为非0元素。在存在异质性时,考虑随机效应模型,采用DerSimonian-Laird法(DL法)进行效应量估计。a.对单样本诊断试验,设meta分析纳入的第i个诊断试验在固定模型下最小二乘法(逆方差法)的权重为WLSi,其灵敏度为Seni,特异度为Spei, Youden指数为Ji,其合并效应量为JDS,设异质性方差为τ2,则,随机效应模型下的权重矿为:单样本Youden指数的合并效应量为:其方差Var(JDS)的估计为:其(1-α)%置信区间为:b.对两独立样本诊断试验,设meta分析纳入的第i个研究的在固定效应模型下最小二乘法权重为WLSi其灵敏度为Seni1,特异度为Spei1, Youden指数为Ji1,;另一个诊断试验的病例组样本量为mi1,对照组样本量为mi0,其灵敏度为Seni2,特异度为Spei2; Youden指数为Ji2,两个Youden指数之差为DJi,其在随机效应模型下的权重为矿,合并效应量为DDS,设异质性方差为产,则其权重为:随机效应模型下两独立样本Youden指数差的合并效应量为:其方差Var(DDs)的估计为:其(1-α)%置信区间为:c.对配对样本诊断试验,设meta分析纳入的第i个研究在固定效应模型下的最小二乘法权重为WLSi,其中一个诊断试验的灵敏度为Seni1,特异度为Spei1, Youden指数为Ji1,;另一个诊断试验的灵敏度为Seni2,特异度为Spei2; Youden指数为Ji2,两个Youden指数之差为DJi,其载随机效应模型下权重为矿,合并效应量DDS,则其权重为:随机效应模型下配对样本Youden指数差的合并效应量为:其方差Var(DDS)的估计为:其(1-α)%置信区间为:Monte Carlo模拟结果显示,在固定效应模型下,Mantel Haenszel法得到的置信区间覆盖率稳定在95%水平周围,效果最为理想;极大似然法得到的置信区间,覆盖率水平略低,但仍处在可以接受的水平内;最小二乘法得到的置信区间覆盖率不理想,且覆盖率随纳入研究数变化的趋势有问题。在随机效应模型下,DerSimonian-Laird法所得的置信区间覆盖率差,且随着纳入研究数的增加而下降。结论:本研究在灵敏度和特异度关联条件下,针对诊断试验的三种常见类型,建立了基于改进Youden指数方法的样本量估计方法和meta分析方法。结果表明,三种样本量估计方法均具有良好的统计性能,可以符合实践的要求。对于Meta分析效应量合并方法,在纳入研究间不存在异质性时,基于Mantel Haenszel法或者极大似然法得到的合并效应量,统计性能良好,均可满足实际应用的要求;在纳入研究存在异质性时,本研究所提出效应量合并方法统计性尚不理想,有待深入研究做进一步研究。
其他文献
摘要:基于研究主导权与研究能力的缺乏,人文社会科学研究者往往面临着各种研究伦理风险。现有研究主要关注实证研究风险与结果呈现风险两个方面,忽视了作为“总开关”的研究问题风险与极具科学性的研究设计风险。文章应用文献荟萃分析法,基于过程视角将研究伦理风险细分为研究问题风险、研究设计风险、实证研究风险和结果公开风险,根据意识与行为、局部與整体的内在关联将其归因于研究者缺乏风险防控意识,受到学术利益驱使,缺
斗杆是挖掘机式冲击器的主要承载部件之一,工作中斗杆结构处于较高频率振动状态,并承受较大振幅的动载荷,很容易出现疲劳破损的现象。为了分析斗杆结构在冲击工作时的振动特
科学认识和准确概括恩格斯思想的历史地位与贡献一直是国内外学界争论不休的话题。为避免出现拔高或贬低恩格斯思想的两种倾向,我们应对其思想进行整体性定位,既注重恩格斯思
广州市耳鼻咽喉头颈外科医院(广州市第十二人民医院)将于2012年3月30日-4月1日在广州市举办国家级继续医学教育项目“激光与微创治疗在耳鼻咽喉头颈外科应用的新进展学习班”,编号:2012.07-01.115(国)。届时,学习班将邀请多名国内外著名专家通过专题讲座、手术演示、手术录像、专题讨论等方式进行教学。学习结束授予I类继续教育学分10分。真诚欢迎相关专业的专家和同道们踊跃参加!
理顺中央和地方财政关系对于充分发挥中央、地方两个积极性,提升国家治理能力具有重要意义。首先回顾了1978年改革开放以来中央和地方财政关系演进的三个阶段,深入剖析了当前
回 回 产卜爹仇贱回——回 日E回。”。回祖 一回“。回干 肉果幻中 N_。NH lP7-ewwe--一”$ MN。W;- __._——————》 砧叫]们羽 制作:陈恬’#陈川个美食 Back to yield
随着我国经济的快速发展,信用卡这一便捷支付型电子卡得到了广泛的推广与使用。但这一现象的背后则是相关犯罪行为的频发,可以说是对人民财产利益与国家金融秩序的严重侵害。
摘要:学术期刊是知识生产的一个特殊部门,它不仅是学者发表成果和学术交流的载体,更承担着学术评价的功能与责任。期刊评价是学术评价的重要组成部分,是推进学术评价伦理化、合理化的路径之一。学术期刊的社会和公共属性决定了期刊评价具有深刻的政治性、行政性、社会性和情感性。期刊评价中的定性和定量评价法是质化和量化研究在评价中的应用,主要体现为文献计量法的影响因子类评价和同行评议的定性评价。文献计量法重视对期刊
回 回 产卜爹仇贱回——回 日E回。”。回祖 一回“。回干 肉果幻中 N_。NH lP7-ewwe--一”$ MN。W;- __._——————》 砧叫]们羽 制作:陈恬’#陈川个美食 Back to yield
目的 研究产后胎盘植入诊断中超声造影的临床应用效果。方法 随机选取2015年12月至2017年12月本院收治的疑似产后胎盘植入患者85例,具实施超声造影检查,详细观察病灶的边界、