Skip to content
Back to formatted view

Raw Message

Message-ID: <Pine.LNX.4.64.1102251433090.32269@orpheus.qimr.edu.au>
Date: 2011-02-25T04:35:37Z
From: David Duffy
Subject: Do I need random effects?
In-Reply-To: <4D671CE7.5090603@waikato.ac.nz>

On Fri, 25 Feb 2011, Murray Jorgensen wrote:

> Then I would like to do:
>
> m0 = lmer(mincr ~ icfac)
> anova(m0,m1)
>
> But lmer won't let me fit models with no random effects.

make a dummy variable eg
data(Pastes)
dummy <- rep(1,nrow(Pastes))
m1 <- lmer(strength ~ (1|batch), Pastes, REML=F)
m0 <- lmer(strength ~ (1|dummy), Pastes, REML=F)
anova(m1, m0)

-- 
| 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