#include <boost/accumulators/statistics/roc_curve.hpp>

namespace accumulators = boost::accumulators;

void example_roc()
{
	roc_distribution<double> roc;

	double ds=0;

	accumulators::accumulator_set<
		double
		, accumulators::stats<
			accumulators::tag::roc_curve
		>
    > acc(
		accumulators::roc_curve_num_cells = 10
		,accumulators::sample =ds
	);



	boost::lagged_fibonacci607 rng;
    boost::normal_distribution<> mean_sigma(0,1);
    boost::variate_generator<boost::lagged_fibonacci607&, boost::normal_distribution<> > normal(rng, mean_sigma);

	int n=100;
    for (std::size_t i=0; i<n; ++i)
    {
		double tmpd=normal()+1;
		double fd=1;
		roc(  tmpd,true);
		acc(tmpd,accumulators::roc_curve_positive_or_negative=fd);
    }
    for (std::size_t i=0; i<n; ++i)
    {
		double tmpd=normal()-1;
		double fd=0;
		acc(tmpd,accumulators::roc_curve_positive_or_negative=fd);   
    }
 	std::cout << accumulators::roc_curve(acc); 

}