Як вертикально скласти два графіки з однаковою шкалою x, але різною шкалою y в R?

Привітання,

В даний час я роблю наступне в 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))

Фрагмент резюме.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
...

Останні два рядки описують потрібну мені інформацію, і результат нагадує наступне: Синя лінія - це ентропія в байтах артефакту, який мене цікавить. Зелені лінії відображають ентропію змін.

Тепер у цьому графіку він працює добре, оскільки масштабів немає великої різниці. Але в мене є інші графіки, де зелені лінії стають такими маленькими, що їх неможливо побачити.

Я шукав рішення, включаючи дві речі:

1. Щоб перемістити зелені вертикальні лінії до другого графіка, трохи нижче першого, із власною віссю y, але розділеною осі x.
2. Щоб надати йому логарифмічну шкалу, оскільки мене більше цікавить "величина", ніж конкретні значення.

Спасибі заздалегідь!

PS Якщо хтось також може сказати мені, як я міг би поставити "незначні кліщі" у шкалі x, що стосуються місяців, я вдячний :-) Якщо це занадто багато питань для однієї посади, я можу поділити їх далі.

Ви можете використовувати par(new=TRUE)для побудови на одному графіку за допомогою двох різних осей y! Це також повинно вирішити вашу проблему.

Далі ви знайдете простий приклад, який побудує дві випадкові нормальні змінні, одна на середній 0, а інша на середній 100 (обидва sd s = 1) в одному графіку. Перший червоним кольором на лівій осі y, другий синім кольором на правій осі у. Потім додаються мітки осі.

Ось вам:

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)

виглядає ось так (згадайте червоний на лівій осі, синій на правій осі):

ОНОВЛЕННЯ:
На основі коментарів я створив оновлену версію свого графіка. Тепер я трохи більше заглиблююсь у функціональність базового графіка, використовуючи par(mar=c(a,b,c,d))для створення більшого поля навколо графіка (необхідного для правої мітки осі), mtextдля показу міток осі та розширеного використання axisфункції:

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)

Як бачите, це досить прямо вперед. Ви можете визначити позицію своїх даних ylimу plotфункції, а потім atза допомогою axisфункції вибрати, яку клітинку осі ви хочете бачити. Крім того, ви можете навіть надати мітки для осі кліщів (досить корисно для номінальної осі x) за labelsдопомогою axisфункції (зроблено тут на правій осі). Щоб додати мітки осей, використовувати mtextз atдля вертикального позиціонування ( lineдля горизонтального позиціонування).

Переконайтеся в тому , щоб перевірити ?plot, ?par, ?axisі ?mtextдля отримання додаткової інформації.
Великими веб-ресурсами є: Швидкий R для графіків: 1 , 2 та 3 .

Я думаю, ви можете отримати те, що хочете, використовуючи ggplot2. Використовуючи наведений нижче код, я можу створити:

Очевидно, що такі речі, як кольори ліній, можна змінити на все, що завгодно. На осі x я вказав основні лінії по роках і другорядні лінії по місяцях.

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