intersect multi-dimensional arrays in R -



intersect multi-dimensional arrays in R -

i interested in intersecting multiple lists of identifiers, making table number of overlaps between pairs of lists ('rk' vs 't'). have vague thought sapply way go still stuck after searching , reading tutorials.

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")

edit: thought maybe it'd best grouping 2 lists of lists , seek sapply/mapply suggested

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

but i'm r beginner , appreciate guidance on looping , using apply functionals. see intersections rk1 (but not rk2, should 1 2 3 2)

using 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"

or utilize mapply (basic thought @richard scriven's comment)

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" update

to convert res matrix

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

if need length output in matrix,

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

arrays r sapply mapply

Comments

Post a Comment

Popular posts from this blog

java Multi query from Mysql using netbeans -

java Multi query from Mysql using netbeans - i have java code retrieve info mysql appropriate code in netbeans.the code below. want multi query such illustration finding total number of products , on.could help me please. in advance import java.sql.*; public class javamysqlselectexample { public static void main(string[] args) { seek { string mydriver = "com.mysql.jdbc.driver"; string myurl = "jdbc:mysql://localhost/sample"; class.forname(mydriver); connection conn = drivermanager.getconnection(myurl, "root", "mypasscode1"); string query = "select * products"; statement st = conn.createstatement(); resultset rs = st.executequery(query); while (rs.next()==true) { int code = rs.getint("code"); string productname = rs.getstring("product_name"); string producttype = rs.getstring("product_type"); system.out.format(" %d, %s, %s\n", code, pr
Read more

c# - DotNetZip fails with "stream does not support seek operations" -

c# - DotNetZip fails with "stream does not support seek operations" - i using dotnetzip in c# unzip stream follows: public static void unzipfromstream(stream stream, string outdir) { //omit seek grab block using (zipfile zip = zipfile.read(stream)){ foreach (zipentry e in zip){ e.extract(outdir, extractexistingfileaction.overwritesilently); } } } stream obtained using webclient client = new webclient(); stream fs = client.openread(url); however, got next exception exception during extracting zip stream system.notsupportedexception: stream not back upwards seek operations. @ system.net.connectstream.get_position() @ ionic.zip.zipfile.read(stream zipstream, textwriter statusmessagewriter, encoding encoding, eventhandler`1 readprogress) on server side(asp.net mvc 4), returning filepathresult or filestreamresult both caused exception. should obtain stream differently on client side? or how create server homec
Read more

c++ - StartServiceCtrlDispatcher don't can access 1063 error -

c++ - StartServiceCtrlDispatcher don't can access 1063 error - i write on c++ service , have error, , don't can fixed it. my service main function int main (int argc, tchar *argv[]) {dword f; for(int i=0;i<113;i++) f = getlasterror(); servicepath = lptstr(argv[0]); outputdebugstring(_t("my sample service: main: entry")); //installservice(); service_table_entry servicetable[] = { {service_name, (lpservice_main_function) servicemain}, {null, null} }; // startservice(); if(!startservicectrldispatcher(servicetable)) { f = getlasterror(); //addlogmessage("error: startservicectrldispatcher"); } else if( memcmp(argv[argc-1],"install",7)) { installservice(); } else if( memcmp(argv[argc-1],"remove",6)) { removeservice(); } else if( memcmp(argv[argc-1],"start",5)) { startservice(); } else if( memcmp(argv[argc-1],"stop",4))
Read more