diff --git a/inst/scripts/phytools2.R b/inst/scripts/phytools2.R
index ca2f4cf..1390a13 100644
--- a/inst/scripts/phytools2.R
+++ b/inst/scripts/phytools2.R
@@ -107,7 +107,7 @@ input_matrix <- input_matrix[tree$tip.label,]
 ## ------------------------------------------------------------------------------------------------------------------------
 system.time({
     fit <- fitMk(
-        tree = tree, x = input_matrix, model = 'ER', pi = 'fitzjohn',
+        tree = tree, x = input_matrix, model = 'ER', pi = 'estimated',
         lik.func = 'pruning', logscale = TRUE
     )
 })