/usr/include/boost/accumulators/statistics
Edit: /usr/include/boost/accumulators/statistics/pot_tail_mean.hpp (7802B)
///////////////////////////////////////////////////////////////////////////////
// pot_tail_mean.hpp
//
// Copyright 2006 Daniel Egloff, Olivier Gygi. Distributed under the Boost
// Software License, Version 1.0. (See accompanying file
// LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
#ifndef BOOST_ACCUMULATORS_STATISTICS_POT_TAIL_MEAN_HPP_DE_01_01_2006
#define BOOST_ACCUMULATORS_STATISTICS_POT_TAIL_MEAN_HPP_DE_01_01_2006
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
namespace boost { namespace accumulators
{
namespace impl
{
///////////////////////////////////////////////////////////////////////////////
// pot_tail_mean_impl
//
/**
@brief Estimation of the (coherent) tail mean based on the peaks over threshold method (for both left and right tails)
Computes an estimate for the (coherent) tail mean
\f[
\widehat{CTM}_{\alpha} = \hat{q}_{\alpha} - \frac{\bar{\beta}}{\xi-1}(1-\alpha)^{-\xi},
\f]
where \f$\bar[u]\f$, \f$\bar{\beta}\f$ and \f$\xi\f$ are the parameters of the
generalized Pareto distribution that approximates the right tail of the distribution (or the
mirrored left tail, in case the left tail is used). In the latter case, the result is mirrored
back, yielding the correct result.
*/
template
struct pot_tail_mean_impl
: accumulator_base
{
typedef typename numeric::functional::fdiv::result_type float_type;
// for boost::result_of
typedef float_type result_type;
pot_tail_mean_impl(dont_care)
: sign_((is_same::value) ? -1 : 1)
{
}
template
result_type result(Args const &args) const
{
typedef
typename mpl::if_<
is_same
, tag::weighted_peaks_over_threshold
, tag::peaks_over_threshold
>::type
peaks_over_threshold_tag;
typedef
typename mpl::if_<
is_same
, tag::weighted_pot_quantile
, tag::pot_quantile
>::type
pot_quantile_tag;
extractor const some_peaks_over_threshold = {};
extractor const some_pot_quantile = {};
float_type beta_bar = some_peaks_over_threshold(args).template get<1>();
float_type xi_hat = some_peaks_over_threshold(args).template get<2>();
return some_pot_quantile(args) - this->sign_ * beta_bar/( xi_hat - 1. ) * std::pow(
is_same::value ? args[quantile_probability] : 1. - args[quantile_probability]
, -xi_hat);
}
// make this accumulator serializeable
template
void serialize(Archive & ar, const unsigned int file_version)
{
ar & sign_;
}
private:
short sign_; // if the fit parameters from the mirrored left tail extreme values are used, mirror back the result
};
} // namespace impl
///////////////////////////////////////////////////////////////////////////////
// tag::pot_tail_mean
// tag::pot_tail_mean_prob
//
namespace tag
{
template
struct pot_tail_mean
: depends_on, pot_quantile >
{
/// INTERNAL ONLY
///
typedef accumulators::impl::pot_tail_mean_impl impl;
};
template
struct pot_tail_mean_prob
: depends_on, pot_quantile_prob >
{
/// INTERNAL ONLY
///
typedef accumulators::impl::pot_tail_mean_impl impl;
};
template
struct weighted_pot_tail_mean
: depends_on, weighted_pot_quantile >
{
/// INTERNAL ONLY
///
typedef accumulators::impl::pot_tail_mean_impl impl;
};
template
struct weighted_pot_tail_mean_prob
: depends_on, weighted_pot_quantile_prob >
{
/// INTERNAL ONLY
///
typedef accumulators::impl::pot_tail_mean_impl impl;
};
}
// pot_tail_mean(with_threshold_value) -> pot_tail_mean
template
struct as_feature(with_threshold_value)>
{
typedef tag::pot_tail_mean type;
};
// pot_tail_mean(with_threshold_probability) -> pot_tail_mean_prob
template
struct as_feature(with_threshold_probability)>
{
typedef tag::pot_tail_mean_prob type;
};
// weighted_pot_tail_mean(with_threshold_value) -> weighted_pot_tail_mean
template
struct as_feature(with_threshold_value)>
{
typedef tag::weighted_pot_tail_mean type;
};
// weighted_pot_tail_mean(with_threshold_probability) -> weighted_pot_tail_mean_prob
template
struct as_feature(with_threshold_probability)>
{
typedef tag::weighted_pot_tail_mean_prob type;
};
// for the purposes of feature-based dependency resolution,
// pot_tail_mean and pot_tail_mean_prob provide
// the same feature as tail_mean
template
struct feature_of >
: feature_of
{
};
template
struct feature_of >
: feature_of
{
};
// So that pot_tail_mean can be automatically substituted
// with weighted_pot_tail_mean when the weight parameter is non-void.
template
struct as_weighted_feature >
{
typedef tag::weighted_pot_tail_mean type;
};
template
struct feature_of >
: feature_of >
{
};
// So that pot_tail_mean_prob can be automatically substituted
// with weighted_pot_tail_mean_prob when the weight parameter is non-void.
template
struct as_weighted_feature >
{
typedef tag::weighted_pot_tail_mean_prob type;
};
template
struct feature_of >
: feature_of >
{
};
}} // namespace boost::accumulators
#endif