为什么需要做meta分析

群体分层是GWAS研究中一个比较常见的假阳性来源.

也就是说,如果数据存在群体分层,却不加以控制,那么很容易得到一堆假阳性位点。

当群体出现分层时,常规手段就是将分层的群体独立分析,最后再做meta分析。

1.如何判断群体是否分层

随后画出PC1和PC2在不同群体的散点图,观察群体之间是否明显分开,如果明显分开,说明群体分层了,需要独立做关联分析,最后再做meta分析

2.如何做meta分析

这里推荐metal软件做meta分析,理由是简单、易上手。

2.1 下载metal

进入下载链接

metal 提供了三个版本的,分别是Linux,macOS, Windows系统;请自行选择。

这里提供Linux系统的命令:

wget http://csg.sph.umich.edu/abecasis/metal/download/Linux-metal.tar.gz

2.2 解压metal

解压用到的命令如下:

tar -zxvf Linux-metal.tar.gz

2.3 meta分析前的数据准备

假定需要进行meta分析的文档分别为DGI_three_regions.txt 和 magic_SARDINIA.tbl

DGI_three_regions.txt 的内容如下:

magic_SARDINIA.tbl 的内容如下:

那么在meta分析前需要准备一个metal.txt文档,metal.txt文档的内容如下:

解释一下,这个txt文档是什么意思。

这一部分指的是MARKER对应的是DGI_three_regions.txt文档的SNP列名;

WEIGHT对应的是DGI_three_regions.txt文档的SNP列名;

其他的以此类推;

第二个文件的准备方法也是一样的。

2.4 meta分析

很简单的一个命令行就搞定了

metal metal.txt

2.5 结果解读

meta分析后会生成两个文件,分别是 METAANALYSIS1.TBL 和 http://METAANALYSIS1.TBL.info

METAANALYSIS1.TBL 是meta分析的结果文档;

内容如下:

P-value 即为meta后的关联分析P值;

METAANALYSIS1.TBL.info 是meta分析的说明文档,比如 Marker 指的是什么。

其内容如下:

如有问题请移步微信公众号“bio生物信息”或“orange_milk_sugar”向我提问。

