Fitting GLMs with Missing Data in Both Responses and Covariates

Fits generalized linear models (GLMs) when there is missing data in both the response and categorical covariates. The functions implement likelihood-based methods using the Expectation and Maximization (EM) algorithm and optionally apply Firth’s bias correction for improved inference. See Pradhan, Nychka, and Bandyopadhyay (2025) , Maiti and Pradhan (2009) , Maity, Pradhan, and Das (2019) for further methodological details.