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.
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.Galili@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 <danipaty13@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]]
http://r.789695.n4.nabble.com/file/n4574990/invgrosor.csv invgrosor.csv Thanks for your answer well, I uploaded a file with mi data, hope you can solve my problem since i'm using R for my thesis -- View this message in context: http://r.789695.n4.nabble.com/Problem-with-Tukey-test-tp4573945p4574990.html Sent from the R help mailing list archive at Nabble.com.
HELLO WELCOME TO R here is your code: 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() If your original data.frame or your raw data have "NA's" you must to write this part AnovaModel.2 <- aov(area ~ trat, na.action=na.omit, , data=apilados) I hope that work in this way suerte!!! Alex Ponce -- View this message in context: http://r.789695.n4.nabble.com/Problem-with-Tukey-test-tp4573945p4575738.html Sent from the R help mailing list archive at Nabble.com.
thank you soooo much i just tried the code you posted and worked perfectly -- View this message in context: http://r.789695.n4.nabble.com/Problem-with-Tukey-test-tp4573945p4575943.html Sent from the R help mailing list archive at Nabble.com.