-4

時系列変数のセットが 2 つあります。

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

Ref の要素で設定された WSN の要素を表すプロットが必要でした。つまり、{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]])
}

ここで、リストを取得し、前述のプロパティを含むプロットを 2 つずつ表示する自動関数が必要です。ここで、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 に答える