Hard-coded exp priors on IOU parameters

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

Hide file tree

Showing 2 changed files with 1 addition and 7 deletions.
diff --git a/src/init.c b/src/init.c
@@ -3624,7 +3624,7 @@ void PHYREX_Set_Default_Migrep_Mod(int n_otu, t_phyrex_mod *t)
   t->max_sigsq         = 1.E+3;
   t->prior_param_sigsq = 10.0;
   assert(t->n_dim > 0);
-  for(int i=0;i<t->n_dim;++i) t->sigsq[i] = 1.E-1;
+  for(int i=0;i<t->n_dim;++i) t->sigsq[i] = 10.;
 
   t->nu = 1.0E-0;
 

diff --git a/src/phyrex.c b/src/phyrex.c
@@ -897,12 +897,6 @@ void PHYREX_XML(char *xml_filename)
   MCMC_Randomize_Veloc(mixt_tree);
   MCMC_Randomize_Contmod(mixt_tree);
 
-//   mixt_tree->mmod->ou_theta = 0.1;
-//   mixt_tree->mmod->ou_mu[0] = 0.1;
-//   mixt_tree->mmod->ou_mu[1] = 0.1;
-//   mixt_tree->mmod->sigsq[0] = 1.E-3;
-//   mixt_tree->mmod->sigsq[1] = 1.E-3;
-
   Set_Ignore_Root(YES,mixt_tree);
   Set_Bl_From_Rt(YES,mixt_tree);