c++ - 使用 boost 查找数据集的峰度

Question

我有一个数据 vector ，我想找出数据集的峰度。我想用 Boost 这样做，这是我目前所拥有的(不可编译):

#include <boost/math/distributions.hpp>
using namespace std;

int main()
{
    vector<double> a;
    a.push_back(-1);
    a.push_back(0);
    a.push_back(1);

    cout << "Kurtosis:"<< kurtosis(a) << endl;
    return 0;
}

为什么这行不通？我的编译器给我错误:“[...]\main.cpp|28|error: 'kurtosis' was not declared in this scope|”

Answer 1

对于一个你没有包括 kurtosis 的标题:

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

即使您这样做了，正如您所见，它不适用于直接的 vector，您可能想要做的是使用 accumulator_set 和更多 header 。

这是一个使用accumulator_set 的最小示例，它展示了解决问题的两种方法:

#include <boost/math/distributions.hpp>
#include <boost/accumulators/accumulators.hpp>
#include <boost/accumulators/statistics/stats.hpp>
#include <boost/accumulators/statistics/mean.hpp>
#include <boost/accumulators/statistics/variance.hpp>
#include <boost/accumulators/statistics/kurtosis.hpp>
#include <iostream>
#include <vector>

using namespace boost::accumulators;

int main()
{
    accumulator_set<double, stats<tag::mean, tag::kurtosis > > acc;
    accumulator_set<double, stats<tag::mean, tag::kurtosis > > acc2;

    acc(2) ;
    acc(3) ;
    acc(4) ;

    std::cout << mean(acc) << " " << kurtosis(acc) << std::endl ;

    std::vector<double> v1 ;

    v1.push_back(2);
    v1.push_back(3);
    v1.push_back(4);

    acc2 = std::for_each( v1.begin(), v1.end(), acc2 ) ;

    std::cout << mean(acc2) << " " << kurtosis(acc2) << std::endl ;
}

这是指向累加器框架的链接 User's Guide .本指南有一些很好的例子。

上一个thread找到了一种使用 vector 的方法，尽管它根本不是直截了当的，而且我无法让它工作。