兰州seo整站优化服务商,活动推广,公众号开发信息开发者密码是什么,用什么制作网站目前不少文章用到了孟德尔随机化meta分析#xff0c;今天咱们也来介绍一下#xff0c;孟德尔随机化meta其实主要就是meta分析的过程#xff0c;提取了孟德尔随机化文章的结果#xff0c;实质上就是个meta分析#xff0c;不过多个孟德尔随机化随机化的结果合并更加加强了结…目前不少文章用到了孟德尔随机化meta分析今天咱们也来介绍一下孟德尔随机化meta其实主要就是meta分析的过程提取了孟德尔随机化文章的结果实质上就是个meta分析不过多个孟德尔随机化随机化的结果合并更加加强了结果的可靠性。有部分人可能对meta分析不是很了解咱们今天先来介绍一下meta分析基础为下一讲孟德尔随机化meta分析做准备。 R语言进行进行meta分析咱们就做最基本的部分就行不必搞得太复杂。咱们先导入数据和R包。很多R包都能做咱们随便选个meta包就可以了。
library(meta)
bc-read.csv(E:/r/test/senlintu1.csv,sep,,headerTRUE)
names(bc)这是一个很简单的数据stud表示研究名称a为实验组阳性人数b为实验组总人数c为对照组阳性人数d为对照组总人数公众号回复森林图数据1可以获得数据。 咱们先来看下函数格式event.e就是实验组阳性人数n.e实验组总人数event.c对照组阳性人数n.c对照组总人数data就是你的数据studlab填入其他的项目method这里选Inverse倒方差的方法就可以了sm这里填入结果类型如果你需要的是OR的结果就填入OR
metabin(event.e, n.e, event.c, n.c, data,studlab paste(), sm, method Inverse)先生成meta分析结果就按顺序把结果填进去就可以了
out-metabin(event.ea, n.eb,event.cc,n.cd,databc,smOR,studlab paste(study),method Inverse)直接填进去也是可以的这样方便点
out-metabin(a,b,c,d,databc,smOR,studlab paste(study),method Inverse)解析结果
summary(out)我简单解析一下这类教程网络多得是可以百度一下。先看异质性如果有异质性应需用随机效应模型否则可选用固定效应模型。我们这里选择随机效应模型。 我这里异质性I89.9%挺大的我们可以使用剔除法来观察剔除单个观察后的异质性我这里随便举例假如我剔除第一个 先加个ID
bc$id-1:13删除第一项我们可以使用亚组函数subset来控制使得subsetid2就可以了
out-metabin(a,b,c,d,databc,smOR,studlab paste(study),method Inverse,commonF,subsetid2)
summary(out)可以看到异质性较前下降了一点接下来就可以绘制森林图了
out-metabin(a,b,c,d,databc,smOR,studlab paste(study),method Inverse,commonF)forest(out)所有结果都可以在图上找到结果主要是看最后的综合结果我这里的OR是2.09.接下来做偏倚检查主要是漏斗图和Egger法 漏斗图比较理想的漏斗图应该是散点对称分布
funnel(out)Egger法P值大于0.05说明没有偏倚还给出了参考文献
metabias(out,method.biasEgger)这里还有一个peters法参考文献是篇JAMA的文章应该也是蛮靠谱的
metabias(out,method.biaspeters)这里简单的介绍了一下meta分析的操作肯定没有别人专门将meta分析的详细如果想进行孟德尔随机化meta分析最好去补一补meta分析的知识。
