Bfastmonitor/Spatial: rmse extra band

[nephentes7]

Is there an easy way to get the RMSE (of the residuals) out, as an extra band (returnLayers) in Bfastmonitor? Many thanks!!!

[janverbesselt]

Hi, thanks for the question. Can you provide a reproducible (simple) example of what you want? rmse of the residuals of the model fit? it can certainly be derived see e.g. this example - and look at the bfastmonitor() code itself - for details - on how the current calculations are done. Try to look at this and see.

library(bfast)
library(zoo)
NDVIb <- as.ts(zoo(som$NDVI.b, som$Time))
plot(NDVIb)
monb <- bfastmonitor(NDVIb, start = c(2010, 13))
monb
bfastmonitor