Hi all,
I've been struggling with the problem below for a few days now and am totally stumped. I am trying to use essGene to extract the data for essential member genes in my significantly differentially expressed gene sets and then geneData for the output and visualization of expression data
for these selected genes. I followed the vignette closely and have pretty much just replicated the code, but every time I run the code below I get the following error: Error in [.data.frame(gData, , icol) : undefined columns selected
Does anyone have any ideas what might be going on or how to fix this?
Many thanks,
Sarah
CODE:
cnts.kegg.p.CvsAC<- gage(cnts.norm, gsets = kegg.gs, ref = c(7:11), samp = c(13:17), compare ="unpaired", same.dir = TRUE)
rownames(cnts.kegg.p.CvsAC$greater)[1:3]
gs=uniqueunlistkegg.gs[rownames(cnts.kegg.p.CvsAC$greater)[1:3]]))
essData=essGene(gs, cnts.norm, ref = c(7:12), samp = c(13:17), compare ="unpaired")
head(essData, 10)
ref1=7:12
samp1=13:17
for (gs in rownames(cnts.kegg.p.CvsAC$greater)[1:3]) {
outname = gsub(" |:|/", "_", substr(gs, 10, 100))
geneData(genes = kegg.gs[[gs]], exprs = essData, ref = ref1, samp = samp1, outname = outname, txt = T, heatmap = T, Colv = F, Rowv = F, dendrogram = "none", limit = 3, scatterplot = T)
}