是否可以要求 R 从一组二元变量中随机选择两个变量来计算比例？

Question

我在数据集中有几个二进制变量，我想计算两个变量的所有组合中“1 1”的比例。例如（a1=1 和 a2=1）的比例。我可以每次手动指定两个变量来运行代码，但是我的数据中有超过 10 个变量，因此至少会有 45 个组合。有没有办法让 R 自动配对变量为我计算？

structure(list(a1 = structure(c(2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L), 
                              .Label = c("0", "1"), class = "factor"), 
               a2 = structure(c(1L,1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L), 
                              .Label = c("0", "1"), class = "factor"),
               a3 = structure(c(1L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 1L,1L, 1L), 
                              .Label = c("0", "1"), class = "factor"), 
               a4 = structure(c(1L,2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), 
                              .Label = c("0","1"), class = "factor"), 
               a5 = structure(c(2L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L), 
                              .Label = c("0", "1"), class = "factor")), row.names = 180:190, class = "data.frame")

Answer 1

我们可以创建一个函数来获取sample的列名选择两个，然后对数据进行子集化，检查两个列是否等于1，得到mean

f1 <- function(dat) {
    nm1 <- sample(names(dat), 2, replace = FALSE)
    setNames(mean(dat[,nm1[1]]== 1 & dat[,nm1[2]] == 1), paste(nm1, collapse="_"))
   }

f1(df1)
# a3_a5 
#   0 

---

如果我们想要所有的组合

f1 <- function(dat) {
       combn(names(dat), 2, FUN = function(nm) {
              nm1 <- paste(nm, collapse="_")
              setNames(mean(dat[, nm[1]] ==1 & dat[, nm[2]] == 1), nm1)},
  simplify = FALSE) 
   }

f1(df1)
#[[1]]
#a1_a2 
#    0 

#[[2]]
#a1_a3 
#    0 

#[[3]]
#a1_a4 
#    0 

#[[4]]
#     a1_a5 
#0.09090909 

#[[5]]
#     a2_a3 
#0.09090909 

#[[6]]
#a2_a4 
#    0 

#[[7]]
#     a2_a5 
#0.09090909 

#[[8]]
#a3_a4 
#    0 

#[[9]]
#a3_a5 
#    0 

#[[10]]
#a4_a5 
#    0