Displaying 2 results from an estimated 2 matches for "pharv".
Did you mean:
pharm
2012 Aug 22
3
Question concerning anova()
...inear mixed effect models, derived from updating the original one. To
compare these, I want to use anova(). I therefore do the following (not reproducible - just to
illustration purpose!):
dat <- loadSPECIES(SPECIES)
subs <- expression(dead==FALSE & recTreat==FALSE)
feff <- noBefore~pHarv*year # fixed effect in the model
reff <- ~year|plant # random effect in the model, where year is the
corr <- corAR1(form=~year|plant) # describing the within-group correlation structure
#
dat.lme <- lme(
fixed = feff, # fixed effect...
2012 Aug 21
2
define subset argument for function lm as variable?
...ad==FALSE & recTreat==FALSE
which obviously does not work as the expression is evaluated immediately. Is is it possible to do
what I want to do here, or do I have to go back to use
dat <- subset(dat, dead==FALSE & recTreat==FALSE)
?
dat <- loadSPECIES(SPECIES)
feff <- height~pHarv*year # fixed effect in the model
reff <- ~year|plant # random effect in the model, where year is the
dat.lme <- lme(
fixed = feff, # fixed effect in the model
data = dat,
random = reff,...