共分散行列

ここの説明を参考に
http://stattrek.com/matrix-algebra/covariance-matrix.aspx

列ベクトルで表現するか、行ベクトルで表現するかで結果が変わってくる。
MATLABのcov()関数のデフォルトは

試行回数1 カテゴリ１値　カテゴリ２値　カテゴリ３値
試行回数2 カテゴリ１値　カテゴリ２値　カテゴリ３値
試行回数3 カテゴリ１値　カテゴリ２値　カテゴリ３値
試行回数4 カテゴリ１値　カテゴリ２値　カテゴリ３値

ぽい。
教科書的な定義とは逆かも？

A=[
    90    60    90
    90    90    30
    60    60    60
    60    60    90
    30    30    30
];
N=size(A,1);
mu = ones(N)*A/N;
Va = ((A-mu)'*(A-mu))/N;

B = [
    90    90    60    60    30
    60    90    60    60    30
    90    30    60    90    30
    ];
N=size(B,2);
mu = B*ones(N)/N;
Vb = ((B-mu)*(B-mu)')/N;

Ca = cov(A,1);
Cb = cov(B',1);

Va
Vb
Ca
Cb

>>
Va =
   504   360   180
   360   360     0
   180     0   720
Vb =
   504   360   180
   360   360     0
   180     0   720
Ca =
   504   360   180
   360   360     0
   180     0   720
Cb =
   504   360   180
   360   360     0
   180     0   720