Computes rates of change from SPDs
Function for computing rates of change between abutting user-defined time-blocks.
spd2rc( spd, breaks = NULL, backsight = NULL, changexpr = expression((t1/t0)^(1/d) - 1) )
spd |
Summed Probability Distribution obtained using the |
breaks |
A vector giving the breakpoints between the time-blocks. |
backsight |
A single numeric value defining the distance in time between the focal year and the backsight year for computing the rate of change. |
changexpr |
An expression defining how the rate of change is calculated, where |
When the argument breaks is supplied the function aggregates the summed probability within each time-block and compared them across abutting blocks using the expression defined by changexpr. When the argument backsight is provided he expression is based on the comparison between the summed probability of each year and the associated backsight year.
An object of class spdRC.
## Not run: data(emedyd) caldates <- calibrate(x=emedyd$CRA, errors=emedyd$Error, normalised=FALSE, calMatrix=TRUE) bins <- binPrep(sites=emedyd$SiteName, ages=emedyd$CRA, h=50) emedyd.spd <- spd(caldates,bins,timeRange=c(16000,9000),runm=100) emedyd.gg <- spd2rc(emedyd.spd,breaks=seq(16000,9000,-1000)) emedyd.gg2 <- spd2rc(emedyd.spd,backsight=10) plot(emedyd.gg) plot(emedyd.gg2) ## End(Not run)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.