I don’t understand why when adding multiple fixed effects in different orders changes my output results. Shouldn’t the output be the same if I run each fixed effect individually as opposed to all at once since the models test them individually as well?

I kept on receiving error when I fit the lmer() function: “Error in initializePtr() : function ‘cholmod_factor_ldetA’ not provided by package ‘Matrix’”.
I have tried multiple times based recommendations from the previous post “Error in initializePtr() : function ‘cholmod_factor_ldetA’ not provided by package ‘Matrix'” when using lmer function in R

I couldn’t find any similar questions because all the questions I found deal with correcting for multiple comparisons within a single lmer model.

I would like to know how to get the overall wald test and its p-value in the glmer model. For example: credec ~ passexp + soctie + recexp + peer + soctie * recexp+ (1 | region), data = data, family = binomial

I have been using the modular lme4 functions to customize some of our problems and I noticed that the control argument seems to be ignored in the modular functions or I am just not placing it in the right place. Here a minimal example based on the rpubs from https://rpubs.com/bbolker/groupmembers