Nothing
plot.bfastmonitor <- function(x, na.rm = TRUE, main = TRUE, ylab = "Data", ...)
{
if(isTRUE(main)) main <- if(is.na(x$breakpoint)) {
"No break detected"
} else {
sprintf("Break detected at: %i(%i)", floor(x$breakpoint),
round((x$breakpoint - floor(x$breakpoint)) * frequency(x$data)) + 1)
}
y <- if(is.null(dim(x$data))) x$data else x$data[,1L]
if(na.rm) y <- na.omit(as.zoo(y))
plot(y, type = "n", main = main, ylab = ylab, ...)
lines(window(y, end = x$history[2]), col = "black")
lines(window(y, start = x$history[1], end = x$history[2]),
col = "darkgreen", type = "p", pch = 19, cex = 0.5)
lines(window(y, start = x$monitor[1]), col = "red")
points(window(y, start = x$monitor[1]), col = "red", pch=19, cex=0.5) # new
test_pred <- predict(x$model, newdata = x$tspp)
test_pred <- zoo(test_pred, x$tspp$time, frequency = frequency(y))
lines(test_pred, col = "blue", lwd = 1.5)
abline(v = x$monitor[1], lty = 2, col = "black", lwd = 1)
abline(v = x$breakpoint, lty = 2, col = "red", lwd = 2)
legend("bottomleft", bty = "n",
c("Historical data", "New data", "Stable history", "Fit based on stable history", "Start of the Monitoring period", "Time of detected break"),
lty = c(1, 1, NA, 1, 2, 2),
col = c("black", "red", "darkgreen", "blue", "black", "red"),
pch = c(NA, NA, 19, NA, NA, NA)
)
invisible(x)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.