# NOT RUN {
require(lattice)
xyplot(circumference ~ age, Orange, groups = Tree, type = c("g", "b"),
auto.key = list(space = "right", lines = TRUE), aspect = "xy",
xlab = "Age (days since 1968/12/31)", ylab = "Circumference (mm)")
# }
# NOT RUN {
m1 <- nlmer(circumference ~ SSlogis(age, Asym, xmid, scal) ~ Asym|Tree,
Orange, verbose = TRUE,
start = c(Asym = 190, xmid = 730, scal = 350))
.Call("mer_optimize", m1, 1L, 1L, PACKAGE = "lme4")
print(m1)
ranef(m1)
# }
Run the code above in your browser using DataLab