在R中绘制时间序列的预测值和原始分数

我现在使用R中的“vars”包来检查两个时间序列之间的相互关系。具体来说，我们的数据有66个时间点。我将其分为测试样本（1-60次观察）和保留样本（61-66次观察）。我想将所有66个观测值的预测值与相同标度（从1到66）中所有66个观测值的原始分数绘制在同一个图中，以比较模型拟合。但我没有这样做的PAR和布局功能。如果您能给我一些关于它的指示，我们将不胜感激

下面是我的R代码：

library("vars")
setwd("c:$temp")  
filename<-"data.txt"
full<-read.table(filename,header=TRUE,sep="\t")
env<-full[1:60,]
varlag1<-VAR(env,p = 2,type = "const");
summary(varlag1)
plot(varlag1)
predict<-predict(varlag1,n.ahead=6,ci=0.95)
list(predict)
raw_v1<-full[1:66,1]
plot(predict,names="v1",lwd=3)
par(new=TRUE)
plot(as.ts(raw_v1),lwd=1)
raw_v2<-full[1:66,2]
plot(predict,names="v2",lwd=3)
par(new=TRUE)
plot(as.ts(raw_v2),lwd=1)

库（“vars”）
setwd（“c:$temp”）
filename它看起来不像你想要的那么漂亮，但我想你是在找这样的东西吧

欢迎来到SO。除了看起来是重复的，没有测试数据。这是应该提供的（无用的空格被编辑掉了）。有时绘图例程在内部使用
par
设置并覆盖您的工作。“未能这样做”并不是特别具有描述性。你得到什么了吗？