recover one line

beiko-lab · Feb 21, 2024 · ebd4bea · ebd4bea
1 parent c137dd3
commit ebd4bea
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/R/EstimateCCM.R b/R/EstimateCCM.R
@@ -113,6 +113,7 @@ EstimateCCM <- function(profiles, phytree, ip=0.1, pen=0.5,  ...){
                     invGAMMA %*% liks[des2, ]
             v <- v.l * v.r
             comp[anc] <- sum(v)
+	    liks[anc, ] <- v/comp[anc]
         }
         nonNeg_comp <- ifelse(comp <= 0, 1e-13, comp)
         dev <- -1 * sum(log(nonNeg_comp[-TIPS]))