No testable fixed effect were included in the model

[yasamanrezvani]

Hi,

I have bulk RNA-seq data with one fixed effect (Condiiton) and 2 random effects (patients and batch Info). I am running dream to perform DE analysis.

param <- SnowParam(4, "SOCK", progressbar = TRUE)
formula <- ~ Condition + (1 | patient) + (1 + batch)
voomObj <- voomWithDreamWeights(dgeObject, formula, metadata, BPARAM=param)
fit <- dream(voomObj, formula, metadata)

However, running the dream function returns a warning No testable fixed effects were included in the model. Running topTable() will fail

How can I solve this issue? Any help is much appreciated!