- rainbow::fboxplot(series_fds, "bivariate", "hdr", plotlegend=FALSE)
+ if (plot_bivariate)
+ rainbow::fboxplot(series_fds, "bivariate", "hdr", plotlegend=FALSE)
+}
+
+#' @title Functional boxplot on filaments
+#'
+#' @description Draw the functional boxplot on filaments obtained by \code{computeFilaments}
+#'
+#' @param data Object return by \code{getData}
+#' @param indices Indices as output by \code{computeFilaments}
+#'
+#' @export
+plotFilamentsBox = function(data, indices, ...)
+{
+ past_neighbs_indices = head(indices,-1)
+ plotFbox(data, function(i) i %in% past_neighbs_indices, plot_bivariate=FALSE)
+ par(new=TRUE)
+ # "Magic" found at http://stackoverflow.com/questions/13842560/get-xlim-from-a-plot-in-r
+ usr <- par("usr")
+ yr <- (usr[4] - usr[3]) / 27
+ plot(data$getSerie(tail(indices,1)), type="l", lwd=2, lty=2,
+ ylim=c(usr[3] + yr, usr[4] - yr), xlab="", ylab="")