BMlmSE() with fixed effects?

[yukitakahashi1]

Hi! I'm now preparing my experiment and referring to this nice SOP.
Could we use BMlmSE() function with fixed effects? It only accepts lm() objects, so with factor() we get very sparse X and BMlmSE() stops working.

[linstonwin]

Thanks for your kind comment! I'm not actively involved in the SOP anymore and also haven't been following the literature on this issue recently, but you might want to take a look at the clubSandwich package and the paper by James Pustejovsky and Beth Tipton:

https://www.jepusto.com/software/clubsandwich/

All best wishes,
Winston

[yukitakahashi1]

Thank you very much Winston - I will check them!
Sincerely,
Yuki