Ryan Welch
Ryan Welch
I think there's some history here on why that behavior wasn't really an issue in the past. Back when METAL was written, association results usually did not have REF or...
Indeed sounds like a bug. Will take a look at patching in your change shortly. Does your organization accept hotfixes, or do they require that it be a new "release"...
Any chance we could use your data for debugging purposes?
We could maybe setup an SFTP or a Dropbox file request (as long as the data is less than a few TB.) Shoot me an email at [email protected] and we'll...
For regular tab-delimited input, most of the time is spent iterating over the file, and translating rsID to chr/pos. The SQLite lookups for each variant are expensive in very large...
Seems like it is trying to compute a z-statistic using ninv() from a p-value of 0, which isn't possible, so it throws an error. https://github.com/statgen/METAL/blob/ab5d19155aa72cef5ca5d2876aea0dd6b9fcba0d/libsrc/MathStats.cpp#L234 Is there an entry in...
Could be the `%.2g` can't handle a very small value beyond double precision, so it writes out 0. What's the smallest p-value seen in any of your studies? Are any...
METAL looks to be compiled with a constant value of 1E-1000 as the minimum p-value accepted when using long double precision: https://github.com/statgen/METAL/blob/ab5d19155aa72cef5ca5d2876aea0dd6b9fcba0d/libsrc/MathStats.cpp#L229-L234 However on most modern systems, a long double...
Are you using the newly produced executable from the recompiled build? Maybe your script is still using `metal` on your $PATH. You could modify the error message to double check...