O código a seguir avalia a semelhança entre duas séries temporais:
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")
Aqui, gam
é usado para avaliar como a temperatura em Nova York e Miami varia da temperatura média (de ambos os locais) em diferentes momentos do dia. O problema que tenho agora é que preciso incluir um termo de interação que mostre como a temperatura de cada local varia ao longo do dia em diferentes dias do ano. Eventualmente, espero exibir todas essas informações em um gráfico (para cada local). Então, para Miami, espero ter um gráfico que mostre como a temperatura varia em relação à média em diferentes horários do dia e em diferentes épocas do ano (gráfico 3D?)