Dear all
I try to add a line, (vertical in my case) but at present time I am
only partially successful.
fm1 <- lme(Orthodont)
# standard plot
plot(augPred(fm1, level = 0:1, length.out = 2))
#plot with vertical but without points and fitted lines
plot(augPred(fm1, level = 0:1, length.out = 2),
panel=function(v,...) {
panel.abline(v=10)}
)
# plot with vertical but without fitted lines
plot(augPred(fm1, level = 0:1, length.out=2),
panel=function(x,y,...) {
panel.xyplot(x,y,...)
panel.abline(v=10)}
)
# plot with vertical and with all points (fitted lines are drawn as
points)
plot(augPred(fm1, level = 0:1),
panel=function(x,y,...) {
panel.xyplot(x,y,...)
panel.abline(v=10)}
)
I am probably somewhere close but I have no clue, which parameter I
shall modify to get measured points, fitted lines and vertical lines
in panels together.
Please help
Thank you
Best regards.
Petr Pikal
petr.pikal at precheza.cz