Расчет площади под непрерывным графиком плотности
У меня есть две кривые плотности, построенные с использованием этого:
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)
Я хотел бы вычислить площадь под каждой кривой для данного диапазона Avg.Position. Вроде как pnorm для нормальной кривой. Есть идеи?