为什么需要做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. 【Android 逆向】Android 系统文件分析 ( /system/ 系统命令和系统应用数据目录 | /system/app/ 系统应用目录 | sys Linux 系统内核文件目录 )
  2. 实战ELK(5) Logstash 入门
  3. 微信小程序页面间传递json数据
  4. 软件构造学习笔记-第四周
  5. java 事物_java事务
  6. CentOS英文版下安装中文语言包
  7. ToPILImage
  8. API接口应该如何设计?
  9. [Ubuntu Error]鼠标指针闪烁的解决方法
  10. 架构猿怒了,放大招:详解DoDAF
  11. Java基于spingboot+vue智慧小区物业管理系统源码 智慧小区管理系统源码
  12. docker 停止、启动、删除镜像指令操作总结
  13. 体育教学硕士毕业论文题目
  14. Matconvnet测试
  15. 对模型进行DP处理(DP-SGD+DP-Logits)
  16. 系统测试之功能测试:测试用例的设计步骤——从登陆开始说起
  17. 1132: 数字字符统计 C语言
  18. html中取颜色快捷键,ps中填充颜色的快捷键是什么?
  19. 大数据研究的若干科学问题——徐宗本
  20. 虚拟环境中cuda版本与系统cuda版本不同有什么影响

热门文章

  1. 2020年mac上最值得使用的五大录屏软件
  2. maven配置tomcat版本
  3. SQL Server(2019)数据库----数据库和表的建立
  4. c语言程序设计必备单词32个,c语言编程必背单词-20210324071350.docx-原创力文档
  5. 群晖6.1安装php3.6_黑群晖 DSM6.1安装简明教程(附下载)
  6. Java常用算法手册-01算法概述
  7. python3GUI--刷屏器(附源码)
  8. 零起点学习Linux系列培训视频-寒冰作品
  9. VS2010使用教程(使用VS2010编写C语言程序)
  10. 13G311-1 混凝土结构加固构造 免费下载