R中配置文件之间的马氏距离

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 分的相似性。我觉得马氏距离是我需要的度量，我检查了一些包（HDMD、StatMatch）但找不到合适的函数。

Answer 1

一种方法是创建差异分数矩阵，然后根据差异分数计算马氏距离。

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

R中配置文件之间的马氏距离

Mahalanobis distance between profiles in R

r

distance

profiles

mahalanobis