ERROR with METAL · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

ERROR with METAL

Good day,

I am busy performing a meta-analysis with my own GWAS results with summary statistics from a paper that has done a similar study to mine with different cohorts. The paper has their MarkerNames as rsIDs while my own summary statistics are in chr:bp:a1:a2 format (eg 10:87086125:A:G ). I converted over the rsIDs to match my own study MarkerNames. However, when I run my METAL analysis I get the following error

###########################################################################
## Running second pass analysis to evaluate heterogeneity...
## Processing file '/home/dominicz/GWAS_SCA/metal/SUDSMETAL_N_with_header'
## ERROR: Input file has changed since analysis started
## WARNING: Heterogeneity analysis failed...

Has anyone run into this issue or know how to fix it?

Thanks!

gwas metal meta-analysis error

• 2,370 views • link

updated 5.4 years ago by bjyxszd •

• written 5.5 years ago by d.s.zimmerman •

I don't know METAL.

https://github.com/statgen/METAL/blob/master/metal/Main.cpp#L1476 error occurs when the metadata 'header' is different from a 'history' file.

• 0 views • link

written 5.5 years ago by Pierre Lindenbaum

Hi there,

I see, I am not quite sure what a history file is though? I don't see that included in the quick start guide. Thank you for your reply though!

• 0 views • link

written 5.5 years ago by d.s.zimmerman •

Hi,

I got the same error and came across your pose. Did you find out why?

Thanks!

• 0 views • link

written 5.4 years ago by bjyxszd •

0 answers

No answers yet.

Log in to answer this question.

More posts like this

METAL ERROR: The command you issued could not be processed ...
written 2.2 years ago by JJ •

Hello, the meta-analysis seems to work in metal but I keep seeing this error which the log file does not indicate what the error is …
problems with Meta-analysis in METAL
written 3.3 years ago by ek •

Hi everyone, I'm trying to run an inverse variance weighted meta-analysis in METAL. I used the below script: #Metal script for running an inverse SE …
Help in replicating LDSC heritability estimates
written 3.3 years ago by Jerry •

Hi, I am trying to replicate the heritability estimates based on the [insomnia GWAS][1] summary statistics using LDSC. However, I have encountered a problem as …
Error in using the MetaSubtract package
written 3.8 years ago by Khaleesi95 •

Hi guys, I am working with the MetaSubtract package ([ https://rdrr.io/cran/MetaSubtract/][1]) in R. This package enables to subtract from the summary statistics of a meta-analysis …
(LDSC) - Munge Error - ValueError: could not convert string to float: OR
written 4.7 years ago by d.s.zimmerman •

Hi everyone, I am trying to munge some data for later use in ldsc and I come into this error: /home/expcard/Projects/GWAS_SCA/GWAS_NTR/LDSC/ldsc/munge_sumstats.py \ --sumstats /home/dominicz/LDSC/SCAMILIFELINESforMETALnoSNPFinalLDSCmunge.txt --N …
LDSC -Genetic correlation, some results all have NA valuies
written 5.3 years ago by d.s.zimmerman •

Hi there, I am performing a genetic correlation analysis and receive NA results for some of the traits I'm testing but works for others. This …
How to get RSIDs of indels from GWAS summary statistics
written 5.5 years ago by Fede_Santos95 •

Hi, I have a GWAS summary statistics file with SNPs and indels positions, and I'd like to get their rsIDs. My .txt file has the …
LDSC - IOError: Not a gzipped file
written 5.5 years ago by d.s.zimmerman •

Hi All, I am trying to perform genetic correlation with a few summary stastics using ldsc (https://github.com/bulik/ldsc) but I keep receiving the same error. I …
DESeq2: "model matrix is not full rank"
written 5.7 years ago by biogrin •

Hi everyone, My lab has differentiated iPSCs cell lines and I need to do a bioinformatic analysis to try to understand how close they are …
what is fair comparison between results generated from one again other pipeline ??
written 10.8 years ago by unique379

Dear all, I have interesting discussion. I have developed my own pipeline which analyze NGS seq to functional enrichment. So, I have tested my pipeline …