search for: tukeyhsd

Hello R People: Here's the Saturday night goofy question. I would like to see the code for TukeyHSD function and I tried the following: > getAnywhere("TukeyHSD") A single object matching ?TukeyHSD? was found It was found in the following places package:stats namespace:stats with value function (x, which, ordered = FALSE, conf.level = 0.95, ...) UseMethod("TukeyHSD") &...

Dear r-help members. I would like to compare species numbers of moths between eight different forests (each sampled for six nights). I would like to do a nested anova to compare species numbers between forests and nights. For more site specific details I wanted to do a Tukey test (TukeyHSD). Unfortunately this test is not working (message: "nicht anwendbare Methode f?r 'TukeyHSD' auf Objekt der Klasse "c('anova', 'data.frame" ) /(message: "method not appropriate for TukeyHSD for objects of the class "anova", "data.frame")....

Dear R-developeRs, Attached follows a patch against svn 34959 that adds the printing of p-values to the TukeyHSD.aov function in stats package. I also updated the corresponding documentation file and added a 'see also' reference to the simint function of the multcomp package. As it was already brought up in a previous thread [1] in R-help, one can obtain the adjusted p-values using the multco...

...an prove I've done this before, but I recently installed Rexcel (and it was easiest to reinstall R and some other bits to make it work) and now its no longer working. Before I would do an ANOVA and a tukey post-hoc like this: >data1.aov=aov(result~factor1*factor2, data=data1) then... >TukeyHSD(summary(data1.aov)) and it would give me a nice tukey table of all the pairwise comparisons, now though it gives me: >in UseMethod("TukeyHSD") : no applicable method for "TukeyHSD" Has something changed? Does TukeyHSD no longer accept aov results? Is there a package I...

Dear R-Helpers, I was calling the TukeyHSD function and not getting confidence intervals or p-values. It turns out this was caused by missing data and the fact that I had previously turned on R Commander (Rcmdr). John Fox knew that Rcmdr sets na.action to na.exclude, which causes the problem. If you have this problem, you can either exit Rc...

Hello, I am trying to determine the most appropriate way to run post-hoc comparisons on my lme model. I had originally planned to use Tukey HSD method as I am interested in all possible comparisons between my treatment levels. TukeyHSD, however, does not work with lme. The only other code that I was able to find, and which also seems to be widely used, is glht specified with Tukey: summary(glht(model, linfct=mcp(Treatment="Tukey"))) Out of curiosity, I ran TukeyHSD and the glht code for a simple ANOVA and found...

Hello, When plotting the results of a TukeyHSD multiple comparisons procedure with an ANOVA (lm) object, an extra line appears in the confidence intervals that contain 0. For example (this is straight from the TukeyHSD helpfile): > summary(fm1 <- aov(breaks ~ wool + tension, data = warpbreaks)) > TukeyHSD(fm1, "tension", or...

I would like to do a one way anova and then a tukeyHSD. I have three vectors A,B and C. In a previous help message, I was told to do the following for the anova: y = c(A,B,C) group = factor(rep(a:3,c(7,9,13))) #provided there a 7 elements in A,9 in B and 13 in C and then anova(lm(y~group)) Looking at the tukeyHSD method it looks like it wants...

pairwise.t.test is returning NAs when one of the samples only has one entry, while TukeyHSD returns results (maybe not trustworthy or believable, but results). I stumbled on this because I did not realize one of my samples only had one entry while most of the others had several hundred, so I realize this is not a desirable situation. I'm really just curious about the difference be...

Hello everyone, So I ran an anova with aov and then I want to run post-hoc comparisons but keep receiving this message : > no applicable method for "TukeyHSD" Here is my code: > d<-read.table("d.txt") > d > Obs subj Hand Gaze RT > 1 1 s1 1 1 401.4 > 2 2 s2 1 1 363.3...... > summary(ano <- aov(RT~(Hand*Gaze)+Error(subj/(Hand*Gaze)),data=d )) This seems to work fine, but then I use &g...

Why do I see this error? > library(stats) > require(stats) [1] TRUE > > tHSD <- TukeyHSD(aov) Error in TukeyHSD(aov) : no applicable method for "TukeyHSD" In case it helps: > aov Call: aov(formula = roi ~ (Cue * Hemisphere) + Error(Subject/(Cue * Hemisphere)), data = roiDataframe) Grand Mean: 8.195069 Stratum 1: Subject Terms: Residuals Sum of S...

Dear R community, I'm trying to understand this behavior of TukeyHSD. My goal is to obtain defensible, labelled multiple comparisons of an interaction term. Firstly, if I plot the TukeyHSD from the model that calculates its own interactions, then the y-axis labels appear to be reflected on their median when compared to the text output of the TukeyHSD statement....

...shouldn't matter. Any suggestions? Thanks Amy > summary(aov(EtoH~treat_code, mydata)) Df Sum Sq Mean Sq F value Pr(>F) treat_code 1 16.44 16.44 11.027 0.001014 ** Residuals 287 427.91 1.49 --- Signif. codes: 0 ?***? 0.001 ?**? 0.01 ?*? 0.05 ?.? 0.1 ? ? 1 > TukeyHSD(aov(EtoH~treat_code, mydata)) Error in TukeyHSD.aov(aov(EtoH ~ treat_code, mydata)) : no factors in the fitted model In addition: Warning message: In replications(paste("~", xx), data = mf) : non-factors ignored: treat_code -- Amy Freitag Ph.D. student, Marine Science and Conservat...

I have a question about TukeyHSD and the glht function because I'm getting different answers when a covariate is included in model for ANCOVA.? I'm using the cabbages dataset in the 'MASS' package for repeatability.? If I include HeadWt as a covariate, then I get different answers when performing multiple compa...

Dear R-help, I have been trying to use the TukeyHSDs function in the "pgirmess" package to quickly extract all significant pairwise comparisons in an aov object. However, it seems that this function isn't working as intended when only the two last populations means being tested are significant. An example of this can be seen below:...

...ctor2 + Error(Subject/(Factor1*Factor2)), data=data), AND result.item=aov(ReadingTime~Factor1*Factor2 + Error(Item/(Factor1*Factor2)), data=data). Both both functions returned significances. Now I am confused as to what methods I should use to do post-hoc analysis and how to use them. I attempted TukeyHSD, but it gave me an error message as shown below. I am not sure how to apply TukeyHSD properly since the help file was not clear enough on that. So I wonder if you can tell me what went wrong with my function. > TukeyHSD(result.subject, c('Factor1','Factor2'), ordered=TRUE) Error...

Hi, I am seeking help with an error when running the example from R Documentation for TukeyHSD. The error occurs with any example I run, from any text book or website. thank you... > plot(TukeyHSD(fm1, "tension")). Error in plot(confint(as.glht(x)), ylim = c(0.5, n.contrasts + 0.5), ...) : error in evaluating the argument 'x' in selecting a method for function '...

Attached follows a patch to a minor change in the plot method of the TukeyHSD class (package stats). Basically it defines main= and xlab= as formal arguments of the plot function, with reasonable default values, passing them to title(), instead of using hard-coded only values for main= and xlab=. This way it's possible to change the title and xlab of the plots created us...

Hello R Developers, Dear R-developers, I ran some standard tests with currently (today morning) compiled R release candidate in Linux R 2.12.1 RC (2010-12-13 r53843). Some of these tests used plot.TukeyHSD function. This worked OK on the screen (X11 device), but PostScript file could not be rendered. The following example had the problem with me: postscript(file="tukeyplot.ps") example(plot.TukeyHSD) dev.off() I couldn't view the resulting file with evince in Linux nor in the standard...

Saludos, tengo una base de datos con tres variables: especie, tratamiento y abundancia. Quiero hacer comparaciones múltiples con Tukey. Ya había usado el comando TukeyHSD(aov(x~y)). La cuestión esque ahora tengo varios niveles del factor especie y varios niveles del factor tratamiento. Si hago TukeyHSD(aov(abundancia~especie:tratamiento)) me salen todas las comparaciones posibles, pero yo sólo estoy interesado en que para cada especie (i.e. para cada nivel del fact...