Perfoming clumping using ieugwasr package
1
1
Entering edit mode
@1b588f00
Last seen 3 months ago
Germany

I am conducting a mendelian randomization study, and as an exposure I have hundreds of triats. The code below is part of a loop that iterates over these triats.


  # Apply Clumping 
  SeqId_data <- ieugwasr::ld_clump_local (
  SeqId_data,
  clump_p = 1,
  clump_r2 = r2,
  clump_kb = 10000,
  plink_bin = genetics.binaRies::get_plink_binary(), 
  bfile= "/path/to/reference/EUR" )

And am getting this output

[1] "1 10001_7"
[1] "Loading dataset with the following SeqId: 10008_43"
PLINK v1.90b6.21 64-bit (19 Oct 2020)          www.cog-genomics.org/plink/1.9/
(C) 2005-2020 Shaun Purcell, Christopher Chang   GNU General Public License v3
Logging to /tmp/RtmphPuzvw/file12b415f4a3e.log.
Options in effect:
  --bfile /home/a/alayash/TestFunctions/Sig1/EUR
  --clump /tmp/RtmphPuzvw/file12b415f4a3e
  --clump-kb 10000
  --clump-p1 1
  --clump-r2 0.099
  --out /tmp/RtmphPuzvw/file12b415f4a3e

132236 MB RAM detected; reserving 66118 MB for main workspace.
8550156 variants loaded from .bim file.
503 people (0 males, 0 females, 503 ambiguous) loaded from .fam.
Ambiguous sex IDs written to /tmp/RtmphPuzvw/file12b415f4a3e.nosex .
Using 1 thread (no multithreaded calculations invoked).
Before main variant filters, 503 founders and 0 nonfounders present.
Calculating allele frequencies... done.
8550156 variants and 503 people pass filters and QC.
Note: No phenotypes present.
Warning: No significant --clump results.  Skipping.
Error in file(file, "rt") : cannot open the connection
In addition: Warning message:
In file(file, "rt") :
  cannot open file '/tmp/RtmphPuzvw/file12b415f4a3e.clumped': No such file or directory

Any idea what might be the problem?

Thank you

LinkageDisequilibrium genom linka GenomicDistributionsData • 656 views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 1 day ago
United States

Plink said

Warning: No significant --clump results.  Skipping.

And then R said it can't open the results file from running Plink. Which won't exist, because Plink already said there were none.

Login before adding your answer.

Traffic: 662 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6