\begin{Schunk}
\begin{Sinput}
 pm1 <- predictOMatic(m2, divider="std.dev.", n=5 )
 print(pm1, digits=3)
\end{Sinput}
\begin{Soutput}
$x1
             fit    x1       x2         x3           x4 x5
(m-2sd) 15.90005 30.14 50.46198 0.01624596 -0.005435391  A
(m-sd)  17.95194 40.22 50.46198 0.01624596 -0.005435391  A
(m)     20.00383 50.30 50.46198 0.01624596 -0.005435391  A
(m+sd)  22.05572 60.38 50.46198 0.01624596 -0.005435391  A
(m+2sd) 24.10761 70.46 50.46198 0.01624596 -0.005435391  A

$x2
             fit    x2       x1         x3           x4 x5
(m-2sd) 16.13013 30.48 50.30392 0.01624596 -0.005435391  A
(m-sd)  18.06718 40.47 50.30392 0.01624596 -0.005435391  A
(m)     20.00424 50.46 50.30392 0.01624596 -0.005435391  A
(m+sd)  21.94130 60.45 50.30392 0.01624596 -0.005435391  A
(m+2sd) 23.87836 70.44 50.30392 0.01624596 -0.005435391  A

$x3
             fit    x3       x1       x2           x4 x5
(m-2sd) 20.28135 -1.88 50.30392 50.46198 -0.005435391  A
(m-sd)  20.14272 -0.93 50.30392 50.46198 -0.005435391  A
(m)     20.00408  0.02 50.30392 50.46198 -0.005435391  A
(m+sd)  19.86544  0.97 50.30392 50.46198 -0.005435391  A
(m+2sd) 19.72681  1.92 50.30392 50.46198 -0.005435391  A

$x4
             fit    x4       x1       x2         x3 x5
(m-2sd) 20.36103 -2.07 50.30392 50.46198 0.01624596  A
(m-sd)  20.18322 -1.04 50.30392 50.46198 0.01624596  A
(m)     20.00542 -0.01 50.30392 50.46198 0.01624596  A
(m+sd)  19.82761  1.02 50.30392 50.46198 0.01624596  A
(m+2sd) 19.64980  2.05 50.30392 50.46198 0.01624596  A

$x5
             fit x5       x1       x2         x3           x4
A (70%) 20.00463  A 50.30392 50.46198 0.01624596 -0.005435391
B (30%) 20.07446  B 50.30392 50.46198 0.01624596 -0.005435391

attr(,"flnames")
[1] "x1" "x2" "x3" "x4" "x5"
\end{Soutput}
\end{Schunk}
