[R-sig-ME] Help with glmmPQL using Scottish Lip Cancer data

[David Duffy]

On Mon, 7 Jun 2010, Corey Sparks wrote:

> glmmPQL(Observed~offset(log(Expected))+aff,family=poisson, data=scot,
> random=~1|District, correlation=spcor, niter=200, na.action="na.omit")

> nlminb problem, convergence error code = 1
> message = function evaluation limit reached without convergence (9)

When I added  control=lmeControl(opt="optim") to the above, it does 
converge.

FWIW, David Duffy.


-- 
| David Duffy (MBBS PhD)                                         ,-_|\
| email: davidD at qimr.edu.au  ph: INT+61+7+3362-0217 fax: -0101  /     *
| Epidemiology Unit, Queensland Institute of Medical Research   \_,-._/
| 300 Herston Rd, Brisbane, Queensland 4029, Australia  GPG 4D0B994A v