R: Adding Regimes to a Hansen Model

surfaceForward {surface}

R Documentation

Adding Regimes to a Hansen Model

Description

Carries out the forward phase of SURFACE's stepwise AIC routine, adding regime shifts to a Hansen model. addRegime performs one step of this analysis, and is called repeatedly by surfaceForward. At each step, the delta-AICc of each possible shift placement (i.e. branch) is calculated, and an updated Hansen model is returned with one shift added. This process is iterated until the model stops improving beyond a threshold delta-AICc

Usage

surfaceForward(otree, odata, starting_list=NULL, starting_shifts=NULL, 
exclude=0,  aic_threshold=0, max_steps=NULL, save_steps=FALSE, 
filename="temp_out_list.R", verbose=FALSE, plotaic=FALSE, 
error_skip=FALSE, sample_shifts=FALSE, sample_threshold=2)
addRegime(otree, odata, oldshifts, oldaic, oldfit, alloldaic=NULL, 
exclude=NULL, aic_threshold=0, verbose=FALSE, plotaic=FALSE, 
error_skip=FALSE, sample_shifts=FALSE, sample_threshold=2)

Arguments

`otree`	Phylogenetic tree in `ouchtree` format
`odata`	Data frame with rownames corresponding to `otree@labels`
`starting_list`	An optional list which may containing either a partially completed analysis (which can be built upon instead of starting over), or a custom starting model created with `startingModel`, which may include some pre-specified shifts
`starting_shifts`	An optional named character vector of shifts that are required to be in the Hansen model, which will be passed to `startingModel` when the initial model is built
`exclude`	Optionally, the proportion of the worst models (AICc scores for each shift point) to exclude in the current round (defaults to zero; values greater than 0.5 are not recommended)
`aic_threshold`	Change in AICc needed to accept a candidate model as a sufficient improvement over the previous iteration of SURFACE. Defaults to zero, meaning any improvement in the AICc will be accepted; more stringent thresholds are specified using negative values of `aic_threshold`
`max_steps`	Maximum number of regimes to allow to be added (assuming the model improvement continues to exceed `aic_threshold`)
`save_steps`	A logical indicating whether to save the current iteration of the model at each step (overwriting previous iterations) to a file `filename` (defaults to `FALSE`)
`filename`	Name of the file to save progress to at each step, if `savesteps=TRUE`
`verbose`	A logical indicating whether to print progress (defaults to `FALSE`)
`plotaic`	A logical indicating whether to plot AICc values of candidate models at each step (defaults to `FALSE`)
`error_skip`	A logical indicating whether to skip over any candidate model that produces an error message (this is rare, but can cause an entire analysis to abort; defaults to `FALSE`)
`sample_shifts`	A logical indicating whether to sample from among the best models at each step (those within `sample_threshold` of the best AICc), rather than always selecting the best candidate model (defaults to `FALSE`)
`sample_threshold`	Number of AICc units within which to sample among candidate models that are close to as good as the best model at each step (defaults to 2, but only used if `sample_shifts=TRUE`)
`oldshifts`	Any shifts present in the previous iteration of the Hansen model
`oldaic`	AICc value for the Hansen model from the previous iteration
`oldfit`	Previous fitted Hansen model
`alloldaic`	AICc values for each tested shift point in the previous iteration

Details

Can be time-consuming, as many likelihood searches are carried out at each iteration. Depending on the number of traits and taxa and the number of regimes that are fitted, surfaceForward can take anywhere from minutes to many hours (only tree sizes up to 128 taxa have been tested). Options to manage computation time include adding regimes one at a time with addRegime or using max_steps to perform the analysis several iterations at a time

Value

addRegime returns a list describing one iteration of the forward phase of the SURFACE analysis; surfaceForward returns a list of such lists consisting of each step of the stepwise process

`fit`	The fitted Hansen model selected for improving the AICc most over the previous iteration; consists of a single `hansentree` object if the number of traits `m = 1`, or a list of `hansentree` objects if `m > 1`
`all_aic`	The AICc for each model tested during the iteration (numbered by branch)
`aic`	The AICc of the current Hansen model
`savedshifts`	The shifts present in the current Hansen model; represented as a named character vector of regime shifts (lower-case letters), with names indicating branches containing shifts
`n_regimes`	A two-element vector of the number of regime shifts and the number of distinct regimes in the current model

Author(s)

Travis Ingram

References

Butler, M.A. & King, A.A. (2004) Phylogenetic comparative analysis: a modeling approach for adaptive evolution. American Naturalist 164: 683-695.

Ingram, T. & Mahler, D.L. (2013) SURFACE: detecting convergent evolution from comparative data by fitting Ornstein-Uhlenbeck models with stepwise AIC. Methods in Ecology and Evolution 4: 416-425.

Mahler, D.L., Ingram, T., Revell, L.J. & Losos, J.B. (2013) Exceptional convergence on the macroevolutionary landscape in island lizard radiations. Science 341: 292-295.

Examples

	## Not run: 
data(surfaceDemo)
tree<-surfaceDemo$tree
dat<-surfaceDemo$sim$dat
olist<-convertTreeData(tree,dat)
otree<-olist[[1]]; odata<-olist[[2]]
fwd<-surfaceForward(otree, odata, aic_threshold = 0, exclude = 0)
	
## End(Not run)

[Package surface version 0.5 Index]