adolphi_tiepoint_simmer: Simulate from Adolphi tie-point distributions Returns the...

In eirikmn/bremla: Bayesian Regression Modeling of Layer-Counted Archives

Simulate from Adolphi tie-point distributions Returns the probability distributions for the tie-points given by Adolphi et al. (2018) and Muscheler et al. (2020)

Description

Simulate from Adolphi tie-point distributions Returns the probability distributions for the tie-points given by Adolphi et al. (2018) and Muscheler et al. (2020)

Usage

adolphi_tiepoint_simmer(
  nsims = 10000,
  tieshifts = numeric(5),
  plotdens = FALSE,
  x.ref = NULL,
  plothist = list(compute = FALSE, breaks = 50, col = "orange")
)

Arguments

nsims	Integer. The number of samples to be produced.
tieshifts	numeric which gives the amount each tie-point should be shifted, e.g. to go from describing offset from GICC05 to years before present.
plotdens	boolean that is passed on to adolphiloader. Set TRUE if you want to plot distributions
x.ref	numeric. Gives reference value on the x-axis when distributions are plotted.
plothist	list object describing if and how histogram of samples should be plotted

Value

returns a list containing the pdfs for each tie-point.

Author(s)

Eirik Myrvoll-Nilsen, [e-mail]

Examples

if(inlaloader()){
tieshifts= c(11050,12050,13050,22050,42050)
samples = adolphi_tiepoint_simmer(nsims=2000,tieshifts=tieshifts)
hist(samples[,3],col="orange",freq=0,breaks=50,main="Tie-point 3",xlab="Age (yb2k)")
}

eirikmn/bremla documentation built on Jan. 25, 2025, 4:41 a.m.