tab.summary: tab.summary

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

Description

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.

Usage

1
tab.summary(res, daic = 2, show.rate = FALSE)

Arguments

res

A dataframe obtained from comp.fit.subs function.

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.

Details

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.

Value

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.

Author(s)

Premal Shah, Benjamin Fitzpatrick and James Fordyce.

References

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.

See Also

tab.summary id.subtrees

Examples

 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)

iteRates documentation built on May 2, 2019, 9:12 a.m.