The conf.reg() function calculates confidence regions for contrasts in multivariate factorial designs. In the two-dimensional case, confidence ellipsoids can be plotted via the generic plot() function.
plot()
Useful links