You will be presented with the full text of a biomedical research paper. Please write an abstract for this paper. Your response should include the abstract and no additional text.

Paper text:
macro - evolutionary studies frequently use phylogenetic trees to examine diversification rate variation within and among groups .
variation in diversification rate in a phylogenetic tree can inform evolutionary hypotheses regarding the role of past geological or climatic events , evolutionary novelty , and adaptive and non - adaptive radiations.13 various methods have been developed to examine diversification rate variation , including , but not limited to , those that examine the accumulation of lineages through time,47 tree balance,810 distribution of tree branch lengths,11 and the shape of ordered cladogenic events.12 here , we present the r library iterates , which implements the parametric rate comparison ( prc ) test,13 a new method and approach for identifying rate variation in phylogenetic trees . an in - depth description of the prc , as well as an examination of its statistical power and false positive rates , can be found in the study by shah et al.13 briefly , the prc examines the fit of a distribution of branch lengths extracted from a phylogenetic tree to standard statistical distributions.14 the lengths of terminal branches are treated as censored at the time of sampling .
internal branch lengths and terminal branch lengths are jointly modeled using the censored form of a given distribution .
this approach differs from other approaches aimed at identifying rate heterogeneity in a phylogenetic tree in that it does not attempt to estimate the parameters of a particular model of diversification , rather it simply examines the statistical properties of the distribution of branch lengths . by iterating through subtrees in a tree
, the prc can be used to identify subclades where diversification rate differs from the remainder of the tree ( see fig . 1 in shah et al.13 ) .
the prc can also be used to compare a priori defined groups , and does not require that these groups be monophyletic .
the prc can be used both as a hypothesis testing tool and for exploratory data analysis using functions in iterates library .
the function iterates through subtrees of a phylogenetic tree and compares the distribution of branch lengths in that subtree to the remainder of the tree .
the function has arguments that allow flexibility for the user to govern how 
 comp.subs implements the prc .
an example usage of 
 comp.subs showing some important arguments at their default values is as follows : here , the argument 
 tree is an ultrametric phylogenetic tree of object class 
 phylo .
the argument thr indicates the threshold for the minimum number of branches required for a subtree to be considered in the analysis .
the argument 
 srt determines how the branch that links a subtree to the remainder of the tree is treated in the analysis .
the default is to drop the linking branch from the analysis because there is no way of knowing exactly at what point an inferred rate change might have occurred along the branch ; however , the user has the option to include this linking branch as part of the subtree or part of the remainder of the tree .
there are four different statistical distributions that 
 iterates uses to model the distribution of branch lengths : exponential , weibull , log - normal , and variable rates.13,14 the argument 
 mod.id is used to indicate which of these distributions are used in the prc . the default is to consider only an exponential model . when multiple distributions are included in the analysis , 
 comp.subs will use the akaike information criterion ( aic)15 scores to pick the best - fit model for each subtree vs. remainder of the tree comparison .
there is an expectation that false - positive rates will increase as the tree under scrutiny deviates from a pure - birth model .
a whole - tree randomization test can be used to avoid spurious inferences based on observed statistical significance for each node . here
, the topology of the tree is randomized while retaining the observed branching times . for each randomization ,
the prc test is employed and the number of statistically significant clades is recorded . following many randomizations ,
the observed number of significant clades is compared against the distribution of the number of significant clades from the randomized trees .
this function requires an ultrametric tree of object class 
 phylo and has arguments that allow the user to determine the number of randomizations and the statistical distributions used .
the results of the prc can be visualized on a phylogenetic tree using the function 
 color.tree.plot .
various options are available to the user to illustrate the relative direction , magnitude , and statistical support for a diversification rate change .
the function requires a tree object class phylo and the result object from 
 comp.subs .
figure 1 shows an example of a phylogenetic tree showing regions with rate variation using default settings .
the prc test assumes that taxon sampling is complete , although it is robust to incomplete taxon sampling if taxon sampling is random.13 for situations where there is ambiguity as to what complete taxon sampling might be , for example if there is a recent radiation , or when a tree is based on a higher taxonomic rank ( eg , a genus - level phylogeny ) , the user might choose to trim a particular amount of time from the tree .
for example , a researcher might decide that taxon sampling is complete up until five million years before present in the phylogeny .
the function 
 trimtree can be used to trim a specified amount of time ( or branch length ) from the tips of a tree .
this function returns a list that contains the trimmed tree and a key indicating the taxa from the original tree that have been collapsed to each terminal node of the trimmed tree .
the k - clades prc allows for statistical hypothesis testing of rate variation among a priori defined clades in a tree .
an example usage of 
 comp.fit.subs showing some important arguments at their default values is as follows : here , 
 tree.kclades$subtree is the list of all subtrees provided by 
 id.subtrees .
the argument k indicates the maximum number of different rates to explore , in this case three .
the function will compare all possible models of rate variation , ranging from all subtrees having the same rate , to k number of subtrees having different rates .
as in the function 
 comp.subs , the user can choose any combination of the four available statistical distributions using the argument 
 mod.id .
the results can be summarized using the function 
 tab.summary , which will restrict the returned output to a aic limit determined by the user and the best - fit model for each k. an example of 
 tab.summary and output is as follows : 
 > tab.summary(cfs,daic = 10 ) in this example , the best model is one that groups the first and second subtrees ( subtrees 55 and 27 , respectively ) to share the same rate , and the third subtree ( subtree 3 ) to be modeled as having a separate rate .
