1

識別子の複数のリストを交差させ、リストのペア間のオーバーラップ数を含むテーブルを作成することに興味があります (「rk」と「t」)。私はsapplyが進むべき道であるという漠然とした考えを持っていますが、チュートリアルを検索して読んだ後でも行き詰まっています。

rk1 <- list("YH_sensitive_933","CS_sensitive_1294","YH_sensitive_944","JB_persistent_1224","CS_sensitive_1299","YY_sensitive_922", "YH_sensitive_952","YA_sensitive_949")   
rk2 <- list("YH_sensitive_944","JB_persistent_1224","CS_sensitive_1299","YY_sensitive_922", "YH_sensitive_952","YA_sensitive_949")
t1 <- list("YH_sensitive_933","CS_sensitive_1294","YH_sensitive_944")
t2 <- list("YH_sensitive_944","JB_persistent_1224")
t3 <- list("CS_sensitive_1299","YY_sensitive_922","YH_sensitive_944")
t4 <- list("YH_sensitive_952","YA_sensitive_949")

編集：リストの2つのリストをグループ化し、提案されているようにsapply/mapplyを試すのが最善だと思いました

F <- list(t1,t2,t3,t4)
G <- list(rk1,rk2)
> sapply(mapply(intersect,F,G), length)
[1] 3 2 3 2

しかし、私はRの初心者であり、ループと適用機能の使用に関するガイダンスを本当に感謝しています。しかし、私は rk1 の交点しか見えません (ただし、1 2 3 2 である必要がある rk2 は見えません)。

4

2 に答える 2

1

使用するlapply/sapply

 F <- list(t1, t2, t3, t4)
 G <- list(rk1, rk2)  

res <- do.call(`c`,setNames(lapply(G, function(.y) 
         setNames(sapply(F, `intersect`, .y), paste0("t",1:4))), paste0("rk",1:2)))

 sapply(res, length)
 #rk1.t1 rk1.t2 rk1.t3 rk1.t4 rk2.t1 rk2.t2 rk2.t3 rk2.t4 
 #    3      2      3      2      1      2      3      2 

res$rk1.t1
#[[1]]
#[1] "YH_sensitive_933"

#[[2]]
#[1] "CS_sensitive_1294"

#[[3]]
#[1] "YH_sensitive_944"


intersect(rk1,t1)
#[[1]]
#[1] "YH_sensitive_933"

#[[2]]
#[1] "CS_sensitive_1294"

#[[3]]
#[1] "YH_sensitive_944"

res$rk2.t1
# [[1]]
#[1] "YH_sensitive_944"

 intersect(rk2, t1)
 #[[1]]
 #[1] "YH_sensitive_944"

または、使用できますmapply（@Richard Scrivenのコメントからの基本的なアイデア）

 dat1 <- expand.grid(ls(pattern="^rk"), ls(pattern="^t"),stringsAsFactors=F)
 res1 <- mapply(intersect, mget(dat1[,1]), mget(dat1[,2]))

 res1[[1]]
 #[[1]]
 #[1] "YH_sensitive_933"

 #[[2]]
 #[1] "CS_sensitive_1294"

 #[[3]]
 #[1] "YH_sensitive_944"

アップデート

resをマトリックスに変換するには

  mat1 <- do.call(cbind,lapply(lapply(res, unlist),
                     `length<-`, max(sapply(res, length))))


  mat1
  #       rk1.t1              rk1.t2               rk1.t3             
  #[1,] "YH_sensitive_933"  "YH_sensitive_944"   "CS_sensitive_1299"
  #[2,] "CS_sensitive_1294" "JB_persistent_1224" "YY_sensitive_922" 
  #[3,] "YH_sensitive_944"  NA                   "YH_sensitive_944" 
  #     rk1.t4             rk2.t1             rk2.t2              
  #[1,] "YH_sensitive_952" "YH_sensitive_944" "YH_sensitive_944"  
  #[2,] "YA_sensitive_949" NA                 "JB_persistent_1224"
  #[3,] NA                 NA                 NA                  
  #      rk2.t3              rk2.t4            
  #[1,] "CS_sensitive_1299" "YH_sensitive_952"
  #[2,] "YY_sensitive_922"  "YA_sensitive_949"
  #[3,] "YH_sensitive_944"  NA

Update2

length出力を行列で取得する必要がある場合は、

  resL <- sapply(res,length)
  m1 <-  matrix(resL, nrow=2, byrow=TRUE,
          dimnames=list(paste0("rk", 1:2), paste0("t",1:4)))

  m1
  #   t1 t2 t3 t4
  #rk1  3  2  3  2
  #rk2  1  2  3  2
于 2014-10-17T18:04:58.110 に答える