Skip to content
Back to formatted view

Raw Message

Message-ID: <000a01c31de0$e9ee6360$3d9aee82@uu.se.vaxtbio>
Date: 2003-05-19T08:30:33Z
From: Mikael Niva
Subject: Syntax for random effect in glmmPQL

Dear R-listers

I wonder if someone can help me with the syntax for the random effect in
glmmPQL()? I have a data set with a response variable "y" (counts), two
dependent variables: "treat" (4 levels) and "site" (2 levels). The
latter, I want to use as a random variable. How do I specify this in the
function?

Is it like this: glmmPQL(y~treat,random=~1|site,family=poisson,
data=mydata)? What does the ~1 do?

Thanks in advance, Mikael Niva

********************************************
Mikael Niva
Avd. f?r V?xtekologi, Dept. of Plant Ecology 
EBC
Uppsala Universitet