Entering edit mode
4.3 years ago
anamaria
▴
220
Hello,
I am planing to using ssizeRNA to calculate power for particular number of subjects. I can see from docs: https://cran.r-project.org/web/packages/ssizeRNA/ssizeRNA.pdf
I do have my expression matrix (rows genes, columns subjects) and I also have limma's topTable() top ranked genes document. Is there is a way to generate from this average gene counts an how to do that? I am guessing that dispersion (also required by ssizeRNA) I can calculated using: https://www.rdocumentation.org/packages/edgeR/versions/3.14.0/topics/estimateDisp