| 6996 | |
| 6997 | template <typename URng, typename Iterator, typename Estimator> |
| 6998 | sample resample(URng& rng, int resamples, Iterator first, Iterator last, Estimator& estimator) { |
| 6999 | auto n = last - first; |
| 7000 | std::uniform_int_distribution<decltype(n)> dist(0, n - 1); |
| 7001 | |
| 7002 | sample out; |
| 7003 | out.reserve(resamples); |
| 7004 | std::generate_n(std::back_inserter(out), resamples, [n, first, &estimator, &dist, &rng] { |
| 7005 | std::vector<double> resampled; |
| 7006 | resampled.reserve(n); |
| 7007 | std::generate_n(std::back_inserter(resampled), n, [first, &dist, &rng] { return first[dist(rng)]; }); |
| 7008 | return estimator(resampled.begin(), resampled.end()); |
| 7009 | }); |
| 7010 | std::sort(out.begin(), out.end()); |
| 7011 | return out; |
| 7012 | } |
| 7013 | |
| 7014 | template <typename Estimator, typename Iterator> |
| 7015 | sample jackknife(Estimator&& estimator, Iterator first, Iterator last) { |
no test coverage detected