多柱双向交叉稳定

Question

我有病人诊断的数据。有13列，每列具有不同的疾病类型，它们都存储二进制值(0-无疾病，1-疾病存在)。我想看看这些疾病的共患情况，也就是看看在同样患有type2的病人中存在多少次1型疾病。我的数据(简化)：

ID <- sample(10:50, 20)
type1 <- sample(0:1, 20, replace = T)
type2 <- sample(0:1, 20, replace = T)
type3 <- sample(0:1, 20, replace = T)
type4 <- sample(0:1, 20, replace = T)
type5 <- sample(0:1, 20, replace = T)
type6 <- sample(0:1, 20, replace = T)

data <- cbind.data.frame(ID, type1, type2, type3, type4, type5, type6)

我正在寻找类似于cor输出的东西(这将给我计数而不是频率/相关测量)。

Answer 1

试试这个：

data_matrix <- as.matrix(data[, 2:ncol(data)])
cooccurrence <- crossprod(data_matrix)
diag(cooccurrence) <- 0

这是一个对称矩阵。只取上三角形或下三角形。