在R中創建相關矩陣

Question

我有一個文件，其中包含500行（二進制分數）和120列的矩陣。 該文件是0和1的簡單矩陣。

>file

00010010101010
01001010100101
00101001010001
11110101001010

我正在編寫一個函數，該函數使用特殊的相關公式來查找行之間的此相關。 它需要兩個向量行作為輸入fn（row1，row2）。 例如。 row1和row2並計算此特殊相關性。

例

>fn(file[1,], file[2,])
>0.32

我可以對兩行執行此操作，但是如何為所有行創建500x500的相關矩陣。 有人可以幫忙嗎？ 謝謝。

Answer 1

嘗試這個：

corr.mat <- outer(seq_len(nrow(file)), seq_len(nrow(file)),
                  Vectorize(function(i, j) fn(file[i,], file[j,])))

如果這對於您的需求來說太慢了，也許有更好的方法，但是您必須證明fn應該做什么。

PS： file是R中函數的名稱，應避免將其用於自己的變量。