Hard-coded exp priors on IOU parameters

stephaneguindon · Sep 24, 2024 · c5bcf51 · c5bcf51
1 parent 2a420d9
commit c5bcf51
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 2 deletions.
diff --git a/src/iou.c b/src/iou.c
@@ -160,7 +160,7 @@ phydbl IOU_Prior_Theta(t_tree *tree)
   phydbl lbda,lnP;
 
   lnP = 0.0;
-  lbda = 1.0E+4;
+  lbda = 1.0E+2;
   lnP += Log_Dexp(tree->mmod->ou_theta,lbda);
   return(lnP);
 }
@@ -173,7 +173,7 @@ phydbl IOU_Prior_Mu(t_tree *tree)
   phydbl lbda,lnP;
 
   lnP = 0.0;
-  lbda = 1.0E+3;
+  lbda = 1.0E+0;
   for (int i = 0; i < tree->mmod->n_dim; i++)
       lnP += Log_Dexp(tree->mmod->ou_mu[i], lbda);
   return (lnP);

diff --git a/src/mcmc.c b/src/mcmc.c
@@ -14689,6 +14689,7 @@ void MCMC_PHYREX_IOU_Update_Theta(t_tree *tree)
 
             tree->mcmc->run++;
 
+            PHYREX_Print_MCMC_Stats(tree);
             PHYREX_Print_MCMC_Tree(tree);
             PHYREX_Print_MCMC_Summary(tree);
         }