Meta分析中失安全系数意义的探讨

《Meta分析中失安全系数意义的探讨》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、万方数据·428·Meta分析中失安全系数意义的探讨李宪1王洪源2Meta分析是一类用于比较和综合针对同一科学问题所取得研究结果的统计方法。随着循证医学的发展，Meta分析得到了广泛的应用。在PubMed收录的Meta分析文献已经从1990年的200多篇增加到2005年后的2000多篇【lJ。在我国也同样如此，在CNKI的中国期刊全文数据库中2007年收录以Meta分析为主题的文献已经达到380多篇。由于Meta分析属于观察性研究，在分析时需要进行出版偏性(publicationbias)的估计，这是因为通常情况下有统计学意义的研究结果较无统计学意义(或无效)的研究结果被

2、报告和／或发表的可能性更大，呈现在读者面前的论文大都为“有统计学意义”的论文。因此识别是否存在出版偏性成为了Meta分析的一个重要问题。识别出版偏性可以使用观察漏斗图的对称性、用Begg’s的秩相关法及Egger’s回归法进行出版偏性的检验。1979年Rosenthal【2】提出失安全系数(fail．safeN，Nfs)的概念，人们也用它做为评估出版偏倚的指标，他认为当Meta分析结果有统计学意义时，为排除可能的出版偏倚，可计算最少需要多少个未发表的研究才能使Meta分析的结论逆转。然而也有学者认为用失安全系数评估发表偏倚存在问题【3】，我们认为这种观点是合理的，本文将通

3、过一个病例对照研究的Meta分析随机模拟的样本阐述这种观点。资料与方法i．资料来源采用随机模拟的方法生成病例对照研究Meta分析的研究样本，根据一般研究样本量的分布，确定本研究包括46个独立的病例对照研究，每个研究病例与对照的例数比例均为1：1，在本研究中样本量在30～48、50～98、100～198的各10个研究，样本量在200～498的8个研究，样本量在500～998的5个研究，样本量在1000～2998的3个研究。具体的样本量由程序随机产生。确定本研究总体的病例组暴露率为24％。对照组的暴露率为20％，按各研究的样本量分别从研究总体中随机抽取得到研究样本。本研究采用

4、＆蟠9．1．3软件随机产生研究样本。2．统计方法1．北京积水潭医院病案统计科(100035)2．北京大学公共卫生学院流行病学与卫生统计学系(100083)由于OR近似服从对数正态分布，所以Meta分析时采用1nOR做为效应值，Meta分析时所用的计算公式见文献[4]。分别计算各独立研究效应值是否为“O”的统计量U，根据u值进行分段，删除部分无统计意义的研究，并计算相应的失安全系数【5】：fFTT＼2№0．052【揣J一惫其中u为每一个独立研究效应值是否为“0”检验的U统计量，忌为已收集独立研究的个数。描述失安全系数与未发表的研究个数的关系。使用sTATA软件进行Meta分

5、析。结果1．模拟数据的描述模拟研究的汇总结果见表1。2．对模拟数据各研究的效应值进行齐性检验并计算合并的效应值首先对46个研究的效应值进行齐性检验，Q=42．634，自由度为45，P=0．573，说明46个研究的效应值来自同一总体，在计算合并的效应值时可以采用固定效应模型。46个研究的合并效应值为0．271，其95％置信区间为(0。199，O．343)，研究总体的1nOR=0．2336，在置信区间的范围内。3．对模拟数据进行出版偏性的估计以每个研究的效应值为横坐标，以每个研究效应值标准差为纵坐标绘制Funnelplot，可以看出Funnelplot的对称性较好，说明本研究

6、收集的数据不存在明显的出版偏性。结果见图1。采用Begg’sTest对46个研究的效应值进行检验(P≥0．514)，没有发现明显的出版偏性。采用Egger’st鹤t对46个研究的效应值进行检验(P≥O．192)，没有发现明显的出版偏性。4．失安全系数(fai卜髓fe舢mber，Nfs)与出版偏性的关系假定每个研究在U<1．96，即P>0．05时可能得不到发表。按【厂统计量进行分段，在不同情况下计算失安全系数。从表3中可以看出如果只有U统计量≥1．96时研究文献才能够发表，46个研究中只发表了9个，有37个没有被发表，9个被发表的研究U统万方数据·429·计量的合计为22．

7、3346，此时计算的N必o．05为此时计算的№o．05为715．74。176．47。如果只有u统计量≥1时研究文献才能够当所有的46个研究都被发表，u统计量的合计发表，46个研究中只发表了23个，有23个没有被发为45．7774，此时计算的晰o．05为733．14。表，23个被发表的研究U统计量的合计为44．5749，表1模拟研究的汇总数据模拟研究编号lIloR的标准差Z2统计量U统计量13210611537．5031．251．3200O．2776O．7467O．13850．37222288105391242027．0813．