ppdPlot: R function for plotting Posterior Probability Densities for...

In GmAMisc: 'Gianmarco Alberti' Miscellaneous

R function for plotting Posterior Probability Densities for Bayesian modeled 14C dates/parameters

Description

The function allows plot Posterior Probability Densities with a nice outlook thanks to 'ggplot2'.
It takes as input a dataframe that must be organized as follows (it is rather easy to do that once the data have been exported from OxCal):
-calendar dates (first column to the left);
-posterior probabilities (second column);
-grouping variables (third column), which could contain the names of the events of interest (e.g., phase 1 start, phase 1 end, phase 2 start, phase 2 end, etc).

Usage

ppdPlot(data, lower = min(data[, 1]), upper = max(data[, 1]), type)

Arguments

data	Dataframe containing the data as returned by the OxCal program.
lower	Lower limit of the calendar date axis.
upper	Upper limit of the calendar date axis; if the lower and upper parameters are not provided, the default values will be the earliest and latest calendar dates.
type	Type of plot the user wishes to plot (a: curves outlined by a line; b: curves plotted as solid areas; c: combination of a and b).

Examples

#load a toy dataset
data(radioc_data)

#plot the Posterior Probability Densities for the phases' parameters
ppdPlot(radioc_data, type="a")

#plot the Posterior Probability Densities for the phases' parameters,
# setting different boundaries for the x-axis and using filled curves instead of simple outlines
ppdPlot(radioc_data, -1000, 100, type="b")

GmAMisc documentation built on March 18, 2022, 6:32 p.m.