Additionally, I’m having issues running emmeans(): Error in X[ii, ii, drop = FALSE] %*% y[ii] : non-conformable arguments. Running without factorize = FALSE makes no difference, so I guess it has nothing to do with that. emmeans works fine with lmer output, as a control.
A bit off-topic, but I forgot to center my variables. Now I did and the p-values for treatment changed by a lot. Wasn’t that not suppose to happen? Or is this normal for sum.contrasts?