[R] mixed model nested ANOVA (part two)
Mark Difford
mark_difford at yahoo.co.uk
Sun Feb 24 22:25:28 CET 2008
Hi Stephen,
Slip of the dactylus: lm() does not, of course, take a fixed=arg. So you
need
To recap:
mod.rand <- lme(fixed=y ~ x, random=~x|Site, data=...)
mod,fix <- lm(y ~ x, data=...) ## or
##mod,fix <- lm(formula=y ~ x, data=...)
Bye.
Mark Difford wrote:
>
> Hi Stephen,
>
>>> Also i have read in Quinn and Keough 2002, design and analysis of
>>> experiments for
>>> biologists, that a variance component analysis should only be conducted
>>> after a rejection
>>> of the null hypothesis of no variance at that level.
>
> Once again the caveat: there are experts on this list who really know
> about this stuff, and I am not one of them. Your general strategy would
> be to set up two models with the same fixed effects, one of which doesn't
> have random effects. You then test the two models using
> anova(mod.withRandom, modWithoutRandom).
>
> I haven't tried this using lmer/2(), but with lme() you do this by fitting
> your fixed+random effects model using lme() and your fixed-only effects
> model using lm(). If you are using weights to model heteroskedasticity,
> then it's better to use gls(), as this will accept the same weights
> argument as the call to lme().
>
> Then you simply do anova(lme.model, lm/gls.model). This tells you about
> the significance of your random effects, i.e. whether you need a
> random-effects component.
>
> To recap:
> mod.rand <- lme(fixed=y ~ x, random=~x|Site, data=...)
> mod,fix <- lm(fixed=y ~ x, data=...)
>
> anova(mod.rand, mod.fix)
>
> HTH, Mark.
>
>
> Stephen Cole-2 wrote:
>>
>> First of all thank you for the responses. I appreciate the
>> suggestions i have received thus far.
>>
>> Just to reiterate
>>
>> I am trying to analyze a data set that has been collected from a
>> hierarchical sampling design. The model should be a mixed model
>> nested ANOVA. The purpose of my study is to analyze the variability
>> at each spatial scale in my design (random factors, variance
>> components), and say something about the variability between regions
>> (fixed factor, contrast of means). The data is as follows;
>>
>> region (fixed)
>> Location (random)
>> Site(random)
>>
>> site nested in location nested in region.
>>
>> Also i have read in Quinn and Keough 2002, design and analysis of
>> experiments for biologists, that a variance component analysis should
>> only be conducted after a rejection of the null hypothesis of no
>> variance at that level.
>>
>> I have tried to implement
>> mod1<-lmer(density ~ 1 + (1|site) + (1|location) + (1|region))
>>
>> However, as i understand it, this treats all my factors as random.
>> Plus I do not know how to extract SS or MS from this model.
>>
>> anova(mod1) gives me
>> Analysis of Variance Table
>> Df Sum Sq Mean Sq
>>
>> and summary(mod1) gives me
>> Linear mixed-effects model fit by REML
>> Formula: density ~ 1 + (1 | site) + (1 | location) + (1 | region)
>> AIC BIC logLik MLdeviance REMLdeviance
>> 15658 15678 -7825 15662 15650
>> Random effects:
>> Groups Name Variance Std.Dev.
>> site (Intercept) 22191 148.97
>> location (Intercept) 33544 183.15
>> region (Intercept) 41412 203.50
>> Residual 696189 834.38
>> number of obs: 960, groups: site, 4; location, 4; region, 3
>>
>> Fixed effects:
>> Estimate Std. Error t value
>> (Intercept) 261.3 168.7 1.549
>>
>> from what i understand the variance in the penultimate column are my
>> variance components. But how do i conduct my significance test?
>>
>> I have also tried
>> mod1<-lmer(density ~ region + (1|site) + (1|location))
>>
>> Which i think is the correct mixed model for my design. However once
>> again i do not know how to evaluate significance for the random
>> factors.
>>
>> Thank-you again for any additional advice i receive
>>
>> Stephen Cole
>>
>> ______________________________________________
>> R-help at r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide
>> http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>>
>>
>
>
--
View this message in context: http://www.nabble.com/mixed-model-nested-ANOVA-%28part-two%29-tp15665478p15669608.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list