# NOT RUN {
# Load in data
data(TestChronData)
data(TestRSLData)
# Run through Bchronology
RSLrun <- with(TestChronData, Bchronology(
ages = ages,
ageSds = ageSds,
positions = position,
positionThicknesses = thickness,
ids = id,
calCurves = calCurves,
predictPositions = TestRSLData$Depth
))
# Now run through BchronRSL
RSLrun2 <- BchronRSL(RSLrun, RSLmean = TestRSLData$RSL, RSLsd = TestRSLData$Sigma, degree = 3)
# Summarise it
summary(RSLrun2)
# Plot it
plot(RSLrun2)
# }
Run the code above in your browser using DataLab