likSurface.rateshift plots the likelihood surface (which is sometimes quite rugged).
rateshift(tree, x, nrates=1, niter=10, method="ML", ...)
"plot"(x, ...)
likSurface.rateshift(tree, x, nrates=2, shift.range=NULL,
density=20, plot=TRUE, ...)"phylo".names(x) should contain the species names and match tree$tip.label. For plot method, x is an object of class "rateshift"."ML" (maximum likelihood) or "REML" (restricted maximum likelihood).plot method, these will be passed to plotSimmap. For rateshift, optional arguments include: tol, tolerance; plot & print, logical values indicating whether to plot or print the progress of the optimization (default to FALSE); quiet, logical argument indicating whether to suppress all notifications (defaults to FALSE); minL, numeric value; and fixed.shift, either a vector of fixed shift points, or a logical value.likSurface.rateshift.likSurface.rateshift.likSurface.rateshift. If plot=FALSE then the surface is returned."rateshift", or, in the case of likSurface.rateshift, a likelihood surface for the shift points.
brownie.lite