# NOT RUN {
data(aidscd4)
lmd.bl <- bcreg(cd4.bl ~ 1, data = aidscd4[aidscd4$weekc == 8, ])$lambda
aidscd4$cd4.bl.tr <- (aidscd4$cd4.bl ^ lmd.bl - 1) / lmd.bl
resar <- bcmmrm(outcome = cd4, group = treatment, data = aidscd4,
time = weekc, id = id, covv = c("cd4.bl.tr", "sex"),
cfactor = c(0, 1), structure = "AR(1)", glabel =
c("Zid/Did", "Zid+Zal", "Zid+Did", "Zid+Did+Nev"))
boxplot(resar, xlab = "Treatment", col = 1:4)
# }
Run the code above in your browser using DataLab