getting residuals from dream

Question

[RNAseqer]

I would like to get the residuals for the variables I used in linear mixed effect model I used in dream. The variable I am interested in is continuous if that matters. I'm unsure how to do this or if it is even possible to do with dream? Here is where I am so far:

cpm.tmm.logno <- cpm(tmm.normfact.out, log=FALSE, normalized.lib.sizes = TRUE) 

form.dream <- ~ Disease + BMI  + (1|Batch) + RIN + (1|Sex) 

vobj <- voomWithDreamWeights(cpm.tmm.logno,   form.dream,   demos)

fitmm <- dream(vobj, form.dream,  demos)

toptable.out <- topTable(fitmm, coef=2, number=nrow(y), adjust.method = 'BH',p.value = 1.0)

If I wanted to then see the residuals for BMI, is there a way I could do that from this point?