Description Usage Arguments Details Value Author(s) References See Also

View source: R/surfaceForward.R

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

1 2 3 4 5 6 7 | ```
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)
``` |

`otree` |
Phylogenetic tree in |

`odata` |
Data frame with rownames corresponding to |

`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 |

`starting_shifts` |
An optional named character vector of shifts that are required to be in the Hansen model, which will be passed to |

`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 |

`max_steps` |
Maximum number of regimes to allow to be added (assuming the model improvement continues to exceed |

`save_steps` |
A logical indicating whether to save the current iteration of the model at each step (overwriting previous iterations) to a file |

`filename` |
Name of the file to save progress to at each step, if |

`verbose` |
A logical indicating whether to print progress (defaults to |

`plotaic` |
A logical indicating whether to plot AICc values of candidate models at each step (defaults to |

`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 |

`sample_shifts` |
A logical indicating whether to sample from among the best models at each step (those within |

`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 |

`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 |

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

`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 |

`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 |

Travis Ingram

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.

`surfaceBackward`

, `surfaceSimulate`

, `surfaceTreePlot`

, `surfaceSummary`

, `convertTreeData`

, `startingModel`

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.