Hi,

is it possible to compare this 7 factors (levels) at once?

dds$brain_areal = factor(dds$brain_areal, levels = c("primary and secondary corticies","limbic and association corticies", "archicortex", "thalamus and hypothalamus", "white matter structures", "cerebellar gray matter", "striatum"))

Normally R uses the first and the last factor for comparison. I want to know, if I can perform everything at once?

I want to create Networks for gene expression with wTO and then compare them with CoDiNA..

You could use the likelihood ratio test instead of the Wald test. The former tests a full model (that includes brain_areal) and a reduced model (that does not include brain_areal), to see if the full model fits the data better. See the vignette for more details.