> Why not use the subset option? Something like: > > lm(diff ~ Age + Race, data=data, subset=data$Meno=="PRE") > should do the trick, and be much easier to read! "data$" could be omitted, simply lm(diff ~ Age + Race, data=data, subset=Meno=="PRE")