Description Usage Arguments Author(s) Examples
A function to associate EG and name symbols to rsem-generated-gene-tables, using UCSC as annotation
1 |
my.path |
Name of the path of interest |
path.db |
Name of the pathway to be used for plotting |
type |
Entrez geneid or Symbol |
Raffaele A Calogero
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 | # samples.dir <- paste(find.package(package="thu12feb25"),"/examples/", sep="")
# samples <- dir(samples.dir)[grep(".genes.results",dir(samples.dir))]
# counts <- lapply(paste(samples.dir,samples, sep=""), function(x){
# tmp <- annotatingGenes(filename=x, org="mm9")
# counts.tmp <- tmp$expected_count
# })
# counts.df <- t(data.frame(matrix(unlist(counts), nrow=length(samples), byrow=T)))
# tmp <- annotatingGenes(filename=paste(samples.dir,samples[1], sep=""), org="mm9")
# dimnames(counts.df) <- list(as.character(tmp$Symbol),gsub(".genes.results","",samples))
# counts.df1 <- counts.df[rowSums(counts.df) > 10,]
# counts.df1 <- as.matrix(counts.df1)
#converting in uppercase gene symbols
# dimnames(counts.df1)[[1]] <- toupper(dimnames(counts.df1)[[1]])
# deg <- DEGraph(counts.df1, conditions, panther, type="RNASeq", IDs = "symbol", norm.method="DESeq2")
# names(which(deg$res[,1] <= 0.05))
#3
# names(panther)[which(names(panther)%in%names(which(deg$res[,1] <= 0.05)))]
# plotPath(my.path="mRNA splicing", path.db="panther", type="symbol")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.