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))