# This example may run some time.
# gait data (the first feature)
library(fda)
gait.data.frame <- as.data.frame(gait)
x.gait <- as.matrix(gait.data.frame[, 1:39])
# vector of group labels
group.label.gait <- rep(1:3, each = 13)
# \donttest{
# all FANOVA tests with default parameters
set.seed(123)
fanova1 <- fanova.tests(x.gait, group.label.gait)
print(fanova1)# }
# \dontshow{
# only three tests with non-default parameters
fanova2 <- fanova.tests(x.gait, group.label.gait,
test = c("FP", "GPF", "Fmaxb"),
params = list(paramFP = list(int = c(0.025, 0.975),
B.FP = 3, basis = "b-spline",
criterion = "eBIC",
commonK = "mean",
minK = 5, maxK = 19,
norder = 4, gamma.eBIC = 0.7),
paramFmaxb = 3))
print(fanova2)# }
Run the code above in your browser using DataLab