similar to: TukeyHSD test

hello, I have programmed this function to calculate the Neuman-Keuls test but I have a problem the function return an empty list and I don't know why. summary(fm1) E <- sqrt((summary(fm1)[[1]]["Residuals","Mean Sq"])/length(LR)) lst <- list() lst1 <- list() lst2 <- list() NK <- function (x) { if (length(x) == 2) { Tstudent <- t.test(subset(exple,

I can 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

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

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

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 comparisons using TukeyHSD and the glht function. The difference appears

Hi, I'm wondering if i can use the function "TukeyHSD" to perform the all pairwise comparisons of a "aov()" model with one factor (e.g., GROUP) and one continuous covariate (e.g., AGE). I did for example: library(multcomp) data('litter', package = 'multcomp') litter.aov <- aov(weight ~ gesttime + dose, data = litter) TukeyHSD(litter.aov, which =

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 Rcmdr before calling TukeyHSD or you can set na.action to

I've tried to run a Tukey post-hoc but keep getting this weird error, whether the aov was significant or not. treat_code is a dummy variable, but that 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.

Hello, When I try to use TukeyHSD in the following way it shows the confidence interval corresponding to the last factor only. > throughput.aov <- aov(Throughput~No_databases+Partitioning+No_middlewares+Queue_size,data=throughput) plot(TukeyHSD(throughput.aov)) # I expected here to see the confidence intervals for all factors but see only the last. OTOH this one works but then it is

All - I think I'm being dense, but for the life of me, I can't figure out why I get error message with the code below. I have data that looks like param level perc.surv asin.tran DO 3 0.6864407 0.9764544 DO 3 0.1250000 0.3613671 DO 3 0.8738739 1.2077299 DO 4 0.4615385 0.7468986 DO 4 0.5691057 0.8547258 DO 4 0.8504673 1.1737517 DO 6 0.5494505 0.8349297 DO

Hey, I have this TukeyHSD output from which I would like to extract only the P-values (p adj, last number). The problem is that the test output is a character list. How can I "break" this sentence to separate the Pv? Tukey multiple comparisons of means 95% family-wise confidence level Fit: aov(formula = Fe1$Fe ~ Fe1$genotype) $`Fe1$genotype` diff lwr upr

Hi all, I search the archive for finding a simple solution for using TukeyHSD with a multistratum aov result (a repeated emasure anova). The Question have been asked but I've found no clear answer. res<-aov(y~Fa*Fb+Error(Subject/(Fa*Fb)) ) I think that the problem is that res is an aovlist object instead of the "aov" object required by TukeyHSD. Is there an easy solution to

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 <-

Hello everybody! I?m working with a dataset from eleven field trails on barley fertilization. I use R 1.6.2 (Windows) It is quite easy to fit aov() objects to the dataset. The call: > (l1t4y.aov <- aov(Yield ~ Trial + Treatment, data=led1t4b)) Results in an object with this anova table: Df Sum Sq Mean Sq F value Pr(>F) Trial 10 121423585 12142358 63.499

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. The labels are integers.

Hi, the function TukeyHSD gives incorrect results for balanced incomplete block designs, as the example below shows, but I can only half fix it. There are two problems, 1. It uses model.tables to estimate treatment means, 2. It uses the wrong standard error The first problem can be fixed using dummy.coef, if the lines > TukeyHSD.aov function (x, which = seq(along = tabs), ordered = FALSE,

R Version 2.0.0 Under development (unstable) (2004-05-07) I have been wanting to use TukeyHSD for two and three way aov's, as they are especially simple for students to use correctly. I realize model simplification is usually a preferred methodology, but my disciplinary enertia and simplicity of TukeyHSD is also compelling. However, the lack of names on the comparisons for the higher order

Dear R people, I have a couple of questions about post-doc analyses for 2 by 2 within subjects ANOVA. I conducted a psycholinguistic study that combined a 2 by 2 design and a latin square design. Specifically, I had 32 items each of which generated 4 conditions. Participants saw each of the 32 items only once: 8 in Condition A, 8 in B, 8 in C, and 8 in D. The table below serves as an example.

hello, I need to groupe some means and I wonder how to do to get critical range table for tukey or for Neuman-Keuls I think it's possible to do that with the qtukey() function? thanks. _____________________________________________________________________________ [[alternative HTML version deleted]]

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