DidzisElfertsの答えは良い答えです。時系列を扱うときはxts
パッケージを使用できることに注意してください。特に、そのプロットはベースグラフィックスで行われます。ここに例があります:
monthly.apply
私は自分のオブジェクトを分割するために使用します。
- 私はplot.xtsを使用して時系列をpltします。
まず、ランダムデータを生成します
library(xts)
days <- seq.Date( as.Date("2011-01-01"), as.Date("2011-12-31") ,1)
dat <- xts(rnorm(365),days)
## I use monthly apply to compute months widths.
## no need to give them by hand.
widths <- coredata(apply.monthly(dat,length))
par(bg="lightyellow", mar=c(2,2,2,0))
layout(matrix(c(rep(1,12),2:13),nrow=2,byrow=T),widths=widths*2)
mon <- months(days,abbreviate=T)
plot(dat,main = 'my year time series')
apply.monthly(dat,function(x) {
if(unique(format((index(x)),'%m')) =='01') {#JAN
par(mar=c(2,2,2,0)) ## special case of JAN because it contians y axis
plot(x,main='')
}
else{
par( mar=c(2,0,2,0))
plot(x,main='',ylab='')
}
})

注:JANパネルは2月1日以上です。y軸が含まれています。