I have been using the
rnaseqDTU workflow and the function
testForDEU internally uses the
nbinormLRT. I was wondering what the ramifications are if I wanted to calculate fold change for example for specific timepoints.. Hypothetically could I use the
wald test instead of
LRT in this case to get FCs?
dxd <- estimateSizeFactors(dxd) dxd <- estimateDispersions(dxd, quiet=FALSE) dxd <- testForDEU(dxd, reducedModel=~sample + exon + replicate:exon)