ggplot2: vùng biểu đồ đánh dấu

Question

Tôi đang làm việc với dữ liệu chuỗi thời gian và muốn đánh dấu khu vực biểu đồ bất cứ khi nào một số điều kiện nhất định trở thành sự thật. Ví dụ:

require(ggplot2) 
require(quantmod) 
initDate <- "1993-01-31" 
endDate <- "2012-08-10" 
symbols <- c("SPY") 
getSymbols(symbols, from=initDate, to=endDate, index.class=c("POSIXt","POSIXct")) 
spy<-SPY$SPY.Adjusted 
spy$sma<-SMA(spy$SPY.Adjusted,200) 
spy<-spy[-(1:199),] 
spy<-as.data.frame(spy) 
ggplot(spy,aes(x=index(spy),y=spy$SPY.Adjusted))+geom_line()+geom_line(aes(x=index(spy),y=spy$sma)) 

Đoạn mã trên vẽ lô dữ liệu, nhưng làm cách nào để làm nổi bật phần này khi nào gần trên sma? Câu hỏi này tương tự như How to highlight time ranges on a plot?, nhưng sau đó nó là thủ công. Có chức năng nào trong ggplot2 cho âm mưu có điều kiện không?

Answer 1

Dựa trên mã trong tệp TA.R của gói quantmod, đây là mã sử dụng rle để tìm bắt đầu và kết thúc của hình chữ nhật.

runs <- rle(as.logical(spy[, 1] > spy[, 2])) 
l <- list(start=cumsum(runs$length)[which(runs$values)] - runs$length[which(runs$values)] + 1, 
      end=cumsum(runs$lengths)[which(runs$values)]) 
rect <- data.frame(xmin=l$start, xmax=l$end, ymin=-Inf, ymax=Inf) 

Cùng với đó là một số ggplot2 mã từ các accepted answer cho câu hỏi bạn có liên quan đến:

ggplot(spy,aes(x=index(spy),y=spy$SPY.Adjusted))+geom_line()+geom_line(aes(x=index(spy),y=spy$sma))+geom_rect(data=rect, aes(xmin=xmin, xmax=xmax, ymin=ymin, ymax=ymax), color="grey20", alpha=0.5, inherit.aes = FALSE) 

Và bạn nhận được:

Nếu bạn đảo ngược thứ tự của âm mưu và sử dụng alpha=1 trong geom_rect nó có thể (hoặc có thể không) trông giống như bạn mong muốn:

ggplot(spy,aes(x=index(spy),y=spy$SPY.Adjusted))+geom_rect(data=rect, aes(xmin=xmin, xmax=xmax, ymin=ymin, ymax=ymax), border=NA, color="grey20", alpha=1, inherit.aes = FALSE)+geom_line()+geom_line(aes(x=index(spy),y=spy$sma)) 

---

Vì bạn có một đối tượng xts. Thậm chí bạn có thể không muốn chuyển đổi thành một số data.frame. Đây là cách bạn có thể vẽ nó sử dụng thương hiệu plot.xts phương pháp mới trong gói xtsExtra tạo ra bởi Michael Weylandt như một phần của một mùa hè Google của Bộ luật project.

spy <- as.xts(spy) 
require(xtsExtra) 
plot(spy, screens=1, 
    blocks=list(start.time=paste(index(spy)[l$start]), 
       end.time=paste(index(spy)[l$end]), col='lightblue'),      
    legend.loc='bottomright', auto.legend=TRUE)