posterior3p {EvoPhylo} | R Documentation |

An example dataset of posterior parameter samples resulting from a clock-based Bayesian inference analysis using the skyline fossilized birth–death process (FBD) tree model with Mr. Bayes after combining all parameter (.p) files into a single data frame with the `combine_log`

function. This particular example was produced by analyzing the data set with three morphological partitions from Simões & Pierce (2021).

`data("posterior3p")`

A data frame with 4000 observations on several variables estimated for each generation during analysis. The number of variables depends on parameter set up, but for clock analyses with Mr. Bayes, will typically include the following:

`Gen`

A numeric vector for the generation number

`LnL`

A numeric vector for the natural log likelihood of the cold chain

`LnPr`

A numeric vector for the natural log likelihood of the priors

`TH.all.`

A numeric vector for the total tree height (sum of all branch durations, as chronological units)

`TL.all.`

A numeric vector for total tree length (sum of all branch lengths, as accumulated substitutions/changes)

`prop_ancfossil.all.`

A numeric vector indicating the proportion of fossils recovered as ancestors

`sigma.1.`

,`sigma.2.`

,`sigma.3.`

A numeric vector for the standard deviation of the lognormal distribution governing how much rates vary across characters for each data partition

`m.1.`

,`m.2.`

,`m.3.`

A numeric vector for the rate multiplier parameter for each data partition

`net_speciation_1.all.`

,`net_speciation_2.all.`

,`net_speciation_3.all.`

,`net_speciation_4.all.`

A numeric vector for net speciation estimates for each time bin

`relative_extinction_1.all.`

,`relative_extinction_2.all.`

,`relative_extinction_3.all.`

,`relative_extinction_4.all.`

A numeric vector for relative extinction estimates for each time bin

`relative_fossilization_1.all.`

,`relative_fossilization_2.all.`

,`relative_fossilization_3.all.`

,`relative_fossilization_4.all.`

A numeric vector for relative fossilization estimates for each time bin

`tk02var.1.`

,`tk02var.2.`

,`tk02var.3.`

A numeric vector for the variance on the base of the clock rate for each clock partition

`clockrate.all.`

A numeric vector for the base of the clock rate

Datasets like this one can be produced from parameter log (.p) files using `combine_log`

. The number of variables depends on parameter set up, but for clock analyses with Mr. Bayes, will typically include the ones above, possibly also including an `alpha`

for each partition, which contains the shape of the gamma distribution governing how much rates vary across characters (when shape of the distribution is unlinked across partitions). When using the traditional FBD model rather than the skyline FBD model used to produce this dataset, there will be only one column for each of `net_speciation`

, `relative_extinction`

and `relative_fossilization`

. When using a single morphological partition, different columns may be present; see `posterior1p`

for an example with just one partition.

Simões, T. R. and S. E. Pierce (2021). Sustained High Rates of Morphological Evolution During the Rise of Tetrapods. *Nature Ecology & Evolution* 5: 1403–1414.

`posterior1p`

for an example dataset of posterior parameter samples resulting from an analysis with 1 partition rather than 3.

[Package *EvoPhylo* version 0.3.2 Index]