stata描述性统计结果怎么分析（Stata的操作和结果解读）

云生信学生物信息学 原创不易，请多多关注和转发支持，谢谢！

以下文章来源于尔云间meta分析 ，作者尔云间

尔云间meta分析.

分享meta分析学习方法、医学统计基础知识、meta分析软件操作技巧、SCI论文发表经验等精彩内容。

剪补法基于发表性偏倚造成漏斗图不对称这一假设，采用迭代方法估计缺失研究的数量，其意义并不在于估计出缺失研究的具体数目，在添补一部分研究后，重新进行Meta分析。

如果合并效应量估计值变化不明显，说明发表性偏倚影响不大，结果比较稳健。因此它实际上是一种敏感性分析方法。

剪补法可以通过Stata软件的metatrim命令来实现，下面将介绍剪补法在Stata中实现的过程以及结果解读（有视频教程需要可联系小编）。

1. 数据的读取

采用Stata软件自带的数据，在命令窗口中输入并运行以下代码，可以直接读取数据：

use http://fmwww.bc.edu/repec/bocode/m/metan_example_data

或者运行help metan，打开metan菜单的帮助文件，找到上述链接，点击打开数据。

部分数据截图

2. 发表偏倚检验

首先对数据做发表偏倚检验：

gen logOR=log( OR)

gen selogOR=(log( ORuci)-log( ORlci))/2/1.96

metabias logOR selogOR, egger

结果显示，存在发表偏倚(Egger's test, t=2.29, P=0.034<0.05)，因此有必要采用剪补法评估合并结果的稳定性。

3. 剪补法的操作

metatrim命令格式可以有多种实现形式，可以在命令框中输入help metatrim查看实现的方法。也可参照下图所示，选择对应的参数完成剪补法。

4. 结果解读

第一部分，用metan命令进行Meta分析(合并效应量为logOR)。同时汇报了固定效应模型和随机效应模型的结果，异质性检验，Q= 32.689, P=0.026，采用随机效应模型，效应指标合并结果为logOR=0.106，95%可信区间为(-0.084, 0.296)。

第二部分，显示以随机效应模型，采用Linear法经过4次迭代后，软件估计缺失的研究数目，结果为5。

软件填补了5个虚拟的研究

最后，纳入5篇虚拟研究的数据之后，重新对所有研究进行Meta分析，结果显示异质性检验：Q= 52.453, p= 0.001，采用随机效应模型，所得效应指标合并结果为lOR=0.967，95%CI为(0.790, 1.185)。

也就是说补上5篇研究之后，结果依然没有统计学意义，没有发生逆转，因此合并结果稳健。

Stata的meta分析操作，我们已经分享了很多教程，你学会了吗？如果你还有操作上的问题，联系小编帮你解决。