similar to: tukey or Neuman-Keuls

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 am embarrassed to have to ask this but can anyone tell me of a Tukey multiple comparisons procedure available for R? I have looked through the search page, through the FAQ, and in the index of V&R (1999), and I still can't find such a thing. I see there is a ptukey function and a qtukey function but that is as far as I got. Do I need to roll my own? -- Douglas Bates

Hello, I think that to apply the Neuman-Keuls test under R it's possible with the TukeyHSD function and I know that to do a Neuman-Keuls test I have to sort the means but I can't use the TukeyHSD function I don't understand how to do to apply this function to my data. In the R help they use the TukeyHSD function to a aov() do you think that I have to do teh same? thanks.

bonjour, je vous écris pour savoir si vous aviez trouver réponse pour l'écriture du test de Newman & Keuls sous R car j'en ai aussi besoin merci encore cordialement. --------------------------------- [[alternative HTML version deleted]]

I am using R 1.1 with Redhat 6.2 and RW 1.001 with Win98 (the upkey doesn't work on my IBM either as has been previously reported by others). The function aov doesn't return either the residuals or the residual degrees of freedom for split-plot designs. If you use the following code from Baron and Li's "Notes on the use of R for psycology experiments and questionnaires"

Hi, I'm building Tukey's table using qtukey function. It happens that I can't get the values of Tukey's one degree of freedom and also wanted to eliminate the first column. The program is: Trat <- c(1:30) # number of treatments gl <- c(1:30, 40, 60, 120) # degree freedom tukval <- matrix(0, nr=length(gl), nc=length(Trat)) for(i in

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

Hi, I'm trying to make sense of the options for multiple comparisons options in R. I've found the following options: pairwise.t.test, which provides standard t-tests, with options for choosing an appropriate correction for multiple comparisons TukeyHSD, which provides the usual Tukey test glht(package multcomp), which provides a variety of options >From the help list, it appears

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

Hello, I am interested in parametric multi comparison tests such as Dunnett, Duncan, Tukey, Newman-Keuls, Bonferonni, Scheffe, and non-parametric tests such as Kruskal-Wallis, and Mann-Whitney U. Are there packages that include most of these tests in each category? Many packages exist for an individual test but their outputs vary in great detail (test statistics, p-values, etc.)

Jose - Before implementing SNK and Duncan's, you may want to be aware of some criticisms of these methods: >From Hsu (1996), "Newman-Keuls multiple range test is not a confident inequalities method and cannot be recommended." "Duncan's multiple range test is not a confident inequalities method and cannot be recommended either. In the words of Tukey (1991),

> Tukey divides the points into three groups, not the x and y values separately. > I'll try to get hold of the book for a direct quote, might take a couple of days. Ah well, I can't get it for a week. But the fact that it's often called Tukey's three group line (try a search on *tukey three group line* and you'll get plenty of hits) is pretty much a giveaway. On Mon,

The problem or actual R implementation relies on an assumption that median(x[i] | x[i] <= quantile(x, 1/3)) == quantile(x, 1/6) which reveals not to be true despite very trustful appearance. If we continue with the example of x=y=1:9 then quantile(x, 1/6)=2.5 (here quantile() is taken in C-code sens, not R's one) while median(y[i] | x[i] <= quantile(x, 1/3))=2 On the other sample's

And with "equally spaced" I obviously meant "of equal size". It's getting too hot in the office here... On Wed, May 31, 2017 at 4:39 PM, Joris Meys <jorismeys at gmail.com> wrote: > Seriously, if a method gives a wrong result, it's wrong. line() does NOT > implement the algorithm of Tukey, even not after the patch. We're not > discussing Excel here,

Sorry, I have seen it too late that we had different tab width in the original file and my editor. Here is the patch with all white spaces instead of mixing tabs and white spaces. Serguei. Le 29/05/2017 ? 15:13, Serguei Sokol a ?crit : > Here is an attached patch. > > Best, > Serguei. > > Le 29/05/2017 ? 12:21, Serguei Sokol a ?crit : >> The problem or actual R

Le 31/05/2017 ? 16:39, Joris Meys a ?crit : > Seriously, if a method gives a wrong result, it's wrong. I did not understand why you and others were using term "wrong" based on something that I was considering as just "different" implementation. More thorough reading revealed that I have overlooked this phrase in the line's doc: "left and right /thirds/ of the

Hi everyone, I woudl like to apply a Tukey post hoc after a repeated measure ANOVA. I followed the suggestions that I found in this help -list especially this one: /[R] Tukey HSD (or other post hoc tests) following repeated measures ANOVA You want to use lme() in package nlme, then glht() in the multcomp package. This will give you multiplicity adjusted p-values and confidence intervals.

> On 31 May 2017, at 16:40 , Joris Meys <jorismeys at gmail.com> wrote: > > And with "equally spaced" I obviously meant "of equal size". It's getting > too hot in the office here... We have a fair amount of cool westerly wind up here that I could transfer to you via WWTP (Wind and Weather Transport Protocol). If you open up a sufficiently large pipe,

hello, I wanna know how to create a list of list if it's possible and if it isn't possible how to do without. thanks. _____________________________________________________________________________ [[alternative HTML version deleted]]

Dear R-list, Could anybody tell me (or give me a tip) of how to implement the Duncan distribution in R? I've been trying to make a new and more flexible function for multiple comparison of means: Tukey, SNK and Duncan, from 'aov' objects, like TukeyHSD function. For while, it is running nice (Tukey and SNK), for simple design (completely randomized, randomized block and Latin