R中轮廓之间的Mahalanobis距离

Question

100名受试者的样本回应了两项性格测试。 这些测试的措辞略有不同，但通常是相同的，即它们都衡量相同的4种态度。 因此，我有2个这样的矩阵，每个科目有4个分数：

>test1
subj     A1     A2     A3     A4
   1  -2.14   1.21   0.93  -1.72
   2   0.25   1.17   0.67   0.67

>test2
subj     A1     A2     A3     A4
   1  -1.99   1.11   1.00  -1.52
   2   0.24   1.20   0.71   0.65

我想评估两个测试中的概况的相似性，即每个人的两组4个分数的相似性。 我觉得mahalanobis距离是我需要的量度，我检查了一些包（HDMD，StatMatch），但找不到合适的功能。

Answer 1

一种方法是创建差异分数矩阵，然后计算差异分数上的马哈拉诺比斯距离。

testDiff <- test1 - test2
testDiffMahalanobis <- mahalanobis(testDiff, 
                                   center = colMeans(testDiff),
                                   cov = cov(testDiff))