I am using GLLAMM to fit a relatively simple mixed effects model with a
random intercept for individual (gllamm tracks2 sex distperday2,
i(otterid) link(logit) family(bin) robust adapt). I want to use a model
selection strategy based on AIC, but I am struggling with the
calculation of p as presented in Vaida and Blanchard (2005) or more
recent papers (e.g., Liang et al. 2008). I am hoping for guidance on
the calculation or a do file to automate the process.