第四篇公众号：来自微信 天桥下的卖艺者 零基础说科研，仅为个人学习用，如有侵权，可以删除
吸烟没什么创意，唯一的创意就是加入了MR和meta分析，作者显示介绍吸烟与多种疾病之间的因果关系扔不明确，
第一步：搜索各大数据库的 MR随机化与吸烟的文章，把文章中吸烟和疾病的关系的数据提取出来，纳入标准， 原始全文文章，介绍了吸烟或者终生吸烟的遗传易感性与一种或多种循环、消化、神经和肌肉骨骼系统疾病， 内分泌、代谢和眼部疾病或肿瘤风险的关联结果，一共纳入了385篇文献
排除标准：基于相同或重叠研究样本的重复出版物，以及仅使用单一或少数尼古丁依赖或吸烟行为或数量的工具变量的研究，剔除了相同或重叠研究样本的重复出版物，以及使用单一或少数尼古丁依赖或吸烟行为或数量的工具变量，的研究，作者这里提取了年份、样本量、关系的比值比
经过排除后适合分析的29篇
第二步： 有一部分数据就是芬兰基因研究（FinnGen）作者通过检索没有检索到资料，他就自己来做，他使用了 R6版本中的数据进行孟德尔随机化分析，其中包括 260 405 名芬兰人，但剔除了性别不明确、非芬兰血统、基因型缺失率超过 5%、或杂合度过高（±4 个标准差）的数据。此外，作者还利用 GWAS meta 分析中公开的汇总统计数据，对骨关节炎、痛风和原发性开角型青光眼进行了从头开始的 MR 分析。第二部分提取到的数据应该是27篇，因为最后供56篇文章。
通过流程图咱们可以知道最后作者得到14篇循环疾病的文章，消化疾病8篇，神经系统疾病5篇，肌肉骨骼系统4篇，2篇内分泌，3篇眼科疾病，21篇关于肿瘤的文章。整个过程处理起来还是挺花时间的，毕竟要一篇篇的读，提取数据。
接下来咱们看看作者提供的数据，附表1是作者自己做的孟德尔随机化的结果，它的结果有两个，一个是刚开始吸烟，还有一个是终生吸烟。作者也是根据这两个结果进行meta分析的
接下来就是作者的两个主表了，表2是开始吸烟的人的疾病分析，表3是终生吸烟的人的疾病分析，作者就是根据这两个表来做meta分析的，下面我把数据提取出来跑一下
下面我把数据提取出来跑一下，数据量挺大的我就提取刚开始吸烟的患者循环疾病这部分举个例子，其他疾病的都是一样的。这里我们要注意一下，循环系统是有很多疾病的，就拿心房颤动这个疾病来说，作者的数据很多很大，他是把GWAS meta-analysis、FinnGen这2个数据库的结果进行相加，再来做meta分析，如果你的数据没有这么大，你把每个数据库先分别做，然后再汇总也是可以的。
bc<-read.csv("E:/r/test/smokemeta.csv",sep=',',header=TRUE)names(bc)
library(forestplot)
bc$`OR(95% CI)` <- sprintf("%.2f (%.2f to %.2f)", bc$OR, bc$LB, bc$UB) #生成可信区间
dt1<-as.matrix(bc[,c(1,2,6)]) #生成绘图区间，选择需要的变量 这里是126
dt1<- rbind(c("outcome","Cases","OR (95%CI)"),dt1) #注意 dt这个数据是矩阵没有列名，还有生成一个列名
接下来就可以绘图
forestplot(labeltext=dt1,graph.pos="right", mean=c(NA,bc$OR), lower=c(NA,bc$LB), upper=c(NA,bc$UB), graphwidth =unit(60,"mm"),#设置图片位置和宽度 boxsize =0.2,line.margin =unit(5,"mm"),#对散点和线条进行设置 lineheight =unit(5,"mm"),#设置图形行距 col=fpColors(box ="grey0",lines = "grey0",summary = "grey0"), colgap = unit(1,"mm"),#图形列间距 zero = 1,#参照值 xticks = c(0,1,2))#X轴的定义标签
meta分析的异质性和P值 用stata最简单
metan or 1b ub 直接就可以出来！
R语言进行进行meta分析咱们就做最基本的部分就行，不必搞得太复杂。咱们先导入数据和R包。很多R包都能做，咱们随便选个meta包就可以了。
library(meta)
bc<-read.csv("E:/r/test/senlintu1.csv",sep=',',header=TRUE)
names(bc)
out<-metabin(event.e=a,
n.e=b,event.c=c,n.c=d,data=bc,sm="OR",studlab=paste(study),method="Inverse")
study代表研究名称；a实验阳性人数，b实验总人数，c对照组阳性人数，d对照组总人数
咱们先来看下函数格式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")
out<-metabin(a,b,c,d,data=bc,sm="OR",studlab= paste(study),method = "Inverse")
summary(out)
结果区，罗列了随机效应和固定效应、 I值、 异质性检查
如果heterogeneity 这里P<0.05 存在异质性，选择随机效应模型
异质性用I2表示的，是89。9%比较大的话，可以使用剔除法剔除
先加ID
bc$id<-1:13
删除第一项，使用亚组函数subset控制，subset=id>2out<-metabin(a,b,c,d,data=bc,sm="OR",studlab= paste(study), method ="Inverse",common=F,subset=id>2)
summary(out)
out<-metabin(a,b,c,d,data=bc,sm="OR",studlab= paste(study), method ="Inverse",common=F)forest(out)
metabias(out,method.bias="Egger") #p值>0.05说明没有偏倚
metabias(out,method.bias="peters")