Poniższy kod ocenia podobieństwo między dwoma szeregami czasowymi:
set.seed(10)
RandData <- rnorm(8760*2)
America <- rep(c('NewYork','Miami'),each=8760)
Date = seq(from=as.POSIXct("1991-01-01 00:00"),
to=as.POSIXct("1991-12-31 23:00"), length=8760)
DatNew <- data.frame(Loc = America,
Doy = as.numeric(format(Date,format = "%j")),
Tod = as.numeric(format(Date,format = "%H")),
Temp = RandData,
DecTime = rep(seq(1, length(RandData)/2) / (length(RandData)/2),
2))
require(mgcv)
mod1 <- gam(Temp ~ Loc + s(Doy) + s(Doy,by = Loc) +
s(Tod) + s(Tod,by = Loc),data = DatNew, method = "ML")
Tutaj gam
służy do oceny, jak temperatura w Nowym Jorku i Miami różni się od średniej temperatury (z obu lokalizacji) w różnych porach dnia. Problem, który mam teraz, polega na tym, że muszę dołączyć termin interakcji, który pokazuje, jak temperatura w każdej lokalizacji zmienia się w ciągu dnia w różnych dniach roku. W końcu mam nadzieję wyświetlić wszystkie te informacje na jednym wykresie (dla każdej lokalizacji). Tak więc, dla Miami mam nadzieję, że mam jeden wykres, który pokazuje, jak temperatura różni się od średniej w różnych porach dnia i różnych porach roku (wykres 3D?)