HI, I guess there is a mistake in your code. You should have used "typ" instead of "abun" as "abun" is the dependent variable. summary(fm1 <- aov(breaks ~ wool + tension, data = warpbreaks)) myresults <- TukeyHSD(fm1, "tension", ordered = TRUE) library(agricolae)
HSD.test(fm1,"wool",group=TRUE) #Study: #HSD Test for breaks #Mean Square Error: 134.9578 #wool, means # breaks std.err replication #A 31.03704 3.050609 27 #B 25.25926 1.789963 27 #alpha: 0.05 ; Df Error: 50 #Critical Value of Studentized Range: 2.840532 #Honestly Significant Difference: 6.350628 #Means with the same letter are not significantly different. #Groups, Treatments and means #a A 31.037037037037 #a B 25.2592592592593 A.K. ----- Original Message ----- From: Landi <ent-ar...@gmx.de> To: r-help@r-project.org Cc: Sent: Friday, September 28, 2012 5:41 AM Subject: [R] Anova and tukey-grouping Hello, I am really new to R and it's still a challenge to me. Currently I'm working on my Master's Thesis. My supervisor works with SAS and is not familiar with R at all. I want to run an Anova, a tukey-test and as a result I want to have the tukey-grouping ( something like A - AB - B) I came across the HSD.test in the agricolae-package, but... unfortunately I do not get an output (like here in the answer http://stats.stackexchange.com/questions/31547/how-to-obtain-the-results-of-a-tukey-hsd-post-hoc-test-in-a-table-showing-groupe ) I did it like this: ###### ANOVA anova.typabunmit<-aov(ds.typabunmit$abun ~ ds.typabunmit$typ) summary(anova.typabunmit) summary.lm(anova.typabunmit) ###### post HOC tukey.typabunmit<-TukeyHSD(anova.typabunmit) tukey.typabunmit ###### HSD HSD.test(anova.typabunmit, "abun", group=TRUE) and the ONLY output is this: Name: abun ds.typabunmit$typ I would be very pleased about some ides..:! -- View this message in context: http://r.789695.n4.nabble.com/Anova-and-tukey-grouping-tp4644485.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. ______________________________________________ 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.
