[BioC] DexSeq extracols option
Alejandro Reyes
alejandro.reyes at embl.de
Mon Jul 29 09:42:16 CEST 2013
Dear Ioannis,
Thanks for your message, it is usually nice and easier to track the
problem if you add the output of your sessionInfo and reproducible code,
but I am guessing that you are using a old version of DEXSeq.
The table that the HTML reporter produces a table with information for
each gene, not for each exon. I would be almost sure that the 396 are
the number of genes that have significant exons. As the documentation
indicates, the parameter extraCols accepts a data frame in which the
rownames are genes of your ExonCountSet object.
Let me know if this is the problem,
Best regards,
Alejandro
> Hello,
>
> A (hopefully) simple question.
>
> As I saw in ??DEXSeqHTML, there is an option to add attributes in the HTML
> results (extracols).
>
> I downloaded from biomaRt a set aof attributes to test this (gene name,
> biotype and description), and created a data frame having as a first column
> the exons of my ecs (retaining exon order), geneID and the new attributes.
>
> The first two columns were taken by accessing fData from the ecs.
>
> The resulting data frame naturally had an equal number of rows with the
> exons present in my ecs object (355320).
>
> I tried to set the extracols option in a different number of ways (putting
> as row.names the geneID or exonID, having only one column etc), but I always
> got this error:
>
> Error in data.frame(..., check.names = FALSE) :
> arguments imply differing number of rows: 355320, 396
>
> Where 355320 is the correct number of rows of both ecs exons and of my extra
> attributes data.frame. (I haven't figured out what has 396 rows yet).
>
> I checked the online reference manual and there is no mention of the
> extracols option. I also checked the code of the relevant function and I
> don't see it get used somewhere.
>
> Any ideas?
>
> Is there another way to add gene or exon attributes in the HTML?
>
> Thank you all,
>
> Best regards,
>
> Ioannis
>
More information about the Bioconductor
mailing list