rgcgithub / regenie

regenie is a C++ program for whole genome regression modelling of large genome-wide association studies.
https://rgcgithub.github.io/regenie
Other
183 stars 53 forks source link

Genotype going to -nan or nan #483

Open jerome-f opened 9 months ago

jerome-f commented 9 months ago

Hi,

I am facing an issue with the burden test where the genotype goes to -nan or nan:

     -reading in genotypes, computing gene-based tests and building masks...ERROR: Throw location unknown (consider using BOOST_THROW_EXCEPTION)
Dynamic exception type: boost::wrapexcept<std::domain_error>
std::exception::what: Error in function boost::math::cdf(const complement(normal_distribution<double>&), double): Random variate x is -nan, but must be finite!

The genotypes are from UKBB exome sequencing and I am using the latest v3.4.

joellembatchou commented 8 months ago

Hi,

Can you please retry with the newest v3.4.1?

jerome-f commented 8 months ago

@joellembatchou will try it out and report thanks for looking in to this.