[BioC] Using edgeR with no replicates, estimating dispersion and the pseudo.alt slot
mark.robinson at imls.uzh.ch
Fri Oct 26 15:27:11 CEST 2012
> useful measures of normalised expression
Would the cpm() function suit your needs? It's not identical to $pseudo.alt, but should be similar.
On 25.10.2012, at 17:37, WATSON Mick wrote:
> We have a collaborator with RNA-Seq data with no replicates - and yes, I understand the limitations of this.
> In the egdeR manual, one solution suggested to estimating common dispersion is choosing a value based on experience e.g. "d$common.dispersion <- 0.4". However, this has the disadvantage of not filling up the pseudo.alt slot of the object, which are useful measures of normalised expression. The estimateTagWiseDisp() function then fails.
> It's fairly trivial to edit the estimateCommonDisp() function to accept a user provided common dispersion value, but I am not sure this is entirely valid.
> The question being how to populate the pseudo.alt slot when you are estimating the common dispersion parameter yourself, and whether it is even valid to do so.
> The University of Edinburgh is a charitable body, registered in
> Scotland, with registration number SC005336.
> Bioconductor mailing list
> Bioconductor at r-project.org
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
Prof. Dr. Mark Robinson
Institute of Molecular Life Sciences
University of Zurich
v: +41 44 635 4848
f: +41 44 635 6898
e: mark.robinson at imls.uzh.ch
More information about the Bioconductor