When there is drift, try going back to using the internal node simula…

…tion
beast-dev · Dec 15, 2023 · dab5926 · dab5926
1 parent cf0a9ac
commit dab5926
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/...r/evomodel/treedatalikelihood/preorder/MultivariateConditionalOnTipsRealizedDelegate.java b/...r/evomodel/treedatalikelihood/preorder/MultivariateConditionalOnTipsRealizedDelegate.java
@@ -151,7 +151,7 @@ private void simulateTraitForExternalNode(final int nodeIndex,
 
         } else {
 
-            if (zeroCount == dimTrait) { //  All missing completely at random
+            if (zeroCount == dimTrait && !likelihoodDelegate.getDiffusionProcessDelegate().hasDrift()) { //  All missing completely at random
                 //TODO: This is N(X_pa(j), l_j V_root). Why not N(X_pa(j), V_branch) ?
 
                 final double sqrtScale = Math.sqrt(1.0 / branchPrecision);