Skip to contents

Likelihood of C14 samples for JAGS models use in Age_OSLC14

ModelC14.Rd

A list of models for C14 data to define likelyhood in JAGS models.

Usage

data("ModelC14")

Format

This list contains:

full

a model considering error on calibration curve.

naive

a model not considering error on calibration curve.

References

Reimer PJ, Bard E, Bayliss A, Beck JW, Blackwell PC, Bronl Ramsey C, Buck CE, Cheng H, Edwards RL, Friedrich M, Grootes PM, Guilderson TP, Haflidason H, Hajdas I, Hatte C, Heaton TJ, Hoffmann DL, Hogg AG, Hughen KA, Kaiser KF, Kromer B,Manning SW, Niu M, Reimer RW, Richards DA, Scott EM, Southon JR, Staff RA, Turney CSM, van der Plicht J. 2013. IntCal13 ans Marine13 radiocarbon age calibration curves 0-50000 years cal BP. Radiocarbon 55(4)=1869-1887.

Hogg AG, Hua Q, Blackwell PG, Niu M, Buck CE, Guilderson TP, Heaton TJ, Palmer JG, Reimer PJ, Reimer RW, Turney CSM, Zimmerman SRH. 2013. SHCal13 Southern Hemisphere calibration, 0-50000 years cal BP. Radiocarbon 55(4):1889-1903

Examples

data(Model_AgeC14)
writeLines(Model_AgeC14$full)
#> model{
#>   # vraisemblance
#>   for(i in 1:N){
#>     X[i] ~ dnorm(mu[i], prec[i])
#>     mu[i] <- interp.lin(Age[i], xTableauCalib, yTableauCalib)
#>     Z[i]~dcat(c(0.1,0.9))
#>     err[i] <- interp.lin(Age[i], xTableauCalib, zTableauCalib)
#>     prec[i] <- 1/(alpha[i]^(-Z[i]+2)*(pow(sigma[i],2)+pow(err[i],2)))
#>   }
#>   # a priori
#>   Age[1]~dunif(xbound[1],xbound[2])
#>   invalpha[1]~dgamma(3,4)
#>   alpha[1]<-1/invalpha[1]
#>   for(j in 2:N){
#>     amin[j]<-max(StratiConstraints[1:j,j]*c(xbound[(2*(j-1)+1)],Age[1:(j-1)]))
#>     Age[j]~dunif(amin[j],xbound[2*j])
#>     invalpha[j]~dgamma(3,4)
#>     alpha[j]<-1/invalpha[j]
#>     }
#>   }