Eu tenho duas curvas de densidade plotadas usando isso:

Network <- Mydf$Networks
quartiles <-  quantile(Mydf$Avg.Position,  probs=c(25,50,75)/100)
density <- ggplot(Mydf, aes(x = Avg.Position, fill = Network))
d <- density + geom_density(alpha = 0.2) + xlim(1,11) + opts(title = "September 2010") + geom_vline(xintercept = quartiles, colour = "red")
print(d)

Gostaria de calcular a área sob cada curva para um determinado intervalo de posição média. Mais ou menos como pnorm para a curva normal. Alguma ideia?