reworkhow / JWAS.jl

Julia for Whole-genome Analysis Software
http://QTL.rocks
GNU General Public License v2.0
96 stars 44 forks source link

fix two bug due to datatype (Float32 or Float 64) and the variance module #150

Closed zhaotianjing closed 6 months ago

zhaotianjing commented 6 months ago

fixed two bugs for single-step:

  1. use geno.genetic_variance.val instead of geno.genetic_variance because now it is the Variance struct.
  2. in impute_genotypes(), make sure the data type (Float32 or Float64) of genotypes is unchanged after imputation. Otherwise the impute_genotypes() will only create Float64 genotype matrix, and this will cause an error for BLAS.axpy!(oldAlpha-α[j],x,yCorr) when x is Float64 but yCorr is Float32.