X-Git-Url: https://git.auder.net/?p=talweg.git;a=blobdiff_plain;f=pkg%2Ftests%2Ftestthat%2Ftest-computeFilaments.R;h=7e1cafa7a7317f08878cc3d514769585b4778f4f;hp=355d58d921818d0befa82660d417d77dac51f4dc;hb=c36568faefc97bb417d355e2084398c1ad1acf92;hpb=3ddf1c12af0c167fe7d3bb59e63258550270cfc5 diff --git a/pkg/tests/testthat/test-computeFilaments.R b/pkg/tests/testthat/test-computeFilaments.R index 355d58d..7e1cafa 100644 --- a/pkg/tests/testthat/test-computeFilaments.R +++ b/pkg/tests/testthat/test-computeFilaments.R @@ -1,46 +1,12 @@ context("computeFilaments") -#shorthand: map 1->1, 2->2, 3->3, 4->1, ..., 149->2, 150->3 -I = function(i) - (i-1) %% 3 + 1 - -#MOCK data; NOTE: could be in inst/testdata as well -getDataTest = function(n) -{ - data = Data$new() - x = seq(0,9.5,0.1) - L = length(x) #96 1/4h - s1 = cos(x) - s2 = sin(x) - s3 = c( s1[1:(L%/%2)] , s2[(L%/%2+1):L] ) - #sum((s1-s2)^2) == 96 - #sum((s1-s3)^2) == 58 - #sum((s2-s3)^2) == 38 - s = list(s1, s2, s3) - series = list() - for (i in seq_len(n)) - { - serie = s[[I(i)]] + rnorm(L,sd=0.01) - level = mean(serie) - serie = serie - level - # 10 series with NAs for index 2 - if (I(i) == 2 && i >= 60 && i<= 90) - serie[sample(seq_len(L),1)] = NA - time = as.POSIXct(i*15*60, origin="2007-01-01", tz="GMT") - exo = runif(4) - exo_hat = runif(4) - data$append(time, serie, level, exo, exo_hat) - } - data -} - test_that("output is as expected on simulated series", { data = getDataTest(150) # index 143 : serie type 2 pred = computeForecast(data, 143, "Neighbors", "Zero", - horizon=length(data$getSerie(1)), simtype="endo", h_window=1) + horizon=length(data$getSerie(1)), simtype="endo", local=FALSE, h_window=1) f = computeFilaments(data, pred, 1, limit=60, plot=FALSE) # Expected output: 50-3-10 series of type 2, then 23 series of type 3 (closest next) @@ -63,7 +29,7 @@ test_that("output is as expected on simulated series", # index 142 : serie type 1 pred = computeForecast(data, 142, "Neighbors", "Zero", - horizon=length(data$getSerie(1)), simtype="endo", h_window=1) + horizon=length(data$getSerie(1)), simtype="endo", local=FALSE, h_window=1) f = computeFilaments(data, pred, 1, limit=50, plot=FALSE) # Expected output: 50-10-3 series of type 1, then 13 series of type 3 (closest next)