如何使用c++获得线性回归线的斜率？

我需要获得线性回归的斜率，类似于以下链接中 Excel 函数的实现方式：

http://office.microsoft.com/en-gb/excel-help/slope-function-HP010342903.aspx

是否有 C++ 库或某人创建的简单编码解决方案可以做到这一点？

我已经根据这个公式实现了代码，但是它并不总是给我正确的结果（取自这里http://easycalculation.com/statistics/learn-regression.php) ....

Slope(b) = (NΣXY - (ΣX)(ΣY)) / (NΣX2 - (ΣX)2)
         = ((5)*(1159.7)-(311)*(18.6))/((5)*(19359)-(311)2)
         = (5798.5 - 5784.6)/(96795 - 96721)
         = 13.9/74
         = 0.19 

如果我尝试针对以下向量，我会得到错误的结果（我应该期望 0.305556）： x = 6,5,11,7,5,4,4 y = 2,3,9,1,8,7,5

提前致谢。

这是一个 C++11 实现：

#include <algorithm>
#include <iostream>
#include <numeric>
#include <vector>

double slope(const std::vector<double>& x, const std::vector<double>& y) {
    const auto n    = x.size();
    const auto s_x  = std::accumulate(x.begin(), x.end(), 0.0);
    const auto s_y  = std::accumulate(y.begin(), y.end(), 0.0);
    const auto s_xx = std::inner_product(x.begin(), x.end(), x.begin(), 0.0);
    const auto s_xy = std::inner_product(x.begin(), x.end(), y.begin(), 0.0);
    const auto a    = (n * s_xy - s_x * s_y) / (n * s_xx - s_x * s_x);
    return a;
}

int main() {
    std::vector<double> x{6, 5, 11, 7, 5, 4, 4};
    std::vector<double> y{2, 3, 9, 1, 8, 7, 5};
    std::cout << slope(x, y) << '\n';  // outputs 0.305556
}

您可以添加数学要求的测试（x.size() == y.size() and x不是恒定的）或者，如上面的代码，假设用户会处理这个问题。