Description Usage Arguments Value Author(s) References Examples
View source: R/chronoEvents_Gauss.R
Bayesian modeling for combining Gaussian dates. These dates are assumed to be contemporaneous of the event date. The posterior distribution of the event date is sampled by MCMC algorithm as well as those of all parameters of the Bayesian model as described in Lanos & Philippe (2017).
1 2 3 | chronoEvents_Gauss(M, s, refYear=NULL, measurementsPerEvent, studyPeriodMin,
studyPeriodMax, numberChains = 2, numberAdapt = 10000, numberUpdate = 10000,
variable.names = c("theta"), numberSample = 50000, thin = 10)
|
M |
vector of measurement |
s |
vector of measurement errors |
refYear |
vector of year of reference for ages for coversion into calendar dates |
measurementsPerEvent |
vector containing the number of measurements associated with the first event, then the second ... |
studyPeriodMin |
numerical value corresponding to the start of the study period in BC/AD format |
studyPeriodMax |
numerical value corresponding to the end of the study period in BC/AD format |
numberChains |
number of Markov chains simulated |
numberAdapt |
number of iterations in the Adapt period of the MCMC algorithm |
numberUpdate |
number of iterations in the Update period of the MCMC algorithm |
variable.names |
names of the variables whose Markov chains are kept |
numberSample |
number of iterations in the Acquire period of the MCMC algorithm |
thin |
step between consecutive iterations finally kept |
This function returns a Markov chain of the posterior distribution. The MCMC chain is in date format BC/AD, that is the reference year is 0. Only values for the variables defined by 'variable.names' are given.
Anne Philippe & Marie-Anne Vibet
P. Lanos and A. Philippe. Hierarchical Bayesian Modeling for Combining Dates in Archaeological Context. Journal de la SFdS, Vol. 158 (2) pp 72-88 2017.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 | ### simulated data
# Number of events
Nevt = 3
# number of dates by events
measurementsPerEvent = c(2,3,2)
# positions
pos = 1 + c(0, cumsum(measurementsPerEvent) )
# simulation of data
theta.evt = seq(1,10, length.out= Nevt)
theta.evt[3] <-theta.evt[3] - 3 # stratigraphic inversion
theta = NULL
for(i in 1:Nevt ){
theta = c(theta, rep(theta.evt[i],measurementsPerEvent[i]))
}
s = seq(1,1, length.out= sum(measurementsPerEvent))
M=NULL
for( i in 1:sum(measurementsPerEvent)){
M= c(M, rnorm(1, theta[i], s[i] ))
}
sO2 = 1:Nevt
for (i in 1:Nevt) {
sO2[i]= 1/mean( 1/(s[pos[i]:(pos[i+1] -1)])^2 )
}
MCMCSample = chronoEvents_Gauss( M=M, s=s, measurementsPerEvent=measurementsPerEvent,
studyPeriodMin=-10, studyPeriodMax=30)
plot(MCMCSample)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.