Thank you. Currently my input is RMA normalised microarray data. So I guess I can analyse this directly.
Does GAGE package in R require log-transformed normalized expression values or "just" normalized expression values as an input?
2 answers
What are you using GAGE for (RNA-seq or microarray)?
For RNA-seq you can use log2 normalized data:
cnts=raw_gene_count_data
sel.rn=rowSums(cnts) != 0
cnts=cnts[sel.rn,]
libsizes=colSums(cnts)
size.factor=libsizes/exp(mean(log(libsizes)))
cnts.norm=t(t(cnts)/size.factor)
cnts.norm=log2(cnts.norm+8)
For microarray data, they use RMA/FARMS normalization in the vignette.
Why does the last line add 8 to cnts.norm? I understand something must be added to prevent -infs in the log2 transformed dataset, but why specifically 8?
You may use RMA, FARMS or other normalization methods. Array normalization usually does log(2) transformation too like in RMA and FARMS. If not, it is always advisable to do log2 or log transformation on array or RNA-Seq data for differential expression or pathway analysis.
Log in to answer this question.