Re: [R] a confidence interval from an ANOVA model

2013-01-22 Thread Yvonnick Noel
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

Re: [R] a confidence interval from an ANOVA model

2013-01-21 Thread Eik Vettorazzi
Hi Erin, how about this: dta-read.table(textConnection(trt resp A 23 A 28 A 37 A 30 B 37 B 44 B 31 B 35 C 42 C 47 C 52 C 38),header=T) summary(ano-aov(resp~trt,dta)) ano0-update(ano,.~.+0) confint(ano0) # and the corresponding means are coef(ano0) # or use the effects package library(effects)