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 infere
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(effect
Dear R People:
I have a data set with trt and resp:
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
And I did read.table, constructed an aov object and ran the summary on
the aov object. I also ran the TukeyHSD function. All is well so
far.
My question is: is there a way to
3 matches
Mail list logo