Question

我正在尝试创建一个可以变化的列表，并且可以根据不同的用途进行互换。我希望它的格局为：

P1A1 P1A2 P1A3 P1A4 P1A5 P1A6 
P2A1 P2A2 P2A3 P2A4 P2A5 P2A6
P3A1 P3A2 P3A3 P3A4 P3A5 P3A6

P从1变为num.Patches，A从1变为num.allele。因此，对于上面的示例，num.Patches = 3和num.alleles = 6

我正在尝试使用for循环：

pdum <- matrix(data=NA,nrow=1,ncol=num.Patches)
Adum <- matrix(data=NA,nrow=1,ncol=num.allele*num.Patches+1)  
key2 <- matrix(data=NA,nrow=1,ncol=num.allele*num.Patches+1)

for (i in 1:num.Patches) pdum[1,i] <- matrix(paste("P",i,sep=""))
pdum2 <- as.character(pdum)

for (k in 1:num.Patches){
  for (i in pdum2) {
    for (j in 1:num.allele){
     Adum[1,k+num.allele*(1-j)] <- matrix(paste(i,"A",j,sep=""))
    }
  }
}

我遇到麻烦的是：

Adum[1,k+num.allele*(1-j)] <- matrix(paste(i,"A",j,sep=""))

我不知道如何引用Adum矩阵的每个条目并用特定值填充它。最终目标是创建一个列表，用作较大矩阵的列名，以便可以轻松引用它。

非常感谢你。

Answer 1

要获得1：num.Patches和1：num.alleles的所有组合然后合并，您可以尝试：

combinations <- expand.grid(1:num.Patches, 1:num.alleles)
names <- paste0("P", combinations[,1], "A", combinations[,2])

这将为您提供所需的模式（作为向量）。您可以稍后将它们设置为您创建的某个矩阵的列名。

Answer 2

您正在寻找outer：

num.Patches <- 3
num.alleles <- 6

mm <- outer(seq(num.Patches), seq(num.alleles),
            function(x,y)paste0('P',x,'A',y))

     [,1]   [,2]   [,3]   [,4]   [,5]   [,6]  
[1,] "P1A1" "P1A2" "P1A3" "P1A4" "P1A5" "P1A6"
[2,] "P2A1" "P2A2" "P2A3" "P2A4" "P2A5" "P2A6"
[3,] "P3A1" "P3A2" "P3A3" "P3A4" "P3A5" "P3A6"

如果你想得到一个名字的载体：

as.vector(mm)

在R中为列名创建列表

2 个答案: