合并并调整R中的两个重叠时间序列

Question

我有两个动物园时间序列，其中一个区域的参数相同但来自不同的平台。这两个时间序列的数字略有变化，即使它们的个别趋势是正确的，并且具有重叠的日期。我想将两个时间序列合并为一个连续序列，同时使用重叠日期的数据调整两个序列中的错误。我该怎么办？ 我在下面添加了一些示例数据。

library(ggplot2)
library(gtable)
library(grid)
library(zoo)
library(reshape)
library(reshape2)
##Create zoo objects
x<-as.zoo(as.matrix(cbind(a=1:8,b=2:9)))
y<-as.zoo(as.matrix(cbind(a=2:9,b=3:10)))
##Create dates
CCnb1<-seq(from=as.Date("2004-01-01"),to=as.Date("2004-08-01"),by="1 months")
CCnb2<-seq(from=as.Date("2004-06-01"),to=as.Date("2005-01-01"),by="1 months")
##Index appropriately
index(x)<-CCnb1
index(y)<-CCnb2
####Create dataframes
x1<-as.data.frame(x)
y1<-as.data.frame(y)
##Add date columns
x1$Date=CCnb1
y1$Date=CCnb2
##Melt data frames
x2<-melt(x1, id.vars="Date")
y2<-melt(y1, id.vars="Date")

我已经使用ggplot2包含了一个伪图，看看线条的外观。我的实际时间序列要长得多，而且价值的变化并不像这个例子那么糟糕。

#Plot
NT<-ggplot(x2, aes(x=Date, y=value,colour=variable, group=variable)) +
theme_bw()+ geom_line(size=0.5,colour="grey30")
NTb<-NT + geom_line(data=y2,aes(x=Date, y=value,group=variable))
NTb+facet_wrap(~variable)

Answer 1

差异是不变的：

x$a

所以假设y$a要调整为x$b，而y$b和va <- apply(merge(x$a + coredata(y$a - x$a)[1], y$a), 1, mean, na.rm = TRUE) a <- zoo(va, as.Date(names(va))) vb <- apply(merge(x$b + coredata(y$b - x$b)[1], y$b), 1, mean, na.rm = TRUE) b <- zoo(vb, as.Date(names(vb))) merge(a, b) 则相同：

            a  b
2004-01-01 -3 -2
2004-02-01 -2 -1
2004-03-01 -1  0
2004-04-01  0  1
2004-05-01  1  2
2004-06-01  2  3
2004-07-01  3  4
2004-08-01  4  5
2004-09-01  5  6
2004-10-01  6  7
2004-11-01  7  8
2004-12-01  8  9
2005-01-01  9 10

，并提供：

extract () {
    local pattern="$1"
    local fname="$2"
    pattern="${pattern//\?/[[:digit:]]}"
    pattern="${pattern/\*/([^[:digit:].]+)}"
    [[ $fname =~ $pattern ]]
    echo "${BASH_REMATCH[1]}"
}