Hello,
Apologies if this has been asked before. I used the search function before posting and couldn't find anything on this forum.
I'm exploring mixed models for DoE in an effort to gain deeper understanding on how fixed and random effects contribute to a process output. However, when I try to manually use backwards elimination to simplify my fixed effect regression parameters, I notice that the confidence intervals for my parameter estimates (ex: my intercept, which is important for my purposes) get really wide, which is also reflected by wide and uneven confidence intervals in the profiler.
I would have thought that simpler models would free up more degrees of freedom for better variance / standard error measurements and thus make the confidence intervals smaller, but the opposite appears true. It's worth noting these effects are more pronounced when I increase the complexity of my random effects (ex: concatenating parameters as a random block vs performing variance component analysis on each random factor individually).
See below of an example DoE I have been trying to build and explore using simulated results (not shown):

X1 X2 and X3 would be fixed factors, while the rest are my example random factors.
Any help would be appreciated! So far it just seems like my best course of action would be to retain the full response surface model and not explore any stepwise elimination options at this point.
Best,
Matt