[R] Problem with Tukey test
danipaty13
danipaty13 at hotmail.com
Fri Apr 20 15:21:32 CEST 2012
I'm new using R im trying to do a tukey test, but when i see the results the
p value results in NA im guessing its because i have missing values but im
not sure how to fix it
AnovaModel.2 <- aov(area ~ trat, data=apilados)
> summary(AnovaModel.2)
Df Sum Sq Mean Sq F value Pr(>F)
trat 4 11847 2961.76 9.9905 1.500e-06 ***
Residuals 76 22531 296.46
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
14 observations deleted due to missingness
> numSummary(apilados$area , groups=apilados$trat, statistics=c("mean",
> "sd"))
mean sd n NA
X1 79.29625 15.06667 16 3
X2 77.43400 10.64849 15 4
X3 71.19895 15.85500 19 0
X4 77.80938 24.67359 16 3
X5 45.97200 16.65112 15 4
> TukeyHSD(AnovaModel.2)
Tukey multiple comparisons of means
95% family-wise confidence level
Fit: aov(formula = area ~ trat, data = apilados)
$trat
diff lwr upr p adj
X2-X1 -1.862250 NA NA NA
X3-X1 -8.097303 NA NA NA
X4-X1 -1.486875 NA NA NA
X5-X1 -33.324250 NA NA NA
X3-X2 -6.235053 NA NA NA
X4-X2 0.375375 NA NA NA
X5-X2 -31.462000 NA NA NA
X4-X3 6.610428 NA NA NA
X5-X3 -25.226947 NA NA NA
X5-X4 -31.837375 NA NA NA
--
View this message in context: http://r.789695.n4.nabble.com/Problem-with-Tukey-test-tp4573945p4573945.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list