我想在同一个图中绘制一年中不同季节(分为冷季或暖季)的时间序列数据,并确定这两个时期。我怎样才能做到这一点?
可重复的数据如下:
library(gamair)
data(chicago)
chicago$date<-seq(from=as.Date("1987-01-01"), to=as.Date("2000-12-31"),length=5114)
data.cold <- subset(chicago, quarters(date) %in% c("Q1", "Q4"))
data.warm <- subset(chicago, quarters(date) %in% c("Q2", "Q3"))
下面的这些代码将为两个时期创建两个单独的图:
with(data.cold,plot(date,death,pch=".", ylab= expression("Mortality Count"), main = "Daily Mortality in cold season"))
with(data.warm ,plot(date,death,pch=".", ylab= expression("Mortality Count"), main = "Daily Mortality in warm season"))
下面的代码创建了一个图,但是这两个时期没有明确的界限。我的愿望是在一个情节中为两个时期创建可识别的情节。
with(chicago,plot(date,death,pch=".", ylab= expression("Mortality Count"), main = "Daily Mortality from 1987-2000"))