[R] Tukey HSD not working

PGZC polizaragoza at hotmail.com
Sun Dec 12 19:21:59 CET 2010


Drug	US1	US2	Aptecha
Celebrex	235.54	269.99	121.02
Detrol LA	157.99	190.99	55.3
Flomax	166.00	190.99	93.45

Lipitor	174.99	200.99	137.7

Novaldex	108.6	129.99	22.48

Norvasc	186.66	203.99	161.93
Plavix	107.99	106.99	64.53
Prevacid	117.39	134.99	59.83
Prilosec	115.99	126.99	57.75
Zyrtec	181.1	200.99	58.79

US1=c(235.54,157.99,166,174.99,108.60,186.66,107.99,117.39,115.99,181.10)
US2=c(269.99,190.99,190.99,200.99,129.99,203.99,106.99,134.99,126.99,200.99)
Aptecha=c(121.02,55.30,93.45,137.70,22.48,161.93,64.53,59.83,57.75,58.79)
y=c(US1,US2,Aptecha)
n=rep(10,3)
group=rep(1:3,n)
data = data.frame(y = y, group = factor(group))
fit = lm(y ~ group, data)
anova(fit) 
summary (fit)

-- 
View this message in context: http://r.789695.n4.nabble.com/Tukey-HSD-not-working-tp3084505p3084505.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list