GmAMisc (version 1.1.1)

ppdPlot: 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

# NOT RUN {
#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")

# }

Run the code above in your browser using DataCamp Workspace