# NOT RUN {
data(datnorm)
cand.groups <- subbuild(datnorm, height, labvalue, region, smoker)
fitd <- cbind(datnorm, cand.groups)
subgr <- colnames(cand.groups)
res <- modav(resp = "y", trt = "treat", subgr = subgr, data = fitd,
covars = ~ x1 + x2, fitfunc = "lm")
summary(res)
# }
Run the code above in your browser using DataLab