Description Usage Arguments Details Value Author(s) Examples
View source: R/combination_Gauss.R
Simple Bayesian modeling for combining Gaussian dates with known variance. These dates are assumed to be contemporaneous of the target date and have non identical distributions as the variance may be different for each date. The posterior distribution of the modeling is sampled by a MCMC algorithm implemented in JAGS.
1 2 3 | combination_Gauss(M, s, refYear=NULL, 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 |
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 |
If there are Nbobs measurements M associated with their error s, the model is the following one :
for j in (1:Nbobs), Mj ~ N(theta, sj^2)
theta ~ U(ta, tb)
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
1 2 3 4 5 | data(sunspot)
MCMC = combination_Gauss(M=sunspot$Age[1:10], s= sunspot$Error[1:10], refYear=rep(2016,10),
studyPeriodMin=900, studyPeriodMax=1500, variable.names = c('theta'))
plot(MCMC)
gelman.diag(MCMC)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.