我正在使用dlnm
和lm
来拟合多个时间序列的线性模型。我想在视觉上将拟合值与实际值进行比较,因此我想在一个csv文件中导出它们。但是,我似乎无法使用ts.intersect
绑定它们,因为模型的$fitted.values
不是时间序列。如果将lm
传递给$fitted.values
,则na.action = NULL
suggests的文档将lm
为时间序列。但是,只有在用作输入的值被定义到任何地方时才有可能实现这种情况,但如果crossbasis
应用了非平凡滞后则不是这种情况。
有没有一种方便的方法可以并排获得实际和合适的时间序列?
答案 0 :(得分:0)
我找到了一个合理的解决方法,将fit.values重建为时间序列。 ts.intersect
不适用于zoo
系列,merge
是适当的功能:
actual <- actual.values
fit <- read.zoo(data.frame(names(model$fitted.values),model$fitted.values))
result <- merge(actual,fit,all=FALSE) #default all=TRUE gives union