glmmFEL_pl function

Pseudo-likelihood engine for RSPL/MSPL (Wolfinger-style, simplified R = I)