MRCIEU / gwasglue

Linking GWAS data to analytical tools in R
Other
112 stars 39 forks source link

problem in gwasvcf::query_gwas #34

Open zh-zhang1984 opened 2 years ago

zh-zhang1984 commented 2 years ago

Dear authors, I have the following error reporting, are there any hints for this?

# Extract the instruments for HDL cholesterol
gwasvcf::set_bcftools(
  NULL
)
HDL_exp <- gwasvcf::query_gwas(
  "/Users/zhangzhongheng/Documents/2022/GWAS_sepsis/DATA/ieu-b-109.vcf.gz", 
  pval = 5e-10)
HDL_exp <- gwasglue::gwasvcf_to_TwoSampleMR(
  HDL_exp, type="exposure")
#sepsis outcome 
sepsis_out <- gwasvcf::query_gwas(
  "/Users/zhangzhongheng/Documents/2022/GWAS_sepsis/DATA/ieu-b-4982.vcf.gz", 
  rsid = HDL_exp$SNP
    )
'tools_bcftools' option is not set, using native read which may be substantially slower. See 'set_bcftools' for information.
Note, this is much slower than searching by chromosome/position (e.g. see query_chrompos_file)
starting prefilter
prefiltering 12243372 records
Warning: TRE pattern compilation error 'Out of memory'Error in h(simpleError(msg, call)) : 
  error in evaluating the argument 'i' in selecting a method for function 'extractROWS': Filter 'fil' failed: invalid regular expression 'rs1191478940|rs1269057460|rs386627425|rs2298214|rs1936360|rs61703480|rs2749141|rs113177974|rs1291065947|rs2144734|rs1457642705|rs664755|rs10737472|rs1558013721|rs2298632|rs649458|rs1198443|rs1212034|rs10753556|rs11340914|rs72654647|rs1370916257|rs147266247|rs961976684|rs60659886|rs1203827760|rs34232231|rs34789477|rs34728853|rs113520491|rs34175245|rs12753981|rs148268445|rs12759594|rs71640337|rs35989952|rs71007882|rs34217609|rs60551503|rs71640338|rs12746810|rs12754468|rs35917448|rs35161452|rs12738345|rs1254174435|rs35738294|rs71636769|rs71636770|rs114165349|rs34585496|rs4579782|rs57217461|rs12743862|rs12744003|rs6666121|rs77408691|rs71514277|rs71007886|rs34423276|rs57532083|rs12737407|rs12737912|rs190077222|rs12729444|rs12726081|rs12730386|rs79568175|rs71581091|rs12726287|rs12756154|rs35407375|rs78533482|rs566260414|rs34517168|rs1369054892|rs800