Search
Question: DEXSeqHTML and biomart object
0
11 months ago by
ZheFrench10
ZheFrench10 wrote:

Could you provide a more detailed example of code to use biomart object and extra columns to add for example GeneSymbol in the output of DEXSeqHTML  ?

Not clear to me how you use mart  and attributes parameters in DEXSeqHTML

modified 11 months ago by Alejandro Reyes1.6k • written 11 months ago by ZheFrench10
1
11 months ago by
Alejandro Reyes1.6k
Dana-Farber Cancer Institute, Boston, USA
Alejandro Reyes1.6k wrote:

Hi! The easiest way to do that would be something along these lines:

library(DEXSeq)
library(biomaRt)

dxr <- DEXSeq( dxd )

mart <- useMart( biomart="ensembl", dataset="dmelanogaster_gene_ensembl" )
extraCols <- getBM(attributes=c("ensembl_gene_id", "ensembl_peptide_id"), filter="ensembl_gene_id",
mart=mart, values=geneIDs(dxd) )
extraCols <- extraCols[!duplicated(extraCols$ensembl_gene_id),] rownames(extraCols) <- extraCols$ensembl_gene_id
DEXSeqHTML( dxr, extraCols=extraCols )

1

Ok so this looks like the answer to my question, sorry I didn't see it first.  I am confused, how do you know for sure that the 'ensembl_gene_id' lines up with the 'ensembl_gene_id' in dxr?  I guess because you pull from it  initially with values=geneIDs(dxd)?  Thanks!!

Names are matched based on the rownames of extraCols. Thus, the rownames of extraCols should match what is stored in geneIDs().

Another minimal example:

library(DEXSeq)
data(pasillaDEXSeqDataSet, package="pasilla")
dxr <- DEXSeq( dxd )
extraCols <- data.frame( ids=unique(geneIDs(dxd)), mygeneIDs=sprintf("gene%0.3d", seq_along(unique(geneIDs(dxd)))))
rownames(extraCols) <- extraCols\$ids
DEXSeqHTML( object=dxr, extraCols=extraCols )