如何计算2个给定向量之间的Pearson相关性?

时间:2016-02-29 14:35:25

标签: c# correlation pearson pearson-correlation

我必须在C#中编写代码

您能否在下面给出的示例中逐步解释?

vector 1 : [0.3, 0, 1.7, 2.2]
vector 2 : [0, 3.3, 1.2, 0]
非常非常

这将用于文档聚类

1 个答案:

答案 0 :(得分:6)

我在 Java

上的答案的改编

How to find correlation between two integer arrays in java

代表C#。首先,Pearson Correlation是

http://en.wikipedia.org/wiki/Correlation_and_dependence

假设两个向量(让它们为IEnumerable<Double>)具有相同的长度

  private static double Correlation(IEnumerable<Double> xs, IEnumerable<Double> ys) {
    // sums of x, y, x squared etc.
    double sx = 0.0;
    double sy = 0.0;
    double sxx = 0.0;
    double syy = 0.0;
    double sxy = 0.0;

    int n = 0;

    using (var enX = xs.GetEnumerator()) {
      using (var enY = ys.GetEnumerator()) {
        while (enX.MoveNext() && enY.MoveNext()) {
          double x = enX.Current;
          double y = enY.Current;

          n += 1;
          sx += x;
          sy += y;
          sxx += x * x;
          syy += y * y;
          sxy += x * y;
        }
      }
    }

    // covariation
    double cov = sxy / n - sx * sy / n / n;
    // standard error of x
    double sigmaX = Math.Sqrt(sxx / n -  sx * sx / n / n);
    // standard error of y
    double sigmaY = Math.Sqrt(syy / n -  sy * sy / n / n);

    // correlation is just a normalized covariation
    return cov / sigmaX / sigmaY;
  }

测试:

  // -0.539354840012899
  Double result = Correlation(
    new Double[] { 0.3, 0, 1.7, 2.2 }, 
    new Double[] { 0, 3.3, 1.2, 0 });