Суммарная частота по коэффициенту

Мне нужно узнать совокупную частоту, преобразованную в проценты, факторной непрерывной переменной. Например:

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

Мне нужна кумулятивная частота столбца G по факторам Plot ~ Site , чтобы построить график geom_step ggplot G против d для каждого участка и участка.
Мне удалось вычислить кумулятивную сумму G по коэффициенту:

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

Теперь мне нужно выразить каждый график G в диапазоне [0..1], где 1 составляет G итогов для каждого участка . Я предполагаю, что я должен разделить G на его Plot Gtot , а затем применить к нему новый cumsum . Как это сделать?
Обратите внимание, что я должен построить график этой совокупной частоты с d , а не с G как таковой, так что это неправильный ecdf.
Спасибо.