From c5bcf51b952e809e4aa83b0a0ca3b0b8def629c5 Mon Sep 17 00:00:00 2001
From: Stephane Guindon <stephane.guindon@gmail.com>
Date: Tue, 24 Sep 2024 14:42:55 +0200
Subject: [PATCH] Hard-coded exp priors on IOU parameters

---
 src/iou.c  | 4 ++--
 src/mcmc.c | 1 +
 2 files changed, 3 insertions(+), 2 deletions(-)

diff --git a/src/iou.c b/src/iou.c
index 0266495c..ff43e32b 100644
--- 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
index 6236ad31..f5ee0f41 100644
--- 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);
         }