Redundant gene sets & gagePipe
0
0
Entering edit mode
Luo Weijun ★ 1.6k
@luo-weijun-1783
Last seen 10 months ago
United States
Christian, You may check and remove the potential redundancy using esset.grp function on individual GAGE runs. Such redundancy processing has not been incorporated into the gagePipe procedure as it may not always be desirable for the users. Below is a function I wrote (warning: not fully tested!) to process the output of your gagePipe procedure. The arguments here follow the same definition of gagePipe. HTH Weijun gagePipe.esset.grp <- function(arraydata, dataname = "arraydata", sampnames, gsname = c("kegg.gs", "go.gs"), ref.list, samp.list, compare = "paired", ...) { load(paste(dataname, ".gage.RData", sep = "")) for (gs in gsname) { gsets = evalas.name(paste(gs))) for (i in 1:length(sampnames)) { ref = if (is.list(ref.list)) ref.list[[i]] else ref.list samp = samp.list[[i]] gage.p = evalas.name(paste(sampnames[i], paste(gs, ".p", sep = ""), sep = "."))) outname.up = paste(sampnames[i], paste(gs, ".up", sep = ""), sep = ".") assign(paste(sampnames[i], paste(gs, ".esg.up", sep = ""), sep = "."), esset.grp(gage.p$greater, arraydata, gsets = gsets, ref = ref, samp = samp, compare = compare, test4up = T, output = T, outname = outname.up, make.plot = F, ...)) outname.dn = paste(sampnames[i], paste(gs, ".dn", sep = ""), sep = ".") assign(paste(sampnames[i], paste(gs, ".esg.dn", sep = ""), sep = "."), esset.grp(gage.p$less, arraydata, gsets = gsets, ref = ref, samp = samp, compare = compare, test4up = T, output = T, outname = outname.dn, make.plot = F, ...)) } } save(list = objects(pattern = "esg.up$|esg.dn$"), file = paste(dataname, ".gage.esset.grp.RData", sep = "")) return(invisible(1)) } -------------------------------------------- On Sat, 11/30/13, Christian De Santis <christian.desantis at="" stir.ac.uk=""> wrote: Subject: Redundant gene sets & gagePipe To: "bioconductor at r-project.org" <bioconductor at="" r-project.org="">, "luo_weijun Date: Saturday, November 30, 2013, 8:08 AM Hi Weijun, ? Is there also a quick way to calculate the redundant gene-sets in batch when analysing multiple rounds of GAGE using gagePipe, or can this only be done individually on each of the contrasts? ? Thanks and have a nice weekend, Christian? ? The University of Stirling has been ranked in the top 12 of UK universities for graduate employment*. 94% of our 2012 graduates were in work and/or further study within six months of graduation. *The Telegraph The University of Stirling is a charity registered in Scotland, number SC 011159.
gage gage • 1.0k views
ADD COMMENT

Login before adding your answer.

Traffic: 775 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6