รหัสต่อไปนี้ประเมินความคล้ายคลึงกันระหว่างอนุกรมเวลาสองชุด:
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
ใช้เพื่อประเมินว่าอุณหภูมิที่นิวยอร์กและไมอามีแตกต่างจากอุณหภูมิเฉลี่ย (ของทั้งสองที่) ในแต่ละช่วงเวลาของวันอย่างไร ปัญหาที่ฉันมีอยู่ในตอนนี้คือฉันต้องรวมคำที่ใช้โต้ตอบซึ่งแสดงให้เห็นว่าอุณหภูมิของแต่ละตำแหน่งแตกต่างกันอย่างไรในแต่ละวันในวันที่แตกต่างกันของปี ในที่สุดฉันก็หวังว่าจะแสดงข้อมูลทั้งหมดนี้ในกราฟเดียว (สำหรับแต่ละตำแหน่ง) ดังนั้นสำหรับไมอามี่ฉันหวังว่าจะมีกราฟหนึ่งกราฟที่แสดงว่าอุณหภูมิแตกต่างจากค่าเฉลี่ยในช่วงเวลาต่าง ๆ ของวันและเวลาที่แตกต่างกันของปี (พล็อต 3 มิติหรือไม่)