R help - Feb 2004 - Offset in GLMM

Dear R-list,

I try to adjust GLMM on incidence cancer data. Without random effect, in GLM
the command is, for example with sex effect,
glm(Observed~sex+offset(log(Expected)),family=poisson) because the observed
are Poisson distribued with parameter Expected*incidence rate. But know I
want to introduce random effect (for example spatial effect) and it seems to
me that the "offset" does not work in
glmmPQL(Observed~sex+offset(log(Expected)),random=~1|ED,...,family=poisson)?
I mean it does not work in lme.

what is the problem with my offset? thank you in advance, with best regards,

Erik.

===========================================Erik-André SAULEAU, MD, PhD

SEAIM
Hôpital du Hasenrain
BP 1070
87, Avenue de Altkirch
68051 MULHOUSE Cédex

Tel: 03-89-64-79-95
Fax: 03-89-64-79-71
Mél: sauleauea@ch-mulhouse.fr
Web: http://www.ch-mulhouse.fr

-----------------------------------

Registre des Cancers du Haut-Rhin
9, Rue du Dr Mangeney
BP1370
68070 MULHOUSE

Web: www.arer68.org
===========================================




************************************************************************************
Afin d'éviter toute propagation de virus informatique, et en complément des
dispositifs en place, ce message (et ses pièces jointes s'il y en a) a été
automatiquement analysé par un antivirus de messagerie.
***********************************************************************************


	[[alternative HTML version deleted]]

SAULEAU Erik-Andr? <SAULEAUEA at ch-mulhouse.fr> writes:
> Dear R-list,
> 
> I try to adjust GLMM on incidence cancer data. Without random
> effect, in GLM the command is, for example with sex effect,
> glm(Observed~sex+offset(log(Expected)),family=poisson) because the
> observed are Poisson distribued with parameter Expected*incidence
> rate. But know I want to introduce random effect (for example
> spatial effect) and it seems to me that the "offset" does not
work
> in
>
glmmPQL(Observed~sex+offset(log(Expected)),random=~1|ED,...,family=poisson)?
> I mean it does not work in lme.
> 
> what is the problem with my offset? thank you in advance, with best
> regards,
Did you try function GLMM from package lme4?  I believe that GLMM
allows an offset.