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

namespace accumulators = boost::accumulators;

//using namespace boost::time_series;

void example_roc()
{
	roc_distribution<double> roc;

	double ds=0;

	accumulators::accumulator_set<
		double
		, accumulators::stats<
			accumulators::tag::roc_curve
			<double , accumulators::tag::covariate1>
		>
    > acc(
		accumulators::roc_curve_num_cells = 100
		,accumulators::sample =ds
	//, accumulators::covariate1 = 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;
		roc(  tmpd,true);
    }
    for (std::size_t i=0; i<n; ++i)
    {
		double tmpd=normal()-1;
		roc(  tmpd,false);        
    }
	roc.get_result();
	roc.roc_data;

}