r - 表中行的成对组合？

Question

假设如下表：

X = 

        col1    col2    col3
row1    "A"      "0"     "1"
row2    "B"      "2"     "NA"
row3    "C"      "1"     "2"

我使用下面的代码选择两行的组合：

pair <- apply(X, 2, combn, m=2)

这将返回以下形式的矩阵：

pair = 

 [,1] [,2] [,3]
[1,] "A"  "0"  "1" 
[2,] "B"  "2"  NA  
[3,] "A"  "0"  "1" 
[4,] "C"  "1"  "2" 
[5,] "B"  "2"  NA  
[6,] "C"  "1"  "2" 

我希望迭代对，一次获取两行，即首先隔离 [1,] 和 [2 ,]，然后是 [3,] 和 [4,]，最后是 [5,] 和 [6,]。然后，这些行将作为参数传递给回归模型，即 lm(Y ~ row[i]*row[j])。

我正在处理一个大数据集。有人能建议如何一次迭代矩阵两行，将这些行分配给变量并作为参数传递给函数吗？

谢谢， S;-)

Answer 1

没有必要像这样将矩阵的行相乘，如果你有一个很大的数据集，可能会出现问题。相反，只需为每个实例挑选相关行即可。但预先创建选择很方便，可能是这样的：

xselect <- combn(1:nrow(X),2)

用您的数据进行说明（假设您只使用第 2 列和第 3 列）：

X <- matrix(c("A", "B", "C", 0,2,1,1,NA,2),3,3)
Y <- rnorm(2, 4, 2)

for (i in 1:ncol(xselect))
{
  x1 <- as.numeric(X[xselect[1,i], c(2,3)])
  x2 <- as.numeric(X[xselect[2,i], c(2,3)])
  print(lm(Y ~ x1 * x2))
}

Answer 2

我不确定您到底想对线性模型做什么，但要迭代 X（一次一对行），为每对生成一个因子，然后使用 < code>by

fac <- as.factor(sort(rep(1:(nrow(X)/2), 2)))
by(X, fac, FUN)

其中 FUN 是您想要对 X 中的行对应用的任何函数。