ggplot2: grafik paneli kümesinde yer paylaşımlı denetim grubu satırı

Question

ggplot2 ile yapılmış yığılmış bir alan grafiğim var:

dists.med.areaplot<-qplot(starttime,value,fill=dists,facets=~groupname, 
    geom='area',data=MDist.median, stat='identity') + 
    labs(y='median distances', x='time(s)', fill='Distance Types')+ 
    opts(title=subt) + 
    scale_fill_brewer(type='seq') + 
    facet_wrap(~groupname, ncol=2) + grect #grect adds the grey/white vertical bars 

Şunlara benziyor:

Tüm grafiklere kontrol grafiğinin (sağ alt) profilinin bir katmanı eklemek istiyorum çıktıda (groupname == rowH kontroldür).

Şimdiye kadar elimden gelen en iyi çabalarım:

cline<-geom_line(aes(x=starttime,y=value), 
    data=subset(dists.med,groupname=='rowH'),colour='red') 

dists.med.areaplot + cline 

Koyu mavi bölümün üstünü kaplayan 1 kırmızı çizgi olmak için 3 kırmızı çizgiye ihtiyacım var. Ve her bir panelin üstünü örtmek için o aynı çizgiye (rowH hattı) ihtiyacım var.

Veri yapısı şuna benzer:

> str(MDist.median) 
'data.frame': 2880 obs. of 6 variables: 
$ groupname: Factor w/ 8 levels "rowA","rowB",..: 1 1 1 1 1 1 1 1 1 1 ... 
$ fCycle : Factor w/ 6 levels "predark","Cycle 1",..: 1 1 1 1 1 1 1 1 1 1 ... 
$ fPhase : Factor w/ 2 levels "Light","Dark": 2 2 2 2 2 2 2 2 2 2 ... 
$ starttime: num 0.3 60 120 180 240 300 360 420 480 540 ... 
$ dists : Factor w/ 3 levels "inadist","smldist",..: 1 1 1 1 1 1 1 1 1 1 ... 
$ value : num 110 117 115 113 114 ... 

Kırmızı çizgi her starttime value toplamı olarak hesaplanmalıdır, burada groupname = 'rowH'. Aşağıdaki yolları cline oluşturmaya çalıştım. Her bir hata veya yanlış bir sonuç ortaya çıkar:

#sums the entire y for all points and makes horizontal line 
cline<-geom_line(aes(x=starttime,y=sum(value)),data=subset(dists.med,groupname=='rowH'),colour='red') 

#using related dataset with pre-summed y's 
> cline<-geom_line(aes(x=starttime,y=tot_dist),data=subset(t.med,groupname=='rowH')) 
> dists.med.areaplot + cline 
Error in eval(expr, envir, enclos) : object 'dists' not found 

Düşünceler?

ETA:

'dists' not found ile karşılaştığım sorunun, ilk arsa dists.med.areaplot'ün qplot ile yaratılmış olması gerektiği anlaşılıyor. Bu sorunu önlemek için bir qplot üzerinde yapılamam. Bu, çalışma grafiğinin kodudur:

cline.data <- subset(
     ddply(MDist.median, .(starttime, groupname), summarize, value = sum(value)), 
     groupname == "rowH") 
cline<-geom_line(data=transform(cline.data,groupname=NULL), colour='red') 

dists.med.areaplot<-ggplot(MDist.median, aes(starttime, value)) + 
    grect + nogrid + 
    geom_area(aes(fill=dists),stat='identity') + 
    facet_grid(~groupname)+ scale_fill_brewer(type='seq') + 
    facet_wrap(~groupname, ncol=2) + 
    cline 

bu grafikle sonuçlanır:

Answer 1

Bu Öğrenim R blog yazısı bazı yardımcı olmalıdır:

http://learnr.wordpress.com/2009/12/03/ggplot2-overplotting-in-a-faceted-scatterplot/

Olabilir Bu çalışma arsa kodudur ile ggplot dışında özet hesaplamaya değer.

cline.data <- ddply(MDist.median, .(starttime, groupname), summarize, value = sum(value)) 
cline.data.subset <- subset(cline.data, groupname == "rowH") 

Sonra

last_plot() + geom_line(data = transform(cline.data.subset, groupname = NULL), color = "red")