search for: rectreat

...t arguments for lm() as variables. I managed easily to define the formula arguments in a variable, but I also would like to have the "subset" in a variable. My reasoning is, that I have the subset in the results object. So I wiould 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...

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 structure # dat.lme <- lme( fixed = feff,...