[R] a confidence interval from an ANOVA model

Eik Vettorazzi E.Vettorazzi at uke.de
Mon Jan 21 11:22:56 CET 2013


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)
summary(allEffects(ano))

hth.

Am 20.01.2013 03:11, schrieb Erin Hodgess:
> 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 put together a confidence interval
> from the aov object on a particular mean, please?  It seems like there
> should be.
> 
> Thanks,
> Have a great night/day,
> Sincerely,
> Erin
> 
> 


-- 
Eik Vettorazzi

Department of Medical Biometry and Epidemiology
University Medical Center Hamburg-Eppendorf

Martinistr. 52
20246 Hamburg

T ++49/40/7410-58243
F ++49/40/7410-57790
--
Pflichtangaben gemäß Gesetz über elektronische Handelsregister und Genossenschaftsregister sowie das Unternehmensregister (EHUG):

Universitätsklinikum Hamburg-Eppendorf; Körperschaft des öffentlichen Rechts; Gerichtsstand: Hamburg

Vorstandsmitglieder: Prof. Dr. Martin Zeitz (Vorsitzender), Dr. Alexander Kirstein, Joachim Prölß, Prof. Dr. Dr. Uwe Koch-Gromus



More information about the R-help mailing list