ich würden gerne das Ergebnis des Tukey-HSD als Heatmap darstellen. In dem Beispiel unten werden die Konfidenzintervalle dargestellt, mich interessiert aber "nur" die p-Werte.
- Code: Alles auswählen
set.seed (0)
#Daten erstellen
data <- data.frame(group = rep(c("A", "B", "C"), each = 30),
values = c(runif(30, 0, 3),
runif(30, 0, 5),
runif(30, 1, 7)))
#Die ersten sechs Zeilen anzeigen
head(data)
#einfaktorielles ANOVA-Modell anpassen
model <- aov(values~group, data=data)
#Sehen Sie sich die Modellausgabe an
summary(model)
#Tukey Test durchführen
TukeyHSD(model, conf.level=.95)
#Konfidenzintervalle plotten
plot(TukeyHSD(model, conf.level=.95), las = 2)
Vielen vielen Dank!!!! Ich hab da grad einen Knoten im Hirn...
Viele Grüße
Stefano