metal分析是什么意思_GWAS数据如何做meta分析?相关推荐

  1. 森林图怎么分析_资料|用Stata怎么做Meta分析?看完你就懂啦

    传统的Meta分析软件除了Stata,还有Revman.R.SAS.Meta-Analyst,对于初学者来说,小助理建议还是学Revman,因为菜单操作简单,上手比较快,但是如果你对统计软件或者计算机 ...

  2. 熟练掌握R语言的Meta分析全流程和不确定性分析,并结合机器学习等方法讲解Meta分析在文献大数据的延伸应用

    Meta分析是针对某一科研问题,根据明确的搜索策略.选择筛选文献标准.采用严格的评价方法,对来源不同的研究成果进行收集.合并及定量统计分析的方法,最早出现于"循证医学",现已广泛应 ...

  3. 如何用python爬取数据_如何使用python爬取知乎数据并做简单分析

    原标题:如何使用python爬取知乎数据并做简单分析 一.使用的技术栈: 爬虫:python27 +requests+json+bs4+time 分析工具: ELK套件 开发工具:pycharm 数据 ...

  4. python 抓取解析接口数据_[干货]用python抓取摩拜单车API数据并做可视化分析(源码)...

    原标题:[干货]用python抓取摩拜单车API数据并做可视化分析(源码) 在APP中能看到很多单车,但走到那里的时候,才发现车并不在那里.有些车不知道藏到了哪里:有些车或许是在高楼的后面,由于有GP ...

  5. 做Meta分析的10个原则,建议收藏。

    近年来文献数量显著增加,大量不同主题的Meta分析数量也显著增长.在生物医学领域,Meta分析是多种发表类型中拥有较高被引数的类型.针对Meta分析已经有不少书籍或指南,本文是一篇指导Meta分析的简 ...

  6. 森林图怎么分析_森林图(forest plot)怎么看|meta分析

    Hello, 大家好,我是钟多少玲这是本公众号的第15篇更文第一次拖更,感谢大家后台(当面)催稿先碎碎念几句--十月是繁忙又收获的季节,很幸运期间收到4篇文章accepted的消息,人生中第一.第二篇 ...

  7. 爬取链家北京租房数据并做简单分析

    在一个来北京不久的学生眼中,北京是一个神秘又充满魅力的大城市.它无比美好,但又无时无刻不再觊觎这你薄弱的钱包. 租房是很多人都离不开的硬性需求,这里就对从链家爬取的北京地区房屋出租数据进行一个简单分析 ...

  8. 手工填报的数据能做可视化分析吗?

    手工填报的数据能不能做可视化分析?该怎么做?效果如何? 我们公司既有手工填报的数据,也有存放在其他不同数据系统中的,能不能都做成可视化分析效果? 如果你也有以上疑惑,请继续往下看! 手工填报.数据源不 ...

  9. r语言数据变量分段_R数据分析:用R语言做meta分析

    这里以我的一篇meta分析为例,详细描述meta分析的一般步骤,该例子实现的是效应量β的合并 R包:metafor或meta包,第一个例子以metafor包为例. 1.准备数据集 2.异质性检验 in ...

  10. aspen怎么做灵敏度分析_不会敲代码怎么做智能分析?用奥威商业智能BI

    不会敲代码,怎么学也不会,那怎么做大数据智能分析?说来也简单,用商业智能BI软件就行.拖拉拽就能做报表,点击即可智能分析,即使是没有IT基础,完全不会敲代码的人也能轻松抽取数据做智能分析. 商业智能B ...

最新文章

  1. 八十二、什么是SpringCloud项目,我把它讲清楚了
  2. Clickhouse安装及使用/clickhouse-mysql安装
  3. WebClient UI忽略所有增强的开关
  4. 软件工程——快速掌握面向对象开发方法
  5. git 上的 python package 的安装
  6. 27.MySQL Status
  7. lightning接口_苹果、安卓充电接口有望统一,欧盟重压下,苹果将不得不妥协
  8. flex blazeds java spring_flex+blazeds+java+spring后台消息推送
  9. getchar()用法
  10. 游戏与幼儿教育如何相辅相成
  11. 为什么ppt图形卡配置不正确_PPT常见问题解决方法,PPT检测到图形卡可能配置不正确怎么办?...
  12. Swift编程二十四(不透明类型)
  13. 【弄nèng - SpringCloud】应用篇 —— 整合Zipkin【Finchley版】
  14. Linux 开机提示kernel panic...解决方法
  15. 一些特殊字符,由于编码问题显示不出来
  16. c# 计算圆锥的体积_用C#如何编写程序计算球,圆柱和圆锥的表面积和体积?
  17. 动态视频目标检测和跟踪技术
  18. 速学堂(java)第十章编程题答案(自写)
  19. 计算机考试题画图板,9月全国计算机一级Photoshop考试试题操作题
  20. 建设银行的短信服务原来是这样取消的。

热门文章

  1. 免费的录屏软件Capture
  2. 如何将苹方字体写入html,在网页上使用苹方字体
  3. Dictionary 索引超出数组界限
  4. 小白学 Python 爬虫(18):Requests 进阶操作
  5. ghost还原固态硬盘_解决方法:SSD可以使用Ghost软件吗?最后说清楚了
  6. json格式转换工具包
  7. 分形、混沌理论、集异璧之大成
  8. Android面试心得,已拿到字节跳动,京东,android应用开发基础学堂在线答案
  9. grasshopper for rhino 6下载_福利篇 | Rhino常用插件功能汇总
  10. qt添加控制台语句_向控制台语句添加样式