Jak ustawić pionowo dwa wykresy o tej samej skali x, ale innej skali Y w R?

Pozdrowienia,

Obecnie wykonuję następujące czynności w języku R:

require(zoo)
data <- read.csv(file="summary.csv",sep=",",head=TRUE)
cum  = zoo(data$dcomp, as.Date(data$date))
data = zoo(data$compressed, as.Date(data$date))
data <- aggregate(data, identity, tail, 1)
cum  <- aggregate(cum, identity, sum, 1)
days = seq(start(data), end(data), "day")
data2 = na.locf(merge(data, zoo(,days)))

plot(data2,xlab='',ylab='compressed bytes',col=rgb(0.18,0.34,0.55))
lines(cum,type="h",col=rgb(0,0.5,0))

Snip of summary.csv:

date,revision,file,lines,nclass,nattr,nrel,bytes,compressed,diff,dcomp
2007-07-25,16,model.xml,96,11,22,5,4035,991,0,0
2007-07-27,17,model.xml,115,16,26,6,4740,1056,53,777
2007-08-09,18,model.xml,106,16,26,7,4966,1136,47,761
2007-08-10,19,model.xml,106,16,26,7,4968,1150,4,202
2007-09-06,81,model.xml,111,16,26,7,5110,1167,13,258
...

Ostatnie dwa wiersze przedstawiają potrzebne informacje, a wynik przypomina: Niebieska linia to entropia w bajtach artefaktu, który mnie interesuje. Zielone linie reprezentują entropię zmian.

Teraz na tym wykresie działa dobrze, ponieważ nie ma ogromnej różnicy w skalach. Ale mam inne wykresy, na których zielone linie stają się tak małe, że nie widać.

Rozwiązanie, którego szukałem, obejmowało dwie rzeczy:

1. Aby przenieść zielone pionowe linie na drugi wykres, tuż poniżej pierwszego, z własną osią y, ale wspólną osią x.
2. Aby zapewnić skalę logarytmiczną, ponieważ bardziej interesuje mnie „wielkość” niż określone wartości.

Z góry dziękuję!

PS Jeśli ktoś może mi również powiedzieć, jak mogę umieścić „drobne tiki” w skali X w odniesieniu do miesięcy, doceniam :-) Jeśli są to zbyt wiele pytań do jednego postu, mogę je dalej podzielić.

Możesz użyć par(new=TRUE)do wykreślenia na tym samym wykresie przy użyciu dwóch różnych osi Y! To powinno również rozwiązać twój problem.

Następnie znajdziesz prosty przykład, który wykreśla dwie losowe zmienne normalne, jedną na średniej 0, drugą na średniej 100 (oba sd s = 1) na tym samym wykresie. Pierwszy na czerwono na lewej osi y, drugi na niebiesko na prawej osi y. Następnie dodaje się etykiety osi.

Proszę bardzo:

x <- 1:10
y1 <- rnorm(10)
y2 <- rnorm(10)+100

plot(x,y1,pch=0,type="b",col="red",yaxt="n",ylim=c(-8,2))
par(new=TRUE)
plot(x,y2,pch=1,type="b",col="blue",yaxt="n",ylim=c(98,105))

axis(side=2)
axis(side=4)

wygląda wtedy tak (pamiętaj czerwony na lewej osi, niebieski na prawej osi):

AKTUALIZACJA:
Na podstawie komentarzy stworzyłem zaktualizowaną wersję mojego wykresu. Teraz zagłębiam się nieco w funkcjonalność wykresu bazowego, używając par(mar=c(a,b,c,d))do stworzenia większego marginesu wokół wykresu (potrzebnego do etykiety prawej osi), mtextaby wyświetlić etykiety osi i zaawansowane użycie axisfunkcji:

x <- 1:100
y1 <- rnorm(100)
y2 <- rnorm(100)+100

par(mar=c(5,5,5,5))

plot(x,y1,pch=0,type="b",col="red",yaxt="n",ylim=c(-8,2),ylab="")
axis(side=2, at=c(-2,0,2))
mtext("red line", side = 2, line=2.5, at=0)

par(new=TRUE)
plot(x,y2,pch=1,type="b",col="blue",yaxt="n",ylim=c(98,108), ylab="")
axis(side=4, at=c(98,100,102), labels=c("98%","100%","102%"))
mtext("blue line", side=4, line=2.5, at=100)

Jak widać, jest to całkiem proste. Możesz zdefiniować pozycję swoich danych ylimw plotfunkcji, a następnie użyć atw axisfunkcji, aby wybrać, które oś znacznika chcesz zobaczyć. Co więcej, można nawet dostarczyć etykiety dla kleszczy osi (bardzo przydatnych dla nominalnej oś x) poprzez labelsw axisfunkcji (wykonanej tutaj na prawej osi). Aby dodać etykiety osi, użyj mtextz atdo pozycjonowania pionowego ( linedo pozycjonowania poziomego).

Upewnij się, aby sprawdzić ?plot, ?par, ?axis, i ?mtextaby uzyskać więcej informacji.
Świetne zasoby internetowe to: Quick-R dla wykresów: 1 , 2 i 3 .

Myślę, że możesz dostać to, czego chcesz ggplot2. Korzystając z poniższego kodu, mogę utworzyć:

Oczywiście rzeczy takie jak kolory linii można zmienić na dowolne. Na osi X określiłem linie główne w latach i linie pomocnicze w miesiącach.

require(ggplot2)
t = as.Date(0:1000, origin="2008-01-01")  
y1 = rexp(1001)
y2 = cumsum(y1)
df = data.frame(t=t, values=c(y2,y1), type=rep(c("Bytes", "Changes"), each=1001))

g = ggplot(data=df, aes(x=t, y=values)) +
  geom_line() +
  facet_grid(type ~ ., scales="free") +
  scale_y_continuous(trans="log10") +
  scale_x_date(major="years", minor="months") +
  ylab("Log values")
g