if (FALSE) {
#Load EUROEVOL Data
data(euroevol)
#Subset Danish Dates
denmark <- subset(euroevol,Country=="Denmark")
#Calibrate and Bin
denmarkDates <- calibrate(x=denmark$C14Age,errors=denmark$C14SD)
denmarkBins <- binPrep(sites=denmark$SiteID,ages=denmark$C14Age,h=200) #200 years bin size
#Compute median date for each bin
bm <- binMed(x=denmarkDates,bins=denmarkBins)
#Compute median date for each date
dm <- medCal(denmarkDates)
#Compute SPD
denmarkSPD <- spd(x=denmarkDates,bins=denmarkBins,timeRange=c(10000,4000))
#Plot SPD and barCodes of median dates
plot(denmarkSPD,runm=200)
barCodes(dm,yrng=c(0,0.01))
#Plot SPD and barCodes of median bins in BC/AD
plot(denmarkSPD,runm=200,calendar="BCAD")
barCodes(BPtoBCAD(bm),yrng=c(0,0.01))
}
Run the code above in your browser using DataLab