glmer (lme4), glmmPQL (MASS) and xtmepoisson (Stata)

Ben Bolker · 2010-01-04T19:43:22Z

lshtm.ac.uk> writes: > I'm trying to specify a generalized linear mixed model in R, > basically a Poisson model to describe monthly > series of counts in different regions. > My aim is to fit subject-specific curves, > modelling a non-linear trend for each region through random > effects for linear splines components (see Durban et al, > Stat Med 2005, or " Semiparametric regression" > by Ruppert et al, 2003). > > I use the command 'glmmPQL' in the MASS package and

Ben Bolker

Mon, Jan 4, 2010 11:43 AM

<Antonio.Gasparrini <at> lshtm.ac.uk> writes:

1. I would suggest continuing this conversation on the
r-sig-mixed-models at lists.r-project.org , which is specialized
for this kind of question.

  2. I don't think that a specific replacements for pdDiag/pdIdent
are on their way any time soon, but you can 
get a diagonal structure for the random effects: see p. 16 of

 (broken URL so gmane doesn't complain about long lines ...)

http://lme4.r-forge.r-project.org/slides/
  2009-07-21-Seewiesen/5LongitudinalD.pdf
 
for an example.

  Ben Bolker

glmer (lme4), glmmPQL (MASS) and xtmepoisson (Stata)

Thread (2 messages)