Hello,

How could I use the lfcshrink() function when I need to examine the condition effect on genotype III?

results(dds, contrast=list( c("condition_B_vs_A","genotypeIII.conditionB") ))

I tried the below code but I got an error.

lfcShrink(dds, contrast = list(c("condition_B_vs_A", "genotypeIII.conditionB")))

#Error in lfcShrink(dds, contrast = list(c("condition_B_vs_A", "genotypeIII.conditionB")),  : 
 #type='apeglm' shrinkage only for use with 'coef'

First, read over the section on the vignette that addresses this, in the Extended section on shrinkage estimators. Easiest would be to use ashr.