Change or create new make_tday for handling other seasonality than daily
Allow the user to model other seasonalities than daily, e.g. week: make_tday2 <- function (time, kseq, tstep = 3600, units = "hours") { tday <- sapply(kseq, function(k) { tk <- time + k * tstep as.numeric(tk - as.POSIXct(as.character(cut(tk, "week")), tz = "UTC"), unit = "hours") }) nams(tday) <- paste0("k", kseq) return(as.data.frame(tday)) }