读取GWAS结果的vcf格式文件
最近米老鼠收到好多朋友关于openGWAS数据库的vcf文件处理问题，在这里我打算写一期推送帮助大家系统解决一下这个问题。这里我们主要使用的是“vcfR”这个包，它可以快速读取vcf格式的文件。
我们以糖化血红蛋白的GWAS为例（https://gwas.mrcieu.ac.uk/datasets/ieu-b-103/）进行讲解，首先我们打开网页并下载该vcf文件（黄色标记）。
install.packages("vcfR") #安装R包arData <-vcfR::read.vcfR("./ieu-b-103.vcf.gz") #直接读取vcf文件str(arData) #查看vcf数据的结构从上图中我们不难看出，这里的vcf文件主要由三部分内容构成（meta，fix和gt），其中meta是vcf的注释信息。
阅读原文 最新评论
推荐文章
作者最新文章
你可能感兴趣的文章
Copyright Disclaimer: The copyright of contents (including texts, images, videos and audios) posted above belong to the User who shared or the third-party website which the User shared from. If you found your copyright have been infringed, please send a DMCA takedown notice to [email protected]. For more detail of the source, please click on the button "Read Original Post" below. For other communications, please send to [email protected].
版权声明：以上内容为用户推荐收藏至CareerEngine平台，其内容（含文字、图片、视频、音频等）及知识版权均属用户或用户转发自的第三方网站，如涉嫌侵权，请通知[email protected]进行信息删除。如需查看信息来源，请点击“查看原文”。如需洽谈其它事宜，请联系[email protected]。
版权声明：以上内容为用户推荐收藏至CareerEngine平台，其内容（含文字、图片、视频、音频等）及知识版权均属用户或用户转发自的第三方网站，如涉嫌侵权，请通知[email protected]进行信息删除。如需查看信息来源，请点击“查看原文”。如需洽谈其它事宜，请联系[email protected]。