[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