【问题标题】:Unique combination of columns独特的列组合
【发布时间】:2011-07-02 03:58:14
【问题描述】:

这是我的简化数据集:

foo <- data.frame(var1= c(1:10), var2=rep(1:5,2),var3=rep(1:2,5),var4=rep(3:7,2) )    

总共 20 个变量

foo

   var1 var2 var3   var4    ... var20
1     1    1    1      3
2     2    2    2      4
3     3    3    1      5
4     4    4    2      6
5     5    5    1      7
6     6    1    2      3
7     7    2    1      4
8     8    3    2      5
9     9    4    1      6
10   10    5    2      7

我需要获得 3 个变量的唯一组合及其每个时期的总和

即。某样东西

  var1var2var3   var1var3var4   var1var5var18  etc...
1     6               sum
2     6           
3     7          
4     10           
5     11           
6     9         
7     10             
8     13          
9     14          
10    17         

注意 var1var3var5 与 var3var1var5 相同

【问题讨论】:

  • 您是在求解任意数据集还是这个特定数据集?这个可能有优化。此外,您的第一个解决方案数据点似乎应该是 3,而不是 6(第 1 行,第一列)
  • 看来combn(20,3) 在这里可以用来生成所有列的唯一组合。如何将其转化为引用 applyfor loop 中每一列的内容有点超出我的理解。

标签: r matrix combinations


【解决方案1】:

正如@Chase 所建议的那样,combn 可以为您提供您想要的:

nams <- apply( combn(colnames(foo),3), 2, function(z) paste(z, collapse = ''))
cols <- combn( ncol(foo), 3)

tripleSums <- apply( cols, 2, function(z) rowSums(foo[,z]))
colnames(tripleSums) <- nams

> tripleSums
      var1var2var3 var1var2var4 var1var3var4 var2var3var4
 [1,]            3            5            5            5
 [2,]            6            8            8            8
 [3,]            7           11            9            9
 [4,]           10           14           12           12
 [5,]           11           17           13           13
 [6,]            9           10           11            6
 [7,]           10           13           12            7
 [8,]           13           16           15           10
 [9,]           14           19           16           11
[10,]           17           22           19           14

【讨论】:

  • 太好了,这正是我想要的。谢谢
  • @Alex,尝试接受这个答案作为正确答案(点击投票分数下的复选标记)。
猜你喜欢
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
  • 2020-09-03
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
相关资源
最近更新 更多