以下代码评估两个时间序列之间的相似性:
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")
在这里,gam
用于评估纽约和迈阿密的温度与一天中不同时间(两个地点)的平均温度之间的差异。我现在遇到的问题是,我需要包括一个交互项,以显示每个位置的温度在一年中不同日期的一天中如何变化。我最终希望将所有这些信息显示在一张图上(针对每个位置)。因此,对于迈阿密,我希望有一张图表可以显示在一天中的不同时间和一年中的不同时间,温度与平均值的变化情况(3d图?)
2
您可能会发现此问题的答案stats.stackexchange.com/questions/18937/…相关。
—
jbowman