Skip to content

is this an ANOVA ?

1 message · Steven McKinney

#
No, the small p-values indicate that the associated estimate appears to be significantly different from zero.

You can use the package "multcomp" to do multiple comparisons.

     > require("multcomp")
     > lma <- aov(y ~ x, data = MD)
     > lmamc <- glht(lma, linfct = mcp(x = "Tukey"))
     > ci.lma <- confint(lmamc)
     > ci.lma

              Simultaneous Confidence Intervals

     Multiple Comparisons of Means: Tukey Contrasts
     
     
     Fit: aov(formula = y ~ x, data = MD)

     Quantile = 2.667
     95% family-wise confidence level
 

     Linear Hypotheses:
                Estimate lwr     upr    
     B - A == 0  5.0000   2.3330  7.6670
     C - A == 0 10.0000   7.3330 12.6670
     C - B == 0  5.0000   2.3330  7.6670

     > lmacld <- cld(lmamc)
     > plot(lmacld)
     > plot(ci.lma)
Not that I know of, though the multcomp tables and plots yield logically equivalent results and plots.
Writing a few lines of code to accomplish your graph is fairly straightforward.

HTH

Steven McKinney