Message-ID: <56EA19DB.6000606@gmail.com>
Date: 2016-03-17T02:43:39Z
From: Ben Bolker
Subject: Convergence warning message
In-Reply-To: <CALrjt7-C8EdYw3jqv0MC22Fx48C4DGqcUXjnFxw+AwkjHQs5DQ@mail.gmail.com>
To answer this particular question:
The importance of accuracy in the approximation of the integral over
the random effects (for which glmmPQL < Laplace < adaptive Gauss-Hermite
quadrature) depends on the degree to which the sampling distribution of
the conditional modes (approx "BLUPs", i.e. the values associated with
deviations of particular grouping factors from the population average)
are Gaussian. In particular, if you have Bernoulli responses and small
numbers of samples per group (which seems to be the case here -- 1 to 4
observations per individual if I read your data correctly) then yes,
this will be a concern ...
On 16-03-16 04:33 PM, Christopher David Desjardins wrote:
> From what I've read online, glmmPQL is inappropriate with Bernoulli
>> >>trials.
>> >>Is that correct?