search

KonradZych / phenotypes2genotypes

Package for creating genetic map using gene expression data.
6 stars 5 forks source link

Missing values #42

Open DannyArends opened 11 years ago

DannyArends commented 11 years ago

normalmixEM from the mitools package cannot handle miture modeling on vectors with a LOT of NA's

eg.: When fitting 2 normals 50% (+1) should be present When fitting 3 normals 66% (+1) should be present

However the currect strategy is t allow massive NAs in pheno2geno... to solve this we can simply find all the na's, and remove them before modeling. (assuming the data is in V)

idx <- is.na(v)
v <- v[idx]
#Call normalmixEM
Fix the resulting output using the idx vector
KonradZych commented 11 years ago

Using NaN insted of NAs and removing code duplication.