lattice_segmented.R

## Segmented regression plots in lattice
## Example for view purposes only
head(iris)
library(xyplot)
xyplot(Sepal.Length ~ Sepal.Width, iris,
       panel = function(x, y, ...){
    panel.xyplot(x, y, ...)
    le <- lm(Sepal.Length ~ Sepal.Width, iris,
         subset = Sepal.Width <= 3)
    llines(iris$Sepal.Width[iris$Sepal.Width <= 3],
          predict(le))
    ld <- lm(Sepal.Length ~ Sepal.Width, iris,
             subset = Sepal.Width > 3)
    llines(iris$Sepal.Width[iris$Sepal.Width > 3],
          predict(ld), col = 2)
})