Frecuencia acumulada por factor

Question

Tengo que averiguar la frecuencia acumulada, convertida en porcentaje, de una variable continua por factor. Por ejemplo:

data <- data.frame(n = sample(1:12), 
       d = seq(10, 120, by = 10), 
       Site = rep(c("FirstSite", "SecondSite"), 6), 
       Plot = rep(c("Plot1", "Plot1", "Plot2", "Plot2"), 3) 
       ) 

data <- with(data, data[order(Site,Plot),]) 
data <- transform(data, G = ((pi * (d/2)^2) * n)/10000) 

data 
    n d  Site Plot   G 
1 7 10 FirstSite Plot1 0.05497787 
5 9 50 FirstSite Plot1 1.76714587 
9 12 90 FirstSite Plot1 7.63407015 
3 10 30 FirstSite Plot2 0.70685835 
7 5 70 FirstSite Plot2 1.92422550 
11 1 110 FirstSite Plot2 0.95033178 
2 3 20 SecondSite Plot1 0.09424778 
6 8 60 SecondSite Plot1 2.26194671 
10 6 100 SecondSite Plot1 4.71238898 
4 4 40 SecondSite Plot2 0.50265482 
8 2 80 SecondSite Plot2 1.00530965 
12 11 120 SecondSite Plot2 12.44070691 

necesito la frecuencia cumulaive de la columna G por factores Plot~Site con el fin de trazar una ggplot geom_step de G contra d para cada parcela y el sitio.
que he logrado para calcular la suma acumulada de G por el factor por:

data.ss <- by(data[, "G"], data[,c("Plot", "Site")], function(x) cumsum(x)) 
# Gtot 
(data.ss.tot <- sapply(ss, max)) 
[1] 9.456194 3.581416 7.068583 13.948671 

ahora tengo que expresar cada PlotG en el rango [0..1] donde 1 es G tot para cada Plot. Imagino que debería dividir G por su PlotGtot, luego aplicarle un nuevo cumsum. ¿Cómo hacerlo?
Tenga en cuenta que tengo que trazar esta frecuencia acumulada contra d no G sí mismo, por lo que no es un ecdf apropiado.
Gracias.

Answer 1

lo general el uso ddply y transform para hacer este tipo de cosas:

> data = ddply(data, c('Site', 'Plot'), transform, Gsum=cumsum(G), Gtot=sum(G)) 
> qplot(x=d, y=Gsum/Gtot, facets=Plot~Site, geom='step', data=data)