是否可以将 boost 累加器与向量一起使用？

Question

我想使用 boost 累加器来计算向量变量的统计数据。有没有一种简单的方法可以做到这一点。我认为不可能使用最愚蠢的东西：

  using namespace boost::accumulators;
  //stuff...

  accumulator_set<vector<double>, stats<tag::mean> > acc;
  vector<double> some_vetor;
  //stuff 
  some_vector = doStuff();
  acc(some_vector);

也许这是显而易见的，但我还是尝试了。 :P

我想要的是有一个累加器来计算一个向量，该向量是许多向量分量的平均值。有简单的出路吗？

编辑：

我不知道我是否完全清楚。我不想要这个：

 for_each(vec.begin(), vec.end(),acc); 

这将计算给定向量的条目的平均值。我需要的是不同的。我有一个会吐出向量的函数：

 vector<double> doSomething(); 
 // this is a monte carlo simulation;

我需要运行多次并计算这些向量的向量平均值：

  for(int i = 0; i < numberOfMCSteps; i++){
  vec = doSomething();
  acc(vec);
  }
  cout << mean(acc);

我希望mean(acc)本身就是一个向量，其条目[i ] 将是累积向量的条目 [i] 的平均值。

Boost 的文档中有对此的暗示，但没有明确的。而且我有点傻。 ：P

Answer 1

我稍微研究了一下你的问题，在我看来，Boost.Accumulators 已经提供了对 std::vector 的支持。这是我可以在 用户指南的一部分：

另一个例子，其中数字
运算符子库很有用的是
当类型没有定义时
使用它所需的运算符重载
用于一些统计计算。
例如，std::vector 不会重载任何算术运算符，但
使用 std::vector 可能很有用
作为样本或变量类型。这
数字运算符子库定义
必要的运算符重载
boost::numeric::operators
命名空间，被纳入范围
通过累加器框架
使用指令。

验证，文件 boost/accumulators/numeric/functions/vector.hpp 确实包含“简单”解决方案所需的运算符去工作。

我相信你应该尝试：

• 包括
  • boost/accumulators/numeric/function/vector.hpp 在任何其他累加器标头之前
  • 定义BOOST_NUMERIC_FUNCTIONAL_STD_VECTOR_SUPPORT时
  • boost/accumulators/numeric/functional.hpp
• 将运算符纳入 using namespace boost::numeric::operators; 的范围。

只剩下最后一个细节：执行将在运行时中断，因为初始累积值是默认构造的，并且当尝试将大小为 n 的向量添加到空向量时会发生断言。为此，似乎您应该使用以下方式初始化累加器（其中 n 是向量中的元素数量）：

accumulator_set<std::vector<double>, stats<tag::mean> > acc(std::vector<double>(n));

我尝试了以下代码，mean 给了我一个 std::vector 大小为 2：

int main()
{
    accumulator_set<std::vector<double>, stats<tag::mean> > acc(std::vector<double>(2));

    const std::vector<double> v1 = boost::assign::list_of(1.)(2.);
    const std::vector<double> v2 = boost::assign::list_of(2.)(3.);
    const std::vector<double> v3 = boost::assign::list_of(3.)(4.);
    acc(v1);
    acc(v2);
    acc(v3);

    const std::vector<double> &meanVector = mean(acc);
}

我相信这就是您想要的？

Answer 2

我现在还没有设置它来尝试，但如果所有 boost::accumulators 需要的是正确定义的数学运算符，那么您可能可以使用不同的向量类型： http://www.boost.org/doc/libs/1_37_0/libs/numeric/ ublas/doc/vector.htm

Answer 3

那么文档怎么样？

// The data for which we wish to calculate statistical properties:
std::vector< double > data( /* stuff */ );

// The accumulator set which will calculate the properties for us:    
accumulator_set< double, features< tag::min, tag::mean > > acc;

// Use std::for_each to accumulate the statistical properties:
acc = std::for_each( data.begin(), data.end(), acc );