# Loading data
data(fakeexperiment)
data(faketree)
O <- new_aphylo(fakeexperiment[,2:3], tree = as.phylo(faketree))
# Baseline plot (all parameters but Pi)
plot_logLik(O)
# No psi parameter
plot_logLik(O ~ mu_d + Pi + eta)
Run the code above in your browser using DataLab