Fix syntax error in bfast0n and recreate the documentation

DESCRIPTION

@@ -43,4 +43,4 @@ URL: http://bfast.R-Forge.R-project.org/
 LazyLoad: yes
 LazyData: yes
 LinkingTo: Rcpp
-RoxygenNote: 5.0.1
+RoxygenNote: 7.1.0

NAMESPACE

@@ -22,6 +22,7 @@ S3method(residuals,bfast01)
 export(bfast)
 export(bfast01)
 export(bfast01classify)
+export(bfast0n)
 export(bfastmonitor)
 export(bfastpp)
 export(bfastts)

R/bfast0n.R

@@ -21,19 +21,21 @@
 #' while also being able to deal with NA values by excluding them
 #' via \code{bfastpp}.
 #'
+#' @param ... Additional arguments to \code{\link[strucchange]{breakpoints}}.
 #' @inheritParams bfastpp
 #' @inheritParams strucchange::breakpoints
+#' @example examples/bfast0n.r
 #'
 #' @export bfast0n
 
-bfast0n <- function(data, order = 3,
+bfast0n <- function(data, formula=response ~ trend + harmon, order = 3,
                    lag = NULL, slag = NULL, na.action = na.omit,
                    stl = c("none", "trend", "seasonal", "both"),
-                   decomp=c("stlplus", "stl"), sbins=1), ...)
+                   decomp=c("stlplus", "stl"), sbins=1, ...)
 {
   data_pp <- bfastpp(data, order = 3,
                    lag = NULL, slag = NULL, na.action = na.omit,
                    stl = c("none", "trend", "seasonal", "both"),
-                   decomp=c("stlplus", "stl"), sbins=1))
-  return(strucchange::breakpoints(data = data_pp, ...))
+                   decomp=c("stlplus", "stl"), sbins=1)
+  return(strucchange::breakpoints(formula=formula, data = data_pp, ...))
 }

examples/bfast0n.r

@@ -0,0 +1,13 @@
+plot(simts) # stl object containing simulated NDVI time series
+datats <- ts(rowSums(simts$time.series))
+# sum of all the components (season,abrupt,remainder)
+tsp(datats) <- tsp(simts$time.series) # assign correct time series attributes
+plot(datats)
+
+# Detect breaks
+bp = bfast0n(datats)
+
+# Default method of estimating breakpoints
+bp[["breakpoints"]]
+# Custom method
+strucchange::breakpoints(bp, breaks=2)