ich stehe mal wieder völlig auf dem Schlauch. Ich möchte eine ANOVA durchführen und danach als Post Hoc einen Tukey HSD. Am Ende möchte ich eine Übersicht haben, welche Spalten sich von welchen signifikant unterscheiden.
Gegeben sei hier eine Testtabelle:
Code: Alles auswählen
Testtabelle <- matrix(nrow=100, ncol=10)
colnames(Testtabelle) <- c(
paste("Proband", 1:10, sep=""))
rownames(Testtabelle) <- paste("gene", 1:100, sep="")
for (i in 1:100) {
Werte<- rnorm(n = 10, mean = 10, sd = 3)
Testtabelle[i,] <- c(Werte)
}
Wie mache ich das jetzt konkret? Ich habe mit dem tidyr package gespielt (Muss ich die Daten umorgenisieren? Muss ich die Spaltennammen als eine Spalte und die Beobachtungen in die zweite Spalte bringen wobei Zeilennamen dann als Information wegfallen?) und mit den Befehlen "aov" und "TukeyHSD" (Was hat das mit formula und ~ bei 10 Spalten auf sich?), Komme leider zu keinem brauchbaren Ergebnis.
Über ein brauchbares Script anhand der Testtabell würde ich mich sehr freuen, dann könnte ich versuchen es nachzuvollziehen.
Vielen Dank
MfG
Werner