note , based on aic , modeling each subtree as having a separate rate also falls within the group of  best  models .
the 
 tab.summary function can be useful when exploring rate variation across numerous subtrees because the number of possible groupings can get quite large ; for example , when k = 5 , there are 52 different groupings explored .
the function iterates through subtrees of a phylogenetic tree and compares the distribution of branch lengths in that subtree to the remainder of the tree .
the function has arguments that allow flexibility for the user to govern how 
 comp.subs implements the prc .
an example usage of 
 comp.subs showing some important arguments at their default values is as follows : here , the argument 
 tree is an ultrametric phylogenetic tree of object class 
 phylo .
the argument thr indicates the threshold for the minimum number of branches required for a subtree to be considered in the analysis .
the argument 
 srt determines how the branch that links a subtree to the remainder of the tree is treated in the analysis .
the default is to drop the linking branch from the analysis because there is no way of knowing exactly at what point an inferred rate change might have occurred along the branch ; however , the user has the option to include this linking branch as part of the subtree or part of the remainder of the tree .
there are four different statistical distributions that 
 iterates uses to model the distribution of branch lengths : exponential , weibull , log - normal , and variable rates.13,14 the argument 
 mod.id is used to indicate which of these distributions are used in the prc . the default is to consider only an exponential model . when multiple distributions are included in the analysis , 
 comp.subs will use the akaike information criterion ( aic)15 scores to pick the best - fit model for each subtree vs. remainder of the tree comparison .
there is an expectation that false - positive rates will increase as the tree under scrutiny deviates from a pure - birth model .
a whole - tree randomization test can be used to avoid spurious inferences based on observed statistical significance for each node . here
, the topology of the tree is randomized while retaining the observed branching times . for each randomization ,
the prc test is employed and the number of statistically significant clades is recorded . following many randomizations ,
the observed number of significant clades is compared against the distribution of the number of significant clades from the randomized trees .
this function requires an ultrametric tree of object class 
 phylo and has arguments that allow the user to determine the number of randomizations and the statistical distributions used .
the results of the prc can be visualized on a phylogenetic tree using the function 
 color.tree.plot .
various options are available to the user to illustrate the relative direction , magnitude , and statistical support for a diversification rate change .
the function requires a tree object class phylo and the result object from 
 comp.subs .
figure 1 shows an example of a phylogenetic tree showing regions with rate variation using default settings .
the prc test assumes that taxon sampling is complete , although it is robust to incomplete taxon sampling if taxon sampling is random.13 for situations where there is ambiguity as to what complete taxon sampling might be , for example if there is a recent radiation , or when a tree is based on a higher taxonomic rank ( eg , a genus - level phylogeny ) , the user might choose to trim a particular amount of time from the tree .
for example , a researcher might decide that taxon sampling is complete up until five million years before present in the phylogeny .
the function 
 trimtree can be used to trim a specified amount of time ( or branch length ) from the tips of a tree .
this function returns a list that contains the trimmed tree and a key indicating the taxa from the original tree that have been collapsed to each terminal node of the trimmed tree .
the k - clades prc allows for statistical hypothesis testing of rate variation among a priori defined clades in a tree .
this function returns a list containing all the subtrees that might be compared . the function 
 comp.fit.subs implements the k - clades prc test .
an example usage of 
 comp.fit.subs showing some important arguments at their default values is as follows : here , 
 tree.kclades$subtree is the list of all subtrees provided by 
 id.subtrees .
the argument k indicates the maximum number of different rates to explore , in this case three . the function will compare all possible models of rate variation , ranging from all subtrees having the same rate , to k number of subtrees having different rates . as in the function 
 comp.subs , the user can choose any combination of the four available statistical distributions using the argument 
 mod.id .
the results can be summarized using the function 
 tab.summary , which will restrict the returned output to a aic limit determined by the user and the best - fit model for each k. an example of 
 tab.summary and output is as follows : 
 > tab.summary(cfs,daic = 10 ) in this example , the best model is one that groups the first and second subtrees ( subtrees 55 and 27 , respectively ) to share the same rate , and the third subtree ( subtree 3 ) to be modeled as having a separate rate .
note , based on aic , modeling each subtree as having a separate rate also falls within the group of  best  models
. the 
 tab.summary function can be useful when exploring rate variation across numerous subtrees because the number of possible groupings can get quite large ; for example , when k = 5 , there are 52 different groupings explored .
the package 
 iterates provides the library functions required to employ the prc test and the k - clades prc as described by shah et al.13 in the r statistical computing environment ( r development core team 2011).16 the package and dependences are available at the comprehensive r archive network ( http://cran.r-project.org/ ) .
the functions therein will be useful to investigators exploring diversification rate variation across a phylogeny and conducting explicit hypothesis testing among subtrees in a phylogeny .
patterns of diversification rate variation detected in phylogenetic hypotheses are frequently used to infer historical , ecological , and evolutionary processes . 
 the parametric rate comparison ( prc ) is a method for detecting rate variation in trees that models branch lengths as random variables drawn from familiar statistical distributions . 
 iterates is a library of functions for the r statistical computing environment for implementing prc on phylogenetic trees . here 
 , we describe some of the functions in iterates for subtree identification , tree manipulation , applying the prc and k - clades prc analyses , and conducting a whole - tree randomization test .