Entering edit mode
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