reworkhow / JWAS.jl

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

Genotype file problem #154

Closed feasiblee closed 2 weeks ago

feasiblee commented 1 month ago

Hello, what should I do if the genotype file has NA values ​​and then read it?

zhaotianjing commented 1 month ago

Hi, in get_genotypes() function, there is an argument of missing_value, and by default, it is missing_value=9.0. I suggest you replace the NA with a 9.0 before your analysis.

feasiblee commented 1 month ago

Problem solved, thank you

feasiblee commented 1 month ago

Hello, I am doing out=runMCMC(model,phenotypes), and ERROR appears: Some levels in CW:id for individuals of interest are not found in training individuals (IDs with non-missing records).You may delete those rows or replace those values ​​with missing. For this question, my phenotypic data reading is like this: phenotypes = CSV.read(phenofile,DataFrame,delim = ',',header=true,missingstring=["NA"]).