На общем графике со временем по оси X я хотел бы выделить период в несколько конкретных лет.

Как я могу лучше всего это сделать? Моя идея, например, светло-желтая полоса для выделенных лет, конечно же, за сюжетом.

Сюжетный код у меня сейчас:

pdf("temperature_imfs_big_interm5.pdf", width=6, height=8);
par(mfrow=c(temperature$bigEmdIm5$nimf+1,1), mar=c(2,1,2,1))
for(i in 1:temperature$bigEmdIm5$nimf) {
    plot(timeline$big, temperature$bigEmdIm5$imf[,i], type="l", xlab="", ylab="", ylim=range(temperature$bigEmdIm5$imf[,i]), axes=FALSE, main=paste(i, "-th IMF", sep=""))#; abline(h=0)
  axis.POSIXct(side=1, at=tickpos$big)
}
plot(timeline$big, temperature$bigEmdIm5$residue, xlab="", ylab="", axes=FALSE, main="residue", type="l")
axis.POSIXct(side=1, at=tickpos$big)
dev.off();

Где температура $ bigEmdIm5 является выходом эмпирического режима разложения. Данные даны в месяцах, поэтому я хотел бы выделить 01/1950 до 12/1950, например.