2017-09-29 63 views
3

我想完成以下任務,而不必輸入for循環,而是在單個apply()命令中輸入。將一個列表粘貼到一個矢量,重複每個矢量級別的列表

我有一個列表a,我想重複N次,其中N是矢量b的長度,a每次重複粘貼到的b的元件。

到目前爲止,我已經做了以下MWE:

var <- paste("var", 1:4, sep="") 
treat <- c("A","B") 
spec <- paste("sp", 1:3, sep="") 
a <- combn(var, 2, simplify = FALSE)#this 6 times, for each treatment and species 
b <- do.call(paste, c(expand.grid(treat, spec), sep='.')) 
a1 <- lapply(a, paste, b[1], sep='.') 
a2 <- lapply(a, paste, b[2], sep='.') 
a3 <- lapply(a, paste, b[3], sep='.') 
a4 <- lapply(a, paste, b[4], sep='.') 
a5 <- lapply(a, paste, b[5], sep='.') 
a6 <- lapply(a, paste, b[6], sep='.') 
a.final <- c(a1,a2,a3,a4,a5,a6) 
a.final 

這將是最佳的,如果我能a之前粘貼b

請注意,我的出發點是3個向量:var,treatspec,所以請隨時更改此處的任何內容。

回答

4

選項1:我們可以在沒有任何apply()循環的情況下完成此操作。我們unlist()a列表,paste()它複製到b值,然後relist()它基於複製a列表。試試這個:

aa <- relist(
    paste(unlist(a), rep(b, each=sum(lengths(a))), sep="."), 
    rep.int(a, length(b)) 
) 

檢查:

identical(aa, a.final) 
# [1] TRUE 

選項1 ba現在,把b值出門前,只要將參數在paste()電話:

relist(
    paste(rep(b, each=sum(lengths(a))), unlist(a), sep = "."), 
    rep.int(a, length(b)) 
) 

選項2:此選項使用apply()循環。這裏我們用Map()來做一對一的粘貼。

ra <- rep(a, length(b)) 
aa2 <- Map(paste, ra, relist(rep(b, each=sum(lengths(a))), ra), sep = ".") 

檢查:

identical(aa2, a.final) 
# [1] TRUE 

選項2 ba只是交換傳遞給paste()無名Map()參數。

ra <- rep(a, length(b)) 
Map(paste, relist(rep(b, each=sum(lengths(a))), ra), ra, sep = ".") 
+0

就是這樣!非常感謝 – DaniCee

1

力求貼近OP的方法,這可以用嵌套lapply()使用匿名函數來解決:

unlist(lapply(b, function(x) lapply(a, function(y) paste(x, y, sep = "."))), 
     recursive = FALSE) 
[[1]] 
[1] "A.sp1.var1" "A.sp1.var2" 

[[2]] 
[1] "A.sp1.var1" "A.sp1.var3" 

[[3]] 
[1] "A.sp1.var1" "A.sp1.var4" 

... 

[[34]] 
[1] "B.sp3.var2" "B.sp3.var3" 

[[35]] 
[1] "B.sp3.var2" "B.sp3.var4" 

[[36]] 
[1] "B.sp3.var3" "B.sp3.var4" 

注意,b中的a前面粘貼。需要unlist()才能列出頂級列表。

要(在baa.final比較)驗證方式工作:

identical(a.final, 
      unlist(lapply(b, function(x) lapply(a, function(y) paste(y, x, sep = "."))), 
       recursive = FALSE)) 
[1] TRUE 
1

這是從頭開始創建標籤,並返回他們完全不同的方法在36行×2列data.table中,而不是具有長度爲2的36個向量的列表:

library(data.table) 
# cross join of treat, spec, var. Note, full labels will be created in sprintf() below 
DT <- CJ(LETTERS[1:2], 1:3, 1:4) 
# non equi join as replacement of combn() 
DT[DT, on = .(V1, V2, V3 > V3), nomatch = 0L, 
    # create labels 
    .(sprintf("%s.sp%s.var%i", V1, V2, V3), 
    sprintf("%s.sp%s.var%i", V1, V2, x.V3))] 
  V1   V2 
1: A.sp1.var1 A.sp1.var2 
2: A.sp1.var1 A.sp1.var3 
3: A.sp1.var1 A.sp1.var4 
4: A.sp1.var2 A.sp1.var3 
5: A.sp1.var2 A.sp1.var4 
6: A.sp1.var3 A.sp1.var4 
7: A.sp2.var1 A.sp2.var2 
... 
29: B.sp2.var2 B.sp2.var4 
30: B.sp2.var3 B.sp2.var4 
31: B.sp3.var1 B.sp3.var2 
32: B.sp3.var1 B.sp3.var3 
33: B.sp3.var1 B.sp3.var4 
34: B.sp3.var2 B.sp3.var3 
35: B.sp3.var2 B.sp3.var4 
36: B.sp3.var3 B.sp3.var4 
      V1   V2 
相關問題