Jak uzyskać cieniowane tło w Xyplot w R?

Question

przy użyciu xyplot z pakietu lattice, planuję szereg czasowy na przestrzeni wielu lat. Chciałbym dodać zacieniony obszar na niektóre z tych lat, aby wskazać, że ten okres był "wyjątkowy" (np. Wojna).

Proszę, przepraszam, jeśli to trywialne, ale nie mogłem wymyślić, jak to zrobić, więc byłbym szczęśliwy, gdyby ktoś mógł mi pomóc, lub przynajmniej wskazać mi właściwy kierunek. Myślę, że moim głównym problemem jest to, że tak naprawdę nie wiem, jak podejść do tego problemu. Nadal jestem względnie nowy dla R, a zwłaszcza dla lattice.

tu minimalny przykład:

xyplot(rnorm(100) ~ 1:100, type="l", col="black") 

W tym samym wykresie, ja jak kolor tła (od powiedzmy x -values ​​o 45, aż do 65) od dołu do góry kreślenia być zacienionym, powiedzmy, jasnoszarym.

Należy zauważyć, że rozwiązania, które znalazłem do tej pory korzystać base grafiki i polygon -function, ale intencją jest cieniu obszar pod lub nad krzywą, która różni się od tego, co chciałbym robić. Nie chcę "tylko" zacieniać obszaru poniżej linii lub powyżej linii. Zamiast tego chciałbym zacienić całe tło dla danego przedziału czasowego.

Jeśli ktokolwiek mógłby mi pomóc, byłbym bardzo wdzięczny!

Answer 1

Zobacz ?panel.xblocks w pakiecie latticeExtra:

library(latticeExtra) 
x <- 1:100 
xyplot(rnorm(100) ~ x, type="l", col="black") + 
     layer_(panel.xblocks(x, x > 20, col = "lightgrey")) 

Answer 2

Spróbuj:

xyplot( 
    rnorm(100) ~ 1:100, type="l", col="black", 
    panel=function (x,y,...){ 
    panel.rect(xleft=45, xright=65,ybottom=-3, ytop=3,col="grey") 
    panel.xyplot(x,y,...) 
    } 
) 

Funkcja panel.rect() steruje prostokąta i jest odpowiednikiem kraty funkcji rect(). Ma wiele ustawień, które mogą okazać się przydatne. Nazywa się go najpierw, a następnie umieszcza się na nim xyplot(). Wiele osób musi się trochę rozegrać, aby parametry wyglądały tak, jak je lubisz.

Answer 3

trellis.focus("panel", 1, 1) 
grid.rect(x =.55, , y=.5, w = .2, height=6, 
    gp = gpar(fill = "light grey")) 
trellis.unfocus() 

To różni się od rozwiązania @ JohnPaul na kilka sposobów (i myślę, że jego odpowiedź jest lepsza). Wykorzystuje to środek żądanego pasma X do umieszczenia w "natywnych współrzędnych" i oblicza szerokość jako "zakres (xlim)/zakres (pasmo)" i modyfikuje istniejący wykres. funkcja grid.rect to funkcja niższego poziomu pakietów grid używana przez panel.rect. Czasem wydaje mi się to przydatne, gdy integruję panele kratowe wewnątrz systemu xyplot pokonując mnie.