Code: Alles auswählen
MoVe <- matrix(NA, 12, 12)
for (i in 2:12)
for (j in 1:(i-1)) {
MoVe[i,j] <- with(d, wilcox.test(LUFTTEMPERATUR[Monat==i], LUFTTEMPERATUR[Monat==j])$p.value)
MoVe[j,i] <- MoVe[i,j]
}
Code: Alles auswählen
MoVe <- matrix(NA, 12, 12)
ind <- which(lower.tri(MoVe), arr.ind = TRUE)
p <- mapply(FUN=function(i,j) with(d, wilcox.test(LUFTTEMPERATUR[Monat==i], LUFTTEMPERATUR[Monat==j])$p.value),
ind[,"row"], ind[,"col"])
ergebnis <- data.frame(R=ind[,"row"], C=ind[,"col"], p=p)
ergebnis[ergebnis$p<0.05, ]
Code: Alles auswählen
ind <- expand.grid(i=1:12, j=1:12)
ind <- ind[ind$i>ind$j,]
ergebnis <- ind
ergebnis$p <- mapply(FUN=function(i,j) with(d, wilcox.test(LUFTTEMPERATUR[Monat==i], LUFTTEMPERATUR[Monat==j])$p.value),
ind$i, ind$j)
Code: Alles auswählen
wilcox.ij <- function(ij) with(d, wilcox.test(LUFTTEMPERATUR[Monat==ij[1]], LUFTTEMPERATUR[Monat==ij[2]])$p.value)
ergebnis <- as.data.frame(t(combn(1:12, 2)))
ergebnis$p <- combn(1:12, 2, FUN=wilcox.ij)