2

假设我有一个集合 N={1,2,3}，那么我们可以列出它所有的 7 个非空子集。

    n=3 # number of elements in a set
    a=2^n-1 # number of non-empty subsets for that set
    subsets=lapply(1:n, function(x) combn(n, x)) # list all the non-empty subest
    subsets

现在我想将这些子集放入一个矩阵并组织如下：

    if n=3     or in an index matrix:    
    1 0 0      1 0 0
    0 2 0      0 1 0
    0 0 3      0 0 1

    1 2 0      1 1 0
    1 0 3      1 0 1
    0 2 3      0 1 1

    1 2 3      1 1 1

任何人都知道如何编写可以轻松扩展到任何 n (=4, 5, 6...) 的代码？我试过这个：

    subindex=matrix(c(0), nrow=a, ncol=n)

    i=1
      while(i<=a){

       j=n
       b=2^(n-1)
       N=i
          while(N>0){
             if(b<=N) {subindex[i,j]=1}&{N=N-b}
             b=trunc(b/2)
             j=j-1
          }

       i=i+1
       }
      subindex

但是我在第 3 行和第 4 行得到的索引矩阵是错误的。如果 n = 4，那么会有更多错误......有人可以纠正这个或简化这个代码吗？或者只是编写一个全新的代码。非常感谢。

4

1 回答 1

3 
n <- 4
lapply(seq_len(n), function(i)t(combn(n, i, FUN = tabulate, nbins = n)))

# [[1]]
#      [,1] [,2] [,3] [,4]
# [1,]    1    0    0    0
# [2,]    0    1    0    0
# [3,]    0    0    1    0
# [4,]    0    0    0    1
# 
# [[2]]
#      [,1] [,2] [,3] [,4]
# [1,]    1    1    0    0
# [2,]    1    0    1    0
# [3,]    1    0    0    1
# [4,]    0    1    1    0
# [5,]    0    1    0    1
# [6,]    0    0    1    1
# 
# [[3]]
#      [,1] [,2] [,3] [,4]
# [1,]    1    1    1    0
# [2,]    1    1    0    1
# [3,]    1    0    1    1
# [4,]    0    1    1    1
# 
# [[4]]
#      [,1] [,2] [,3] [,4]
# [1,]    1    1    1    1
于 2013-10-28T23:53:35.697 回答