Question: locuszoom on a HPC - LD could not be computed for SNP
0
gravatar for tacrolimus
9 weeks ago by
tacrolimus80
tacrolimus80 wrote:

Dear Biostars,

I am running locuszoom via the command line on a HPC to investigate GWAS hits. I have created a "metal" file from my association results outputted from a GWAS performed in PLINK1.9.

When running a single region I get the error:

"PLINK did not complete successfully. Please check the log file to see the directory with the log file. Warning: LD could not be computed for SNP chr1:234590242. This SNP does not exist in the genotype files for computing LD from 1000G_Nov2014/hg38/EUR.. Warning: recombination rate table 'reocomb_rate' not found in database, skipping recomb rate lookups in region. No usable LD information"

My script is this: locuszoom --metal [input dile] --build hg38 --pvalcol P-Value --no-cleanup --markercol MarkerName --pop EUR --source 1000G_Nov2014 --build hg38 -chr 1 --start 234588000 --end 234590742 --prefix x

The actual input files I can't show as I work in an airlock environment but they are similar to the example from the locus zoom website and are tab separated. Example below:

MarkerName P-value chr4:401141 9.4e-390

The log files from plink give: Error: No valid variants specified by --ld-snp/ --ld-snps/ --ld-snp-list .

I get a plot but it is very sparse compared to what I am expecting.

Has anyone else encountered this problem before?

Many thanks

bash gwas locuszoom hpc • 113 views
ADD COMMENTlink written 9 weeks ago by tacrolimus80
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1026 users visited in the last hour
_