有没有办法让一个spearman的相关矩阵超过两个数组?

时间:2013-10-02 14:25:58

标签: mathdotnet

我有生成N个双打数组的代码,我想创建一个spearmans相关矩阵。是否有功能为我做,或者我必须遍历所有组合并使用Correlation.Spearman()构建我自己的相关矩阵?

1 个答案:

答案 0 :(得分:1)

目前还没有计算相关矩阵的例行程序,但是我已经打开了票证#161来跟踪GitHub。

与此同时,您可以使用以下例程(我在这里使用v3.0.0-alpha5):

Matrix<double> Spearman(double[][] data)
{
    var m = Matrix<double>.Build.DenseIdentity(data.Length);
    for(int i=0; i<data.Length; i++)
    for(int j=i+1; j<data.Length; j++)
    {
        var c = Correlation.Spearman(data[i], data[j]);
        m.At(i,j,c);
        m.At(j,i,c);
    }
    return m;
}

var vectors = new[] {
    new[] { 1.0, 2.0, 3.0, 4.0 },
    new[] { 2.0, 4.0, 6.0, 8.0 },
    new[] { 4.0, 3.0, 2.0, 1.0 },
    new[] { 0.0, 10.0, 10.0, 20.0 },
    new[] { 2.0, 4.0, -4.0, 2.0 }
};

Spearman(vectors);

将返回:

DenseMatrix 5x5-Double
           1            1           -1     0.948683    -0.316228
           1            1           -1     0.948683    -0.316228
          -1           -1            1    -0.948683     0.316228
    0.948683     0.948683    -0.948683            1            0
   -0.316228    -0.316228     0.316228            0            1

更新2013-10-20:

添加到master,可在V3.0.0-alpha6和更新版本中使用:

  • Correlation.PearsonMatrix(params double[][] vectors)
  • Correlation.PearsonMatrix(IEnumerable<double[]> vectors)
  • Correlation.SpearmanMatrix(params double[][] vectors)
  • Correlation.SpearmanMatrix(IEnumerable<double[]> vectors)