R: Function to simulate a key innovation in macro-evolution with...
dd_KI_sim
R Documentation
Function to simulate a key innovation in macro-evolution with the innovative clade decoupling from the diversity-dependent diversification dynamics of the main clade
Description
Simulating a diversity-dependent diversification process where at a given time a new clade emerges with different inherent speciation rate and extinction rate and clade-level carrying capacity and with decoupled dynamics
Usage
dd_KI_sim(pars, age, ddmodel = 1)
Arguments
pars
Vector of parameters:
pars[1] corresponds to lambda_M (speciation rate of the main clade) pars[2] corresponds to mu_M (extinction rate of the main clade) pars[3] corresponds to K_M (clade-level carrying capacity of the main clade)
pars[4] corresponds to lambda_S (speciation rate of the subclade) pars[5] corresponds to mu_S (extinction rate of the subclade) pars[5] corresponds to K_S (clade-level carrying capacity of the subclade) pars[7] tinn, the time the shift in rates occurs in the lineage leading to the subclade
age
Sets the crown age for the simulation
ddmodel
Sets the model of diversity-dependence: ddmodel == 1 : linear dependence in speciation rate with parameter K (= diversity where speciation = extinction) ddmodel == 1.3 : linear dependence in speciation rate with parameter K' (= diversity where speciation = 0) ddmodel == 2 : exponential dependence in speciation rate with parameter K (= diversity where speciation = extinction) ddmodel == 2.1 : variant of exponential dependence in speciation rate with offset at infinity ddmodel == 2.2 : 1/n dependence in speciation rate ddmodel == 2.3 : exponential dependence in speciation rate with parameter x (= exponent) ddmodel == 3 : linear dependence in extinction rate ddmodel == 4 : exponential dependence in extinction rate ddmodel == 4.1 : variant of exponential dependence in extinction rate with offset at infinity ddmodel == 4.2 : 1/n dependence in extinction rate with offset at infinity
Value
out
A list with the following elements:
The first element is the tree of extant species in phylo format
The second element is the tree of all species, including extinct species, in phylo format
The third element is a matrix of all species where
- the first column is the time at which a species is born
- the second column is the label of the parent of the species; positive and negative values only indicate whether the species belongs to the left or right crown lineage
- the third column is the label of the daughter species itself; positive and negative values only indicate whether the species belongs to the left or right crown lineage
- the fourth column is the time of extinction of the species
If the fourth element equals -1, then the species is still extant.
- the fifth column indicates whether the species belong to the main clade (0) or the subclade (1)
The fourth element is the subclade tree of extant species (without stem)
The fifth element is the subclade tree of all species (without stem)
The sixth element is the same as the first, except that it has attributed 0 for the main clade and 1 for the subclade
The seventh element is the same as the Second, except that it has attributed 0 for the main clade and 1 for the subclade
The sixth and seventh element will be NULL if the subclade does not exist (because it went extinct).
Author(s)
Rampal S. Etienne
References
- Etienne, R.S. et al. 2012, Proc. Roy. Soc. B 279: 1300-1309, doi: 10.1098/rspb.2011.1439
- Etienne, R.S. & B. Haegeman 2012. Am. Nat. 180: E75-E89, doi: 10.1086/667574
Examples
dd_KI_sim(c(0.2,0.1,20,0.1,0.05,30,4),10)
Results
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(DDD)
Loading required package: deSolve
Attaching package: 'deSolve'
The following object is masked from 'package:graphics':
matplot
Loading required package: ape
Loading required package: ade4
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/DDD/dd_KI_sim.Rd_%03d_medium.png", width=480, height=480)
> ### Name: dd_KI_sim
> ### Title: Function to simulate a key innovation in macro-evolution with
> ### the innovative clade decoupling from the diversity-dependent
> ### diversification dynamics of the main clade
> ### Aliases: dd_KI_sim
> ### Keywords: models
>
> ### ** Examples
> dd_KI_sim(c(0.2,0.1,20,0.1,0.05,30,4),10)
$tes
Phylogenetic tree with 10 tips and 9 internal nodes.
Tip labels:
t5, t13, t12, t11, t2, t4, ...
Rooted; includes branch lengths.
$tas
Phylogenetic tree with 13 tips and 12 internal nodes.
Tip labels:
t1, t5, t13, t12, t11, t3, ...
Rooted; includes branch lengths.
$L
[,1] [,2] [,3] [,4] [,5]
[1,] 10.00000000 0 -1 3.508286 0
[2,] 10.00000000 -1 2 -1.000000 0
[3,] 7.43616338 -1 -3 1.498568 0
[4,] 6.66910649 2 4 -1.000000 1
[5,] 5.25943390 -1 -5 -1.000000 0
[6,] 4.50641870 -3 -6 3.942816 0
[7,] 1.80083728 4 7 -1.000000 1
[8,] 1.34464309 7 8 -1.000000 1
[9,] 0.69275728 4 9 -1.000000 1
[10,] 0.67940950 9 10 -1.000000 1
[11,] 0.59855553 -5 -11 -1.000000 0
[12,] 0.12729220 -5 -12 -1.000000 0
[13,] 0.02469041 -5 -13 -1.000000 0
$tesS
Phylogenetic tree with 5 tips and 4 internal nodes.
Tip labels:
[1] "t4" "t9" "t10" "t7" "t8"
Rooted; includes branch lengths.
$tasS
Phylogenetic tree with 5 tips and 4 internal nodes.
Tip labels:
[1] "t4" "t9" "t10" "t7" "t8"
Rooted; includes branch lengths.
$tes2
Phylogenetic tree with 10 tips and 9 internal nodes.
Tip labels:
t5, t13, t12, t11, t2, t4, ...
The tree includes a mapped, 2-state discrete character with states:
0, 1
Rooted; includes branch lengths.
$tas2
Phylogenetic tree with 13 tips and 12 internal nodes.
Tip labels:
t1, t5, t13, t12, t11, t3, ...
The tree includes a mapped, 2-state discrete character with states:
0, 1
Rooted; includes branch lengths.
>
>
>
>
> dev.off()
null device
1
>