if (FALSE) {
data(emedyd)
northernlevant <- emedyd[emedyd$Region=="2",]
bins <- binPrep(northernlevant$SiteName, northernlevant$CRA, h=50)
x <- calibrate(northernlevant$CRA, northernlevant$Error, normalised=FALSE)
spd.northernlevant <- spd(x, bins=bins, runm=50, timeRange=c(17000,8000))
plot(spd.northernlevant)
}
Run the code above in your browser using DataLab