Dear r-sig-phylo,
I'd like to be able to estimate the likelihood of a model with a given
Brownian rate parameter under a given data set. For OU models, this can be
accomplished through the ouch package with the flag fit=FALSE. However
there does not seem to be a corresponding option for
I have code for this in my TreEvo package (on R-forge -- see
brownieneglnL in http://r-forge.r-project.org/plugins/scmsvn/viewcvs.php/pkg/R/brownie.R?rev=2root=treevoview=markup
) though note that the package itself is not all cleaned up into
install.packages()-ready form. You could also dig
Hi Carl (also Brian),
This can effectively be accomplished in the geiger function
fitContinuous() by setting the upper lower bounds for the parameter to
both be equal to the value for which you want to calculate the
likelihood. Unfortunately, these are not allowed to be exactly the
same;