Uncertainty in estimate of lambda illustrated by simulating 1000 data sets on the Geopsiza tree under a model with lambda = 0.62. Shows the high uncertainty in estimating lambda at all, despite the power in the estimate at this value.
(Code)
Will fill this out further later. Will be interesting to compare at different base values of lambda used for the simulation.
Have updated the codbase to provide bootstraps of parameters as part of the monte carlo output. Requires another generic function introduction: getParameters (in addition to getLikelihood), will probably make optional to keep method simple and general. See the PMC package code.
Inference in larger trees
Two days later… Dec 1st: I run roughly the same strength of phylogenetic signal but with a much larger (simulated) tree, with 354 tips. Here the estimate is at least centered around the true value and contains none of the zero estimates of the geospiza tree.
