時系列をプロットするときにx軸のラベルを日付に設定する

Question

x=read.table(text=" Qtr1 Qtr2 Qtr3 Qtr4 
2010 1.8 8.0 6.0 3.0 
2011 2.0 11.0 7.0 3.5 
2012 2.5 14.0 8.0 4.2 
2013 3.0 15.2 9.5 5.0", 
    sep="",header=TRUE) 
y<-ts(as.vector(as.matrix(x)),frequency=4,start=c(2010,1)) 
plot.ts(y) 
time<-seq(as.Date("2010/1/1"),length.out=20,by="3 months") 
axis(1, at = time) 

グラフを描画するときにx軸に日付を追加したいのですが、axis(1, at = time)でx軸の日付データを追加できないのはなぜですか？

Answer 1

axis(1, at=time)を呼び出すと、ラベルにx軸をプロットするように、timeのポイントが表示されます。しかし、timeは文字のベクトルであり、数値ではありません。

一般に、axis(1, at=..., labels=...)には、実際のラベルと軸に沿って配置する場所を示します。あなたのケースでは、plot.tsへの呼び出しはx軸の制限を暗黙的に2010と2013.75に設定しているので、atのパラメータにこれらの制限が反映されるはずです。

axisを呼び出すと、ラベルはtimeで、位置は2010, 2010.25, 2010.50 ...、つまりseq(from=2010, to=2013.25, by=0.25)となります。一般的な解決方法は次のとおりです。

plot.ts(y,axes=F) # don't plot the axes yet 
axis(2) # plot the y axis 
axis(1, labels=time, at=seq(from=2010, by=0.25, length.out=length(time))) 
box() # and the box around the plot