如何绘制一条稀疏（有间隙）线，其线段根据 R 中的某些因素着色？

Question

我有一个带有时间序列的data.frame。里面还有NAs，还有一个因素我想用它来突出一行的不同部分。

flow.mndnr <- function(id, start, end) {
  uri <- sprintf("http://maps1.dnr.state.mn.us/cgi-bin/csg.pl?mode=dump_hydro_data_as_csv&site=%s&startdate=%s&enddate=%s", id, start, end)
  dat <- read.csv(url(uri), colClasses=c(Timestamp="Date"))
  rng <- range(dat$Timestamp)
  d <- data.frame(Timestamp=seq(rng[1], rng[2], by='day'))
  merge(d, dat, all.x=TRUE)
}
dat <- flow.mndnr("28062001", as.Date("2002-04-02"), as.Date("2011-10-05"))

我可以无条件地绘制它

library(lattice)
xyplot(Discharge..cfs. ~ Timestamp, dat, type='l', cex=0.5, auto.key=TRUE)

但我尝试引入因子时无法摆脱连接线

xyplot(Discharge..cfs. ~ Timestamp, dat, type='l',
    groups=dat$Discharge..cfs..Quality, cex=0.5, auto.key=TRUE)

与 ggplot2 相同

dat$quality <- dat$Discharge..cfs..Quality
ggplot(dat, aes(x=Timestamp, y=Discharge..cfs.)) +
  geom_path(aes(colour=quality)) + theme(legend.position='bottom')

我尝试了geom_line，但没有成功。我在ggplot2 mailing archive 中读到geom_path 是要走的路。但这对我来说不太管用。

附：为什么 ggplot2 不喜欢名称中的点，所以我不得不使用另一个？

Answer 1

问题在于分组。您可以使用year 跳过这些跳转。做吧：

dat$grp <- format(dat$Timestamp, "%Y")
ggplot(dat, aes(x=Timestamp, y=Discharge..cfs.)) +
    geom_path(aes(colour = quality, group = grp)) + 
    theme(legend.position='bottom')

你得到：

编辑： 详细回答评论：只要不知道要分组的变量，就不能正确分组。如果您在一年中缺少几个月，那么这段代码当然会产生跳跃。在这种情况下，我建议这样做：

dat$grp <- paste(format(dat$Timestamp, "%Y"), format(dat$Timestamp, "%m"))
ggplot(dat, aes(x=Timestamp, y=Discharge..cfs.)) +
    geom_path(aes(colour = quality, group = grp)) + 
    theme(legend.position='bottom')

你明白了：