- # TODO: 3 == magic number
- random_var = matrix(nrow=3, ncol=48)
- for (mc in seq_len(nrow(random_var)))
- {
- random_indices = sample(fdays, length(indices))
- random_series = t( sapply(random_indices, function(i) {
- c( data$getSerie(i), data$getSerie(i+1) )
- }) )
- random_var[mc,] = apply(random_series, 2, sd)
- }
- random_var = apply(random_var, 2, median)
-
- yrange = range(ref_var, random_var)