mIntersect {PhosR} | R Documentation |
A recusive loop for intersecting multiple sets.
mIntersect(x, y, ...) mUnion(x, y, ...)
x, y, ... |
objects to find intersection/union. |
An intersection/union of input parameters
data('phospho_liverInsTC_RUV_sample') data('phospho_L6_ratio') site1 <- gsub('~[STY]', '~', sapply(strsplit(rownames(phospho.L6.ratio), '~'), function(x){paste(toupper(x[2]), x[3], sep='~')})) site2 <- rownames(phospho.liver.Ins.TC.ratio.RUV) # step 2: rank by fold changes tmp <- do.call(cbind, lapply(split(1:ncol(phospho.L6.ratio), gsub('_exp\\d+', '', colnames(phospho.L6.ratio))), function(i){rowMeans(phospho.L6.ratio[,i])})) site1 <- t(sapply(split(data.frame(tmp), site1), colMeans))[,-1] tmp <- do.call(cbind, lapply(split(1:ncol(phospho.liver.Ins.TC.ratio.RUV), gsub( '(Intensity\\.)(.*)(\\_Bio\\d+)', '\\2', colnames(phospho.liver.Ins.TC.ratio.RUV))), function(i){ rowMeans(phospho.liver.Ins.TC.ratio.RUV[,i]) })) site2 <- t(sapply(split(data.frame(tmp), site2), colMeans)) o <- mIntersect(site1, site2)