c# - 求数学网络矩阵的协方差

Question

我正在使用 MathNet Numerics 矩阵并试图找出矩阵的协方差。

我们如何找到矩阵的协方差？

我们有方法可以在 Statistics 命名空间下找到两个 IEnumerable 之间的协方差。

http://numerics.mathdotnet.com/api/MathNet.Numerics.Statistics/Statistics.htm

但我不知道如何使用它来查找矩阵。

例如：在 matlab/octave 中

在 C# 中也是如此。我们如何实施？

Answer 1

我认为 Math.NET 中没有任何方法可以模仿covMatlab/Octave 中的函数，但您可以轻松编写自己的方法：

public static Matrix<double> GetCovarianceMatrix(Matrix<double> matrix)
{
    var columnAverages = matrix.ColumnSums() / matrix.RowCount;
    var centeredColumns = matrix.EnumerateColumns().Zip(columnAverages, (col, avg) => col - avg);
    var centered = DenseMatrix.OfColumnVectors(centeredColumns);
    var normalizationFactor = matrix.RowCount == 1 ? 1 : matrix.RowCount - 1;
    return centered.TransposeThisAndMultiply(centered) / normalizationFactor;
}

Answer 2

您可以使用 Accord.Statistics nuget 包中的 Tools.Covariance() 方法。http://accord-framework.net/docs/html/M_Accord_Statistics_Tools_Covariance.htm
做这样的事情：

using Accord.Statistics;
public void CalculateMatrixCovariance()
{
    var matrix = new double[,] {
        {3,5},
        {9,10}};

    var covMatrix = matrix.Covariance();

    Console.WriteLine(covMatrix[0,0] + " " + covMatrix[0, 1]);
    Console.WriteLine(covMatrix[1,0] + " " + covMatrix[1, 1]);
}