Hi,
I can not reproduce your problem.
Can you use ?dput, and try to create a self contained example?
Here is my attempt:
warpbreaks1 <- warpbreaks
warpbreaks1[1,1] <- NA
fm1 <- aov(breaks ~ tension, data = warpbreaks1)
summary(fm1)
TukeyHSD(fm1)
----------------Contact
Details:-------------------------------------------------------
Contact me: tal.gal...@gmail.com | 972-52-7275845
Read me: www.talgalili.com (Hebrew) | www.biostatistics.co.il (Hebrew) |
www.r-statistics.com (English)
----------------------------------------------------------------------------------------------
On Fri, Apr 20, 2012 at 4:21 PM, danipaty13 <danipat...@hotmail.com> wrote:
> 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.
>
> ______________________________________________
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
