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.

HTH,

Yvonnick Noel
University of Brittany, Rennes
France

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to