Stack Overflow
Stack Overflow

Questions tagged [tukey]

The Tukey test is a single-step statistical test for multiple comparison, named after John Tukey. It is used to find significantly different mean, by comparing all pairs on means based on a studentized range distribution.

The Tukey test is a single-step statistical test for multiple comparison, named after John Tukey. It is used to find significantly different mean, by comparing all pairs on means based on a studentized range distribution.

Assumptions

  1. Independence of observations within and among groups
  2. Normal distribution within group
  3. Equal within-group variance across the groups (Homogeneity of variance)

Also known as:

  • Tukey's Honestly Significant (Difference) test
  • Tukey HSD test
  • Tukey-kramer test/method
165 questions
2
votes
2 answers

How to show Tukey groups using tidyverse and rstatix?

After a Tukey test using rstatix, I would like to add another column to the data frame with the Tukey groups. library(tidyverse) library(rstatix) ggplot(data = iris, aes(x = Species, y = Petal.Length)) + geom_boxplot() Tukey <-…
Daniel Valencia C.
  • 2,159
  • 2
  • 19
  • 38
2
votes
2 answers

stat_pvalue_manually and stat_compare_means with Tukey HSD

I am using ggpubr and ggplot2. First I did a anova test on my 3 clusters and now I would like to include my p-adj values of the TukeyHSD test in my box plot by using stat_compare_means or stat_pvalue_manually since I found out that Tukey test…
Valérie
  • 21
  • 2
2
votes
2 answers

R output tukey results only where factors or interactions are significant (conditional statements)

I'm relatively new to the R language, and using it to analyse my data. I am using a package called 'agricolae' and its HSD.test component to output tukeyHSD results based on the anova proceeding it. I have simplified the code below, but essentially…
luke.aaron
  • 21
  • 1
2
votes
2 answers

Getting tidy output from post hoc tests

Consider this data frame dat1: dat1 <- data.frame(Region = rep(c("r1","r2"), each = 100), State = rep(c("NY","MA","FL","GA"), each = 10), Loc = rep(c("a","b","c","d","e","f","g","h"),each = 5), …
Ryan
  • 1,048
  • 7
  • 14
2
votes
1 answer

How can I extract the compact letter display from a Tukey's multiple comparison using the nparcomp package?

The multcomp package has a nice compact letter display function built in, but I'm using the non-parametric multiple comparison package "nparcomp" which does not appear to have a similar feature. I've noticed there are a couple packages such as…
Dieter Kahl
  • 89
  • 6
2
votes
1 answer

How to show the results of a Tukey test with boxplots showing CLD letters

I have collected data on 216 individuals. I measured the concentration of the same 7 Substances in each individual, represented by Sub1:Sub7. The concentration of these Substances may be different in individuals from different Locations. I am…
Ryan
  • 1,048
  • 7
  • 14
2
votes
1 answer

Tukey test results on geom_boxplot with facet_grid

I'd love to place letters representing the results of a Tukey's HSD on boxplots I make in ggplot. I am aware of other posts (Tukeys post-hoc on ggplot boxplot, Tukey's results on boxplot in R) but in my case, I work with facet_grid and I am not…
Cecile
  • 527
  • 5
  • 22
2
votes
1 answer

How to perform Tukey HSD.test() on list of dataframes?

I would like to perform a Tukey Post Hoc test on a list of dataframes. As outcome I would like to have letters indicating which groups are significantly different from each other. The HSD.test() of the agricolae package does this, but I can't figure…
Marjanne
  • 21
  • 5
2
votes
1 answer

Getting wrong p-values for Tukey test for one-way mixed effect ANOVA

I am trying to replicate a textbook example of a Tukey Test for One-Way Mixed Effect ANOVA (from Statistics, William L. Hays p 581-583) but the p-values I am getting using lme & glht don't make sense The study has repeated measures of four…
Sarah Stanley
  • 21
  • 2
2
votes
0 answers

TukeyHSD results on boxplot after two-way anova

I've had similar code working fine on a one-way ANOVA and all of my commands appear to be working, but the letters aren't plotting. I'm not getting any errors in my console, so any help would be appreciated! My data: >…
Alex
  • 261
  • 2
  • 5
  • 11
2
votes
1 answer

Tukey's results on boxplot in R

I am working on placing letters representing the results of a Tukey's HSD on a boxplot that I've made. I've found a description of how to do it here: r-graph gallery, but I am getting an error after I run the function. Here is my…
Alex
  • 261
  • 2
  • 5
  • 11
2
votes
1 answer

What is the most appropriate method for contrasts in lme (nlme)?

I have a question regarding conducting post-hoc contrasts in a mixed glmm model. Essentially, my data is a series of response variables with fixed factors of 'Subject Type' (three types of subjects), 'Treatment' (four treatments), and 'Year' (two…
choppedpete
  • 29
  • 1
  • 3
2
votes
0 answers

Tukey HSD for three-way interaction term in R

For each dependent variable (time, mass, rate, quotient), I constructed a three-way ANOVA with gender, age, temperature and all interactions as the independent fixed factors. My outputs range from getting a significant two-way interaction term …
user4730784
  • 41
  • 3
2
votes
1 answer

How do I interpret the TukeyHSD output in R? (in relation to the underlying regression model)

I built a simple linear regression model with 'Score' as the dependent variable, and 'Activity' as the independent one. 'Activity' has 5 levels: 'listen' (reference level), 'read1', 'read2', 'watch1', 'watch2'. Call: lm(formula = Score ~…
fannilegoza
  • 23
  • 1
  • 1
  • 4
2
votes
1 answer

Testing mean of multiple series with ANOVA and Tukey HSD test in R

In R, I have four numerical vectors. I want to test if mean of these four vectors are the same or not, and if not, which mean is greater and which mean is smaller (basically I want to rank four vector by mean of them). This post is quite suitable…
Yang Mei Lian
  • 75
  • 1
  • 11
1
2
3
…
10 11