fix bug when sorting genotypes with NA values

TiagoOlivoto · Jun 9, 2021 · c62a1bd · c62a1bd
1 parent 304ef87
commit c62a1bd
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/R/FAI_BLUP.R b/R/FAI_BLUP.R
@@ -231,7 +231,7 @@ fai_blup <- function(.data,
   GID <- DEM[1:nrow(scores), (nrow(scores) + 1):nrow(sd.scores)]
   spatial.prob <- (1/GID)/(replicate(IN, c(as.numeric(apply((1/GID), 1, sum)))))
   ideotype.rank <- lapply(1:IN, function(i) {
-    sort(spatial.prob[, i], decreasing = TRUE)
+    sort(spatial.prob[, i] %>% replace_na(replace = 0), decreasing = TRUE)
   })
   names(ideotype.rank) <- paste("ID", 1:IN, sep = "")
   means.factor <- means[, names.pos.var.factor]