SMR(Summary-data-based Mendelian Randomization)使用记录

本文介绍了如何将GWASsummarydatatxt格式转换为besd文件,包括使用smr工具处理多个besd文件的合并,并详细说明了如何通过beqtl功能获取特定基因的相关数据,生成新的besd文件。

①GWAS summary data普通txt格式转换为besd文件。(SMR | Yang Lab

smr --eqtl-flist my.flist --make-besd --out mybesd

my.flist内容为

 

 

②合并多个besd文件(尤其已有按照染色体分开的besd)

./smr_v1.3.1_linux_x86_64_static --besd-flist my_file.list --make-besd --out ./LBC_BSGS_meta/my_sparse  

my_file.list内容为

path1/my_besd1 #my_besd1代表单个的besd文件,path为它的路径

path2/my_besd2

path3/my_besd3

③#eqtl 获取某些基因的相关数据
./smr_v1.3.1_linux_x86_64_static --beqtl-summary ./cis-eQTL/cis-eQTLs-full_eQTLGen_AF_incl_nr_formatted_20191212.new.txt_besd-dense --genes gene.list --cis-wind 1000 --query 5

评论 3
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值