Because the result is a list, you can get your number like this: integrate(dnorm,0,2)[[1]] You don't have an expression for the probability function? You can do a nonlinear regression to a distribution you belive in, and then integrate that. Or, you can calculate the experimental cdf by simple summation. Best, Adrian