Em um gráfico genérico, com o tempo no eixo x, gostaria de destacar um período de alguns anos específicos.

Como posso fazer isso melhor? Minha ideia é, por exemplo, uma barra amarela clara para os anos destacados, por trás da trama, é claro.

O código da trama que tenho agora:

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();

Onde a temperatura $ bigEmdIm5 é a saída da decomposição do modo emperical. Os dados são em meses, então eu gostaria de destacar 01/1950 até 12/1950, por exemplo.