search for: loadspeci

...like to add a line like: subs <- dead==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...

Hi I am comparing four different linear 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 struc...