首页 技术 正文
技术 2022年11月23日
0 收藏 514 点赞 3,871 浏览 1166 个字

目前GWAS方向发了很多文献,但是并没有一个很完善的R包对这些文献的数据进行汇总。

接下来推荐的这个是最新发表的GWAS数据汇总R包​。看了一下功能齐全,但是数据不是收录的很齐全​。
下面具体讲一下。​

在R上安装

install.packages("remotes")

remotes::install_github("ramiromagno/gwasrapidd")

查询是否有发表关于自身免疫性疾病的文章

library(gwasrapidd)

my_studies <- get_studies(efo_trait = 'autoimmune disease')

查询发表自身免疫性疾病文章的数量

n(my_studies)

获取发表自身免疫性疾病的文章ID

my_studies@studies$study_id

查询发表自身免疫性疾病的文章标题

my_studies@publications$title
查询、下载GWAS目录数据的R包(gwasrapidd)

查看发表自身免疫性疾病文章在pubmed的界面信息

open_in_pubmed(my_studies@publications$pubmed_id)

获取发表身高文献的SNP关联信息

my_associations <- get_associations(study_id = my_studies@studies$study_id)

查询P值小于1e-6的位点

dplyr::filter(my_associations@associations, pvalue < 1e-6) %>% # Filter by p-value
tidyr::drop_na(pvalue) %>%
dplyr::pull(association_id) -> association_ids # Extract column association_id

提取显著信号位点信息

my_associations2 <- my_associations[association_ids]

显示显著信号位点的个数

n(my_associations2)

显示显著信号位点的RS ID, risk allele, 频率

my_associations2@risk_alleles[c('variant_id', 'risk_allele', 'risk_frequency')] %>%

print(n = Inf)

获取含有rs12752552位点的文献

s2 <- get_studies(variant_id = 'rs12752552')

测试了一下,优点是减去了很多信息检索的工作,缺点是这个包刚开发不久,信息还不够齐全

更多详细信息请看gwasrapidd

参考文献:Magno R, Maia A T. gwasrapidd: an R package to query, download and wrangle GWAS Catalog data[J]. bioRxiv, 2019: 643940.

微信扫一扫

支付宝扫一扫

本文网址:https://www.zhankr.net/141788.html

相关推荐
python开发_常用的python模块及安装方法
adodb:我们领导推荐的数据库连接组件bsddb3:BerkeleyDB的连接组件Cheetah-1.0:我比较喜欢这个版本的cheeta…
日期:2022-11-24 点赞:875 阅读:5,075
Educational Codeforces Round 11 C. Hard Process 二分
C. Hard Process题目连接:http://www.codeforces.com/contest/660/problem/CDes…
日期:2022-11-24 点赞:806 阅读:3,505
下载Ubuntn 17.04 内核源代码
zengkefu@server1:/usr/src$ uname -aLinux server1 4.10.0-19-generic #21…
日期:2022-11-24 点赞:565 阅读:4,317
可用Active Desktop Calendar V7.86 注册码序列号
可用Active Desktop Calendar V7.86 注册码序列号Name: www.greendown.cn Code: &nb…
日期:2022-11-24 点赞:730 阅读:4,311
Android调用系统相机、自定义相机、处理大图片
Android调用系统相机和自定义相机实例本博文主要是介绍了android上使用相机进行拍照并显示的两种方式,并且由于涉及到要把拍到的照片显…
日期:2022-11-24 点赞:512 阅读:4,912
Struts的使用
一、Struts2的获取  Struts的官方网站为:http://struts.apache.org/  下载完Struts2的jar包,…
日期:2022-11-24 点赞:671 阅读:3,099
发表评论
暂无评论

还没有评论呢,快来抢沙发~

助力内容变现

将您的收入提升到一个新的水平

点击联系客服

在线时间:8:00-16:00

客服电话

400-888-8888

客服邮箱

ceotheme@ceo.com

扫描二维码

关注微信公众号

扫描二维码

手机访问本站