Skip to content
Prev 303827 / 398503 Next

define subset argument for function lm as variable?

Hi

I want to do a series of linear models, and would like to define the input 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               # 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,                          # random effect in the model
              correlation = corAR1(form=~year|plant), #
              subset = dead==FALSE & recTreat==FALSE, #
              na.action = na.omit
              )
dat.lm <- lm(
             formula =  feff,              # fixed effect in the model
             data = dat,
             subset = dead==FALSE & recTreat==FALSE,
             na.action = na.omit
             )

Thanks,

Rainer