我有不同组组件的时间序列数据.具有各种时间戳(以日期给出)的每个组ID具有低和超响应数据.我想通过facet(ggplot)绘制每个组的时间序列,用于(1)组ID和响应,即(2)Hyper和Hypo响应,以便响应的图像是另一个的顶部.任何帮助表示赞赏.
下面给出了一个演示数据集和我到目前为止所做的工作.
set.seed(1)
tdat <- data.frame(Group = rep(paste0("GroupID-", c("A","B")),
each = 100),
Date = rep(seq(Sys.Date(), by = "1 day", length = 100), 2),
Fitted = c(cumsum(rnorm(100)), cumsum(rnorm(100))),
Signif = rep(NA, 200))
tdat <- transform(tdat, Hyper = Fitted + 1.5, Hypo = Fitted - 1.5)
## select 1 region per Site as signif
take <- sample(10:70, 2)
take[2] <- take[2] + 100
tdat$Signif[take[1]:(take[1]+25)] <- tdat$Fitted[take[1]:(take[1]+25)]
tdat$Signif[take[2]:(take[2]+25)] <- tdat$Fitted[take[2]:(take[2]+25)]
数据框看起来像这样 -
> head(tdat)
Group Date Fitted Signif Hyper Hypo
1 GroupID-A 2017-04-18 …