rimov: Year temperature curves

Description Usage Format Details References Examples

Description

Year temperature curves

Usage

1

Format

A curve_set object with water temperatures in 365 days of the 36 years. The component curve_set[['r']] is a vector of days (from 1 to 365), whereas curve_set[['obs']] contains the water temperatures such that each column gives year temperatures in a year.

Details

The water temperature data sampled at the water level of Rimov reservoir in Czech republic every day for the 36 years between 1979 and 2014.

References

Mrkvička, T., Hahn, U. and Myllymäki, M. (2016) A one-way ANOVA test for functional data with graphical interpretation. arXiv:1612.03608 [stat.ME] (http://arxiv.org/abs/1612.03608)

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
# This is an example analysis of the water temperature data set in Mrkvicka et al. (2017).
data(rimov)
groups <- factor(c(rep(1, times=12), rep(2, times=12), rep(3, times=12)))
# Plot data in groups
subs <- function(group, ...) {
  cset <- rimov
  cset$obs <- rimov$obs[, groups == group]
  plot(cset, ...)
}
for(i in 1:3) subs(i, main=paste("group ", i, sep=""), ylab="Temperature")

system.time(res <- graph.fanova(nsim=2499, curve_set=rimov, groups=groups,
                                summaryfun="means", type="erl"))
plot(res)
system.time(res2 <- graph.fanova(nsim=2499, curve_set=rimov, groups=groups,
                                 summaryfun="contrasts", type="erl"))
plot(res2)

myllym/GET documentation built on Sept. 30, 2018, 5:49 a.m.