Fix bug introduced in calc_het()

- dimnames(pr) for output of calc_genoprob gives chromosome-specific genotype names
kbroman · Aug 24, 2023 · fa8a362 · fa8a362
1 parent 480c535
commit fa8a362
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/R/calc_het.R b/R/calc_het.R
@@ -64,7 +64,7 @@ calc_het <-
     het_col <- vector("list", n_chr)
     geno <- dimnames(probs)[[2]]
 
-    if(any(nchar(geno) != 2)) {
+    if(any(nchar(unlist(geno)) != 2)) {
         stop("calc_het requires genotypes to have two-letter names")
     }