Description Usage Arguments Details Value Author(s) References See Also Examples
This function provides an abridged output of results obtained from the comp.fit.subs
function by restricting the output to a user provided delta AIC threshold.
1 | tab.summary(res, daic = 2, show.rate = FALSE)
|
res |
A dataframe obtained from |
daic |
A value indicating a threshold of delta AIC relative to the best fit model for each k to be included in the output. |
show.rate |
A boolean indicting whether the rate parameters are included in the output. |
This function will provide a reduced output of the results provided by the comp.fit.subs
function by allowing the user to choose a critical delta AIC for each value of k that determines which comparisons are included in the output. The best fit model for each k is included in the output regardless of delta AIC. The show.rate argument indicates whether the rate estimate for each of the subtrees is included in the output.
A dataframe that consists of the following:
k |
The number of groups |
Groups |
the groupings for each subtree numbered as 1 to the number of subtrees indicated. The numbering corresponds to the order in which subtrees are identified by focal. Groups are separated with 'vs.'. |
gi_rate |
The rate for the ith group in the analysis. |
LL |
The log likelihood for the entire model for a grouping scheme. |
AIC |
Akaike information criterion score for the entire model for a grouping scheme. |
AICc |
Akaike information criterion corrected for sample size. |
dAICc |
The delta AIC across all grouping schemes and k values relative to the best fit model. |
Premal Shah, Benjamin Fitzpatrick and James Fordyce.
Shah, P., B. M. Fitzpatrick, and J. A. Fordyce. 2013. A parametric method for assessing diversification rate variation in phylogenetic trees. Evolution 67:368-377.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 | ## Not run:
data(hivtree.newick)
cat(hivtree.newick, file = "hivtree.phy", sep = "\n")
tree.hiv <- read.tree("hivtree.phy") # load tree
unlink("hivtree.phy") # delete the file "hivtree.phy"
idHIV<-id.subtrees(tree.hiv)
plot(idHIV$tree,show.node.label=TRUE)
cfsHIV<-comp.fit.subs(idHIV$subtree,focal=c(153,119,96,5),k=4)
tab.summary(cfsHIV)
tab.summary(cfsHIV,daic=1)
tab.summary(cfsHIV,daic=0.01)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.