代码:
neuro <- system.file('select_snp.csv',package = 'TwoSampleMR')
neuro_exp_dat<- read_exposure_data(
filename=neuro,
clump=T,
sep = ",",
snp_col = "SNP",
beta_col = "beta",
se_col = "se",
eaf_col = "EAF_UKB",
effect_allele_col = "A1",
other_allele_col = "A2",
pval_col = "P"
)
报错:
No phenotype name specified, defaulting to 'exposure'.
API: public: http://gwas-api.mrcieu.ac.uk/
Please look at vignettes for options on running this locally if you need to run many instances of this command.
Clumping zEFqI0, 7759 variants, using EUR population reference
Server code: 503; Server is possibly experiencing traffic, trying again...
Server code: 503; Server is possibly experiencing traffic, trying again...
Server code: 503; Server is possibly experiencing traffic, trying again...
Server code: 503; Server is possibly experiencing traffic, trying again...
Server code: 503; Server is possibly experiencing traffic, trying again...
Server code: 503; Server is possibly experiencing traffic, trying again...
Server error: 503
Failed to retrieve results from server. See error status message in the returned object and contact the developers if the problem persists.
No encoding supplied: defaulting to UTF-8.
Error in mtfrm.default(list(url = "http://gwas-api.mrcieu.ac.uk/ld/clump", :
cannot mtfrm
引用 皆我百晓生 小程序回复内容作答:
在进行clump操作时,报错信息中提到服务器返回503错误,这意味着服务器无法处理请求,可能是由于服务器当前负载过重或正在维护。由于这是一个网络请求,你可以尝试等待一段时间后重新运行代码,看看问题是否解决。
另外,clump操作通常需要提供两个参数:rsID和p值。从你提供的代码来看,只提供了rsID(使用snp_col参数指定)和p值(使用pval_col参数指定),但没有指定有效的两个样本群体。你可能需要提供适当的参数来指定参考样本和受试样本的人口。你可以查阅相关文档或向开发者咨询,以获取更多关于如何正确使用该函数的信息。