# Some of the examples may run some time.
# Canadian weather data set
# There are three samples of mean temperatures for
# fifteen weather stations in Eastern Canada,
# another fifteen in Western Canada, and
# the remaining five in Northern Canada.
library(fda)
data_set <- t(CanadianWeather$dailyAv[,, "Temperature.C"])
k <- 3
gr_label <- rep(c(1, 2, 3), c(15, 15, 5))
# \donttest{
# Tukey's contrast matrix
h_tukey <- GFDmcv::contr_mat(k, type = "Tukey")
# testing without parallel computing
res <- multiFANOVA(data_set, gr_label, h_tukey)
summary(res, digits = 3)# }
# \dontshow{
h_tukey <- GFDmcv::contr_mat(k, type = "Tukey")
res <- multiFANOVA(data_set[, 1:10], gr_label, h_tukey, n_boot = 3)
summary(res, digits = 3)# }
Run the code above in your browser using DataLab