我每隔一小时测量两个固定位置的气温测量值。下面的代码是一组随机数字,用于表示我的数据格式:
set.seed(1)
RandData <- rnorm(8760*2,sd=10)
Locations <- rep(c('UK','France'),each=8760)
Date = seq(from=as.POSIXct("1991-01-01 00:00"),
to=as.POSIXct("1991-12-31 23:00"), length=8760)
Final <- data.frame(Loc = Locations,
Doy = as.numeric(format(Date,format = "%j")),
Tod = as.numeric(format(Date,format = "%H")),
Temp = RandData)
我可以使用以下代码将温度变化绘制为一年中的一天的函数:
require(lattice)
xyplot(Temp~Doy | Loc, data = Final, col = "black", type = "l")
这将显示数据的年度模式。但是,我想做的是制作一天中不同时间温度变化的箱线图。所以,对于上面的例子,我想要两个数字,每个国家一个,每个数字应该由显示 00:00、01:00... 的温度变化的箱形图组成,以此类推,参考 Final$Tod . 如何做到这一点?
非常感谢您的帮助。