[R] a confidence interval from an ANOVA model

Yvonnick Noel yvonnick.noel at uhb.fr
Wed Jan 23 08:12:49 CET 2013

Hello Erin,

You might want to compute *credibility* intervals instead. They have the 
true meaning of intervals within which the *unknown parameter* (not the 
sample statistic) has a fixed probability to lie (confidence intervals 
do not).

The AtelieR package, especially in its "Bayesian inference on several 
means" section, lets you input the data summary (means, standard 
deviation and counts) and automatically provides these intervals on the 
unknown means.

The approach implemented is the one described in:

Neath, A. & Cavanauigh, J. (2006). A bayesian approach to the multiple 
comparisons problem. Journal of Data Science, 4, 131-146.

Besides, it will give you the best constrained model (in terms of equal 
means) for your data, and the probability that this model is the true 
one. It is an elegant solution to the multiple comparisons problem.


Yvonnick Noel
University of Brittany, Rennes

More information about the R-help mailing list