Je cherche à réaliser ce graphique avec le même format avec tidyverse.
Exemple :
set.seed(123)
X <- data.frame(a = rnorm(10), b = rnorm(10), c = rnorm(10,2))
m1 <- min(X[,1:2]) - 0.5
m2 <- min(X[,2:3]) - 0.5
m3 <- max(X[,c(1,3)]) + 0.5
library(reshape2)
df <- melt(X)
pval <- pairwise.wilcox.test(df$value, df$variable)$p.value
par(bty="n")
boxplot(X, las = 1, col = "wheat", bty="l", ylim = range(X) + c(-1, 1))
Merci pour vos conseils éclairés
A. Blanchet