Use the WGCNA package.
To install, use the command
install.packages('WGCNA')
This package requires the installation of the impute package that can be installed using
source("http://bioconductor.org/biocLite.R") biocLite("impute")
An example:
The code used
verboseBarplot(height,group,numberStandardErrors=2,AnovaTest=TRUE, main="2 SE, Anova", ylim=c(0, 2.5))