Jonas Gregorio de Souza jonas.gregorio@gmail.com

install.packages("spDates")

devtools::install_github("jgregoriods/spDates")

library(spDates)
data(neof)
data(centers)
jericho <- centers[centers$Site=="Jericho",]
model <- modelDates(neof, "C14Age", jericho, method="ols")
plot(model)

summary(model)
                   Time-versus-distance  Distance-versus-time
Start date:          9782 +/- 11 cal BP   11115 +/- 21 cal BP
Speed of advance: 0.97 +/- 0.0073 km/yr 0.64 +/- 0.0063 km/yr

rmamodel <- modelDates(neof, "C14Age", jericho, method="rma")
plot(rmamodel)

rmabins <- modelDates(neof, "C14Age", jericho, binWidth=500)
plot(rmabins)

iter <- iterateSites(neof, "C14Age", centers, "Site", binWidths=500)
iter$results

        r p bin  n              site
1: 0.9640 0 500  9           Aswad**
2: 0.9516 0 500  9          Cayönü**
3: 0.9496 0 500  8     Çatal Höyük**
4: 0.9488 0 500  8 Shillourokambos**
5: 0.9469 0 500  9         Jericho**
6: 0.9451 0 500  9     Qermez Dere**
7: 0.9429 0 500 10        Abu Madi**
8: 0.9316 0 500  9     Abu Hureyra**
9: 0.9188 0 500 11        Ali Kosh**

plot(iter$model)

plot(iter$map)

library(rcarbon)
mydates <- read.csv("myfile.csv")
mydates$cal <- calibrate(mydates$C14Age, mydates$C14SD)
mydates$med <- medCal(mydates$cal)