FBplot <-
function(fit, x = NULL, method ='MBD', depth = NULL, plot = TRUE,
prob = 0.5, color = 6, outliercol = 2, barcol = 4, fullout = FALSE,
factor = 1.5, xlim = c(1, nrow(fit)), ylim = c(min(fit) -
0.5 * diff(range(fit)), max(fit) + 0.5 * diff(range(fit))),
...){
#
# functional boxplot for functional data using
# method in Sun and Genton.
#
#
# fit is assumed to be an n-by-p matrix
# n= number of subjects
# p= number points where the function has been evaluated.
#
# rows with missing values are automatically removed.
#
library(fda)
elimna(fit)
fit=t(fit)
res=fbplot(fit, x = NULL, method = method, depth = depth, plot = plot,
prob =prob, color =color, outliercol =outliercol, barcol = barcol,
fullout = fullout,
factor = factor, xlim =xlim, ylim = ylim, ...)
res
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.