Description Usage Arguments Value Author(s) See Also Examples
Structural monitoring method for post-disturbance forest regrowth in a time series
1 2 3 4 |
x |
Time series. Can be an object of type |
change |
Numeric. Date of original disturbance. See |
dates |
Date. Vector with length |
startOffset |
Numeric. Number of (decimal) years by which to offset the start of the regrowth monitoring process. Can also be set to "floor", in which case the |
formula |
Formula to be fit to stable history period (see |
order |
Numeric. Order of the harmonic term of |
history |
Method or pre-determined period to use for history model (see |
level |
see |
w |
Numeric. Number of years before a regrowth flag is allowed to be assigned. |
plot |
Logical. Plot the results? |
ylabs |
Character. Vector of length 2 with y labels of the data time series and MOSUM plots, respectively |
object of class tsreg
Ben DeVries
regSpatial
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | data(ndmi)
plot(ndmi, type='b', cex=0.5)
# breakpoint detection with bfastmonitor
x <- bfastts(ndmi, time(ndmi), type="irregular")
bfm <- bfastmonitor(x, start=c(2005, 1), formula=response~harmon, order=1, plot=TRUE)
reg <- tsreg(ndmi, change = bfm$breakpoint, h = 0.5, plot = TRUE)
print(reg)
reg2 <- tsreg(ndmi, change = bfm$breakpoint, startOffset = "floor", h = 0.5, plot = TRUE)
print(reg2)
reg3 <- tsreg(ndmi, change = bfm$breakpoint, startOffset = "floor", h = 0.5, history='all', plot=TRUE)
print(reg3)
reg4 <- tsreg(ndmi, change = bfm$breakpoint, startOffset = "floor", h = 0.5, history='all', s=0, plot=TRUE)
print(reg4)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.