Compute statistical uncertainty
Currently the ROC function returns the numerical uncertainty of the area under the curve introduced by the integration. A more physical approach would be to take the statistical uncertainty into account introduce by MC stats and k-fold cross validation.