I think the easiest way to change the order of the panels is to change the order of the factor outside the plot call.
dat22 = augPred(fit.nlme3.10)
dat22$.groups = factor(dat22$.groups,
levels = c("Above:12", "Above:21", "Above:35", "Above:43", "Above:15",
"Above:23", "Above:32", "Above:41", "Above:13", "Above:24",
"Above:31", "Above:46", "Below:12", "Below:21", "Below:35",
"Below:43", "Below:15", "Below:23", "Below:32", "Below:41",
"Below:13", "Below:24", "Below:31", "Below:46"))
plot(dat22, layout = c(4,6), xlab="", ylab="", ylim=c(-200,2700))
To change the titles in each strip, you can set the factor.levels
using strip.custom
.
plot(dat22, layout = c(4,6), xlab="", ylab="", ylim=c(-200,2700),
strip = strip.custom(factor.levels = c("M", "M", "M", "M", "FP", "FP", "FP", "FP", "P", "P",
"P", "P","M", "M", "M", "M", "FP", "FP", "FP", "FP", "P",
"P", "P", "P")))
This replaces the original names, though, and I wasn't entirely certain if you wanted to augment the names or change them entirely.