+ given through a correlation structure with the option
+ \code{corStruct}.
+
+ In the default setting (i.e., \code{method = "ML"} and \code{model =
+ "BM"}) the maximum likelihood estimation is done simultaneously on
+ the ancestral values and the variance of the Brownian motion process;
+ these estimates are then used to compute the confidence intervals in
+ the standard way (see the package \pkg{geiger} for a different
+ implementation). If \code{method = "pic"} or \code{"GLS"}, the
+ confidence intervals are computed using the expected variances under
+ the model, so they depend only on the tree.