-4

我有两组时间序列变量。

WSN = c(HMn25_29,HMn28_29,HMn29_29)
Ref = c(suplmentryfor29,suplmentryob29)

我想要有一个图，显示 WSN 的一个元素集和一个 Ref 元素。换句话说，我想可视化 {HMn25_29,suplmentryfor29},{HMn28_29,suplmentryfor29},{HMn28_29,suplmentryfor29},{HMn25_29,suplmentryob29},{HMn25_29,suplmentryob29},{HMn25_29,suplmentryob29}。

预期成绩：

在此处输入图像描述

WSN 的样本数据：

                    dt   temp
1 2007-09-29 00:00:00 -1.874500
2 2007-09-29 01:00:00 -1.693667
3 2007-09-29 02:00:00 -1.766667
4 2007-09-29 03:00:00 -1.928667
5 2007-09-29 04:00:00 -2.132581
6 2007-09-29 05:00:00 -1.752258

参考数据示例：

           HMn29_29.h   V1
1 2007-09-29 00:00:00 -1.8
2 2007-09-29 01:00:00 -1.9
3 2007-09-29 02:00:00 -2.4
4 2007-09-29 03:00:00 -2.3
5 2007-09-29 04:00:00 -2.3
6 2007-09-29 05:00:00 -2.4

此外，我有单独绘制的代码。单独绘制变量的代码：

plot29 <- list(HMn25_29,HMn28_29,HMn29_29,HMn31_29,HMn32_29,suplmentryfor29,suplmentryob29)
Mainn <- c('Node 25','Node 28','Node 29','Node 31','Node 32', 'MetoSwiss Forecast Grid Point', "MetoSwiss Sensing Station")
par(mfrow=c(1,2))
for (i in seq_along(plot29))  {
  plot(plot29[[i]], type = "o", pch = 16, lty = 2, col = "Black", xlab = "Hour  (2007/09/29)" , ylab = "Ambient Tempreture")
  title(Mainn[[i]])
}

现在，我需要一个自动函数来获取列表并显示具有上述属性的图，其中 x 轴表示时间-日期，y 轴表示相同比例的变量的绝对值，例如 -2 ~ 10。此外，每个情节都应包括图例。

4

1 回答 1

0

R 中的 ts.plot 和 plot.ts 函数都允许您将多条线绘制到一个图中，也许您可​​以尝试一下？这是一个简单的例子：

a<-c(1:5)
b<-c(2:6)
ts.plot(cbind(a,b))

另一种方法是在制作新图形之前使用“par(new=TRUE)”，它告诉 R 将新图形放入您上次使用的设备中，而不清除旧图形。这并不整洁，但它可能会给你更多的灵活性。如果你决定走这条路，你会想要在后面的图中省略轴和标签，并手动设置 xlim 和 ylim，否则它们会消失。

于 2013-01-11T01:45:01.020 回答