Выделитесь (заштриховывают) фон графика в определенном диапазоне времени
На универсальном графике, со временем на оси X, я хотел бы выделить период нескольких определенных лет.
Как может я bestly делать это? Моя идея является, например, светло-желтой панелью в течение выделенных лет позади графика, конечно.
Код графика я имею теперь:
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();
Где temperature$bigEmdIm5 является выводом эмпирического разложения режима. Данные находятся в течение многих месяцев, таким образом, я хотел бы выделить 01/1950 до 12/1950, например.
3 ответа
Использование альфа-прозрачности:
x <- seq(as.POSIXct("1949-01-01", tz="GMT"), length=36, by="months")
y <- rnorm(length(x))
plot(x, y, type="l", xaxt="n")
rect(xleft=as.POSIXct("1950-01-01", tz="GMT"),
xright=as.POSIXct("1950-12-01", tz="GMT"),
ybottom=-4, ytop=4, col="#123456A0") # use alpha value in col
idx <- seq(1, length(x), by=6)
axis(side=1, at=x[idx], labels=format(x[idx], "%Y-%m"))
или нанесение выделенной области за линиями:
plot(x, y, type="n", xaxt="n")
rect(xleft=as.POSIXct("1950-01-01", tz="GMT"),
xright=as.POSIXct("1950-12-01", tz="GMT"),
ybottom=-4, ytop=4, col="lightblue")
lines(x, y)
idx <- seq(1, length(x), by=6)
axis(side=1, at=x[idx], labels=format(x[idx], "%Y-%m"))
box()
Вот решение, в котором используется zoo просто потому, что это упрощает подмножество. То же самое можно сделать и со стандартной индексацией:
## create a long monthly sequence and a sub-sequence
months <- seq( as.Date("1950-01-01"), as.Date("2009-12-12"), by="month")
subset <- seq( as.Date("1970-01-01"), as.Date("1979-12-31"), by="month")
## generate some random values
set.seed(42)
values <- cumsum(rnorm(length(months)))
## plot as a zoo object, overlay a gray background and overplot a line in red
library(zoo)
Z <- zoo(values, months)
plot(Z)
rect(xleft=head(subset,1), xright=tail(subset,1),
ybottom=par("usr")[3], ytop=par("usr")[4],
density=NA, col="lightgray")
lines(Z[subset], col='red')
box()
(источник: eddelbuettel.com )
Используя par ("usr") , мы избегаем необходимости явного значения для верхней и нижней отметок области. А индексирование zoo упрощает поиск начальной и конечной точек. Это будет работать одинаково для данных с разным временным разрешением.
You can use the chartSeries() function in quantmod with an xts timeSeries and the addTA() function to add the background highlighting:
addTA(xts(rep(TRUE,length(times)), times), on=-1, col="#333333", border=NA)