[R] StepAIC and boxplot

[David martin]

Hello,
I have done a stepwise analysis to determine the best model fitting my data.

data <- stepAIC( aov (Group)~ . , data=mydata)

data
 > summary(data)
                  Df  Sum Sq Mean Sq F value    Pr(>F)
`a`  1  1.9829 1.98290  11.176 0.0011824 **
`b`   1  2.6606 2.66064  14.996 0.0001967 ***
`c`   1  2.0339 2.03394  11.464 0.0010295 **
Residuals        96 17.0325 0.17742
---

 > data$anova
Stepwise Model Path
Analysis of Deviance Table

Initial Model:
as.integer(Group) ~ `a` + `b` + `c` +
     `d` + `e` + `f`

Final Model:
as.integer(Group) ~ `a` + `b` + `c`


                    Step Df   Deviance Resid. Df Resid. Dev       AIC
1                                            93   16.83679 -164.1604
2     - `d`  1 0.01387782        94   16.85067 -166.0780
3 - `e`  1 0.09713881        95   16.94781 -167.5032
4   - `f`  1 0.08471704        96   17.03252 -169.0046
 >

I would like to do a boxplot of the data so that i visualize the model. 
How can i boxplot the model ('a'+'b'+'c') directly ??

thanks,

david