clump的时候，报错：The query to MR-Base exceeded 300 seconds ，尝试做本地clump

参考了一下前人的经验，这里是参考的原文的链接

TwoSampleMR:local clump(MR-Base exceeded 300 seconds) 包括Windows和Linux R解决办法_Bioinfoer的博客-CSDN博客

下面是我自己跑的代码，在原文的基础上加了一些我的修改和解读

先安装包

devtools::install_github("mrcieu/ieugwasr",force = TRUE)

更新所有的包

library(ieugwasr)

devtools::install_github("explodecomputer/plinkbinr")

library(plinkbinr)

get_plink_exe()

下载bfile，我的GWAS是欧洲人种所以用EUR*，存到本地G盘的EUR_ref/EUR文件夹里

#http://fileserve.mrcieu.ac.uk/ld/1kg.v3.tgz

#解压该文件夹

b <- ld_clump(

dplyr::tibble(rsid=a$rsid, pval=a$p, id=a$id),

#get_plink_exe()

plink_bin = "C:/Users/Thinkpad/Documents/R/win-library/4.0/plinkbinr/bin/plink_Windows.exe",

#欧洲人群参考基因组位置

bfile = "D:/EUR_ref/EUR"

)

上面这一段代码在运行的时候又遇到了下面的报错

按照路径打开了该文件，应该就是clump没有成功，导致生成临时文件失败，在clump成功的文件里面找相应的文件也是没有的，因为这个只是一个临时生成的中间文件，只不过在这一步临时文件生成没有成功。

后面看到上文作者改了源代码

这个是源代码的网址

GitHub - MRCIEU/ieugwasr: R interface to the IEU GWAS database API

把源代码的zip下载到G盘，下载下来之后本地解压命名为ieugwasr

作者改的就是ld_matrix这个文件

用R打开，在相应的位置加上 " --threads 1 "，

然后保存，然后再本地重新安装ieugwasr

setwd("G:/ieugwasr")

devtools::install_local() 安装完成之后

library(ieugwasr)

library(plinkbinr)

get_plink_exe()

trait1.exposure_data_clumped <- ld_clump(

dplyr::tibble(rsid=trait1.exposure_data_1$SNP,

pval=trait1.exposure_data_1$pval.exposure,

id=trait1.exposure_data_1$exposure,

beta_col = trait1.exposure_data_1$beta.exposure,

se_col = trait1.exposure_data_1$se.exposure,

eaf_col = trait1.exposure_data_1$eaf.exposure),

#get_plink_exe()

plink_bin = "C:/Users/Thinkpad/Documents/R/win-library/4.0/plinkbinr/bin/plink_Windows.exe",

#欧洲人群参考基因组位置

bfile = "G:/EUR_ref/EUR"

)

运行之后，还是出现了上述的问题。

后来跟在线数据库的clump比对了一下

在线clump

编辑切换为居中

本地clump

本地clump的数据，自己在 dplyr::tibble列哪些数据才展现哪些数据，比在线的clump的数据要少一些。

然后出现前面：The query to MR-Base exceeded 300 seconds的问题，主要是我把很多exposure的snp放在一个txt文档里面，想一起做clump，后来觉得把exposure 逐个拆开，用for循环来逐个文件做clump。

clump的时候，报错：The query to MR-Base exceeded 300 seconds ，尝试做本地clump相关推荐

jpa Specification fetch查询报错，query specified join fetching, but the owner of the fetched association
jpa Specification fetch查询报错,query specified join fetching, but the owner of the fetched association ...
php数据库报错1046怎么解决,DZ报错MySQL Query Error Errno:1046错误的解决办法
这篇文章主要为大家详细介绍了DZ报错MySQL Query Error Errno:1046错误的解决办法,具有一定的参考价值,感兴趣的小伙伴们可以参考一下,有需要的朋友可以收藏方便以后借鉴. DZ报 ...
python编译器报错：“RecursionError: maximum recursion depth exceeded in comparison”解决方案
python编译器报错:"RecursionError: maximum recursion depth exceeded in comparison"解决方案在使用递归迭代语句 ...
已解决-MacBook重装系统报错：准备安装时发生错误，请尝试重新运行此程序
我是月球挖掘机,今天重新安装了下macbook air的系统,用联网恢复的方式安装的时候,结果总是报错:准备安装时发生错误,请尝试重新运行此程序.网上找了方法用了还是不行,最后发现自己的时间虽然按照设 ...
蜻蜓短视频系统-文件上传之本地上传报错问题解决-蜻蜓Q蜻蜓s蜻蜓系统上传采用本地上传报错Class ‘League\Flysystem\Adapter\Local‘ not found
蜻蜓短视频系统-文件上传之本地上传报错问题解决-蜻蜓Q蜻蜓s蜻蜓系统上传采用本地上传报错蜻蜓系统有本地上传方式和云储存上传方式,云方式比较常用是不会有错的,本地上传方式可能会遇到这样的报错. 报错1 ...
vue3项目使用了element-plus中的el-table表格组件，缩放浏览器网页时，报错如下： ResizeObserver loop limit exceeded at eval (w
错误描述: vue3项目使用了element-plus中的el-table表格组件,缩放浏览器网页时,报错如下: ResizeObserver loop limit exceeded at eval ...
mysql bingip,MySQL报错Ignoring query to other database的真正原因
大家都知道报这个错的原因是忘了用户名前加-u这个参数造成的,比如是这样登录了MySQL # mysql -root -p (注意看,-root前面没有加u) 登录后,不管你执行什么命令,都会提示Ign ...
Mysql删除重复数据并解决You can't specify target table 'xx' for update in FROM clause 报错与 query interrupted报错
哇,这个问题在纠结了N个小时总算解决了! 解决方案: -- 在命令行方式中输入 delete from 表名 where id in (select id from (select id from 表 ...
kafka-connect-jdbc 增加source报错：query may not be combined with whole-table copying settings
kafka版本:0.10.2.0 kafka-connect-jdbc插件:10.3.3 增加一个MySQL表作为源监听数据变化. 配置: {"name": "mysql ...

clump的时候，报错：The query to MR-Base exceeded 300 seconds ，尝试做本地clump

clump的时候，报错：The query to MR-Base exceeded 300 seconds ，尝试做本地clump相关推荐

最新文章

热门文章