A Skew-normal copula-driven GLMM

This paper presents a method for fitting a copula‐driven generalized linear mixed models. For added flexibility, the skew‐normal copula is adopted for fitting. The correlation matrix of the skew‐normal copula is used to capture the dependence structure within units, while the fixed and random effects coefficients are estimated through the mean of the copula. For estimation, a Monte Carlo expectation–maximization algorithm is developed. Simulations are shown alongside a real data example from the Framingham Heart Study.