\begin{Schunk}
\begin{Sinput}
 library(rockchalk)
 example(predictOMatic)
 ## Or create a fresh example for yourself with:
 dat <- genCorrelatedData2(1000, 
      means = c(100,200,100, 40), 
      sds = c(10, 20, 20, 5), 
      rho = 0.2, 
      beta = c(0.1, 0.3, 0.3, 0.1, -0.14, 0, 0.025), 
      stde = 140)
 dat$xcat1 <- factor(sample(c("No","Yes","Yes"), 1000, replace = TRUE))
 m1 <- lm(y ~ x1 + x2 + x3 + x4 + xcat1, data = dat)
 summary(m1)
 predictOMatic(m1)
 predictOMatic(m1, divider = "std.dev.")
 predictOMatic(m1, divider = "seq", n = 7)
\end{Sinput}
\end{Schunk}
