首页 > 其他分享 >本地clump data

本地clump data

时间:2023-09-20 17:48:08浏览次数:34  
标签:plink -- data ieugwasr 本地 clump

在使用TwoSampleMR分析样本的时候,时长遇到clump_data报错,原因是需要联网,其实可以在本地clump数据,使用ieugwasr包里面的ld_clump_data函数
,在安装TwoSampleMR的时候会自动安装ieugwasr包,如果需要最新的包,可以使用:
remotes::install_github("mrcieu/ieugwasr")
该函数说明如下:

关键参数:

dat是需要处理的snp信息,至少要有两列,一列必须叫rsid,内容是rs编号,另一列必须叫pval,内容是snp在GWAS中的p值。
clump_kb ,clump_r2,clump_p按需设置
bfile指的是参考基因组数据,参考链接https://ctg.cncr.nl/software/magma
plink_bin 指的是plink的可执行文件位置,如/usr/bin/plink,可以使用命令sudo apt install plink1.9下载1.9版本的,2.0版本的没有本地clump功能

本地clump关键代码:

expo_clump <- ieugwasr::ld_clump_local(dat = dat,
                         clump_r2 = 0.001,
                         clump_p = 1,
                         clump_kb = 10000,
                         plink_bin = "/usr/bin/plink1.9",
                         bfile = "/home/ref/g1000_eur/g1000_eur"  #该参数第二个g1000_eur指的是g1000_eur.bed等文件,前面哪个g1000_eur是文件夹
)

其实也可以:

plink1.9 --bfile /home/ref/g1000_eur/g1000_eur --clump dat --clump-kb 10000 --clump-p1 1 --clump-r2 0.001 --noweb --out /home/gwas/clump_data

标签:plink,--,data,ieugwasr,本地,clump
From: https://www.cnblogs.com/missed-forest/p/17717931.html

相关文章

  • DataGrip 2023:多功能的数据库管理软件
    DataGrip2023是由JetBrains开发的一款功能强大的数据库管理工具,它旨在提供一个集成的开发环境,方便开发人员管理和操作各种类型的数据库。DataGrip2023支持多种数据库系统,包括MySQL、PostgreSQL、Oracle、SQLServer等,它具有直观的用户界面,使用户能够轻松地连接到数据库服务器,进......
  • IDEA 本地启动项目时候提示:错误:无法加载主类StartMain
    这也许是IDEA的bug.最好的办法就是清理IDEA项目缓存。具体步骤如下:一、检查使用包版本是否正确 二、清除缓存重启 ......
  • Handler dispatch failed; nested exception is java.lang.NoSuchMethodError: com.fa
    报错:Handlerdispatchfailed;nestedexceptionisjava.lang.NoSuchMethodError:com.fasterxml.jackson.databind.ObjectMapper.canSerialize(Ljava/lang/Class;Ljava/util/concurrent/atomic/AtomicReference;)Zjar包冲突,找到对应的jar包删除......
  • 迁移docker的data-root流程
    在部署jumpserver的时候采用了官方推荐的docker部署方式,但是由于前期规划不是很周全,导致docker运行一段时间后,磁盘空间不足,导致jumpserver运行异常。未解决这个问题,考虑将docker的data-root目录进行迁移,迁移到一个空间较为充裕的磁盘目录下Docker的data-root目录是用于存储Docker容......
  • 本地缓存无冕之王Caffeine Cache
    本文已收录至GitHub,推荐阅读......
  • How to get a variable data type in Python 3 All In One
    HowtogetavariabledatatypeinPython3AllInOnetypeofinjstype(var)&isinstance(var,type)#!/usr/bin/envpython3#mixlistlt=[1,2,{'k':'v'},{1,'str'}]dt=dict()for[i,item]inenumerate(lt):......
  • 生产环境,测试,本地
    固定IP地址通常是由宽带运营商提供的,而不是路由器给的。在宽带拨号上网的情况下,路由器通常会使用动态IP地址,即每次连接时都会自动获取一个新的IP地址。而固定IP地址则是由宽带运营商提供给用户使用的,用户需要在路由器设置中手动设置固定IP地址。在某些情况下,如果宽带运营商提供了静......
  • 记录--纯前端如何实现录屏并保存视频到本地
    这里给大家分享我在网上总结出来的一些知识,希望对大家有所帮助作为一个资深的切图仔,我们难免会需要用到把自己写的页面的一些功能通过视频的方式分享给别人。还有一个场景,就是当我们面试的时候,我们需要把我们的屏幕分享给面试官看,那么这些都是怎么实现的呢?那么接下来我们就通......
  • Datainside数据分析,基于大数据分析学生成绩综合评价
    Datainside是一种基于大数据分析的学生成绩综合评价方法,通过对海量学生成绩数据进行深度挖掘和分析,为学生的学习表现提供全面、客观的评价。以下是对Datainside数据分析学生成绩综合评价的详细描述:1.数据收集:Datainside首先需要收集大量学生成绩数据,包括不同科目、不同年级、不......
  • 16G内存+CPU本地部署ChatGLM2/Baichuan2推理(Windows/Mac/Linux)
    概述本文使用chatglm.cpp对中文大语言模型(LLM)进行量化与推理,支持ChatGLM2-6B、Baichuan2-13B-Chat等模型在CPU环境16G内存的个人电脑上部署,实现类似ChatGPT的聊天功能。支持的操作系统包括Windows、MacOS、Linux等。其中,量化过程需要临时使用一台内存较大的服务器。4bit量化后......