Я использую функцию сильфона для сравнений между уровнями группы с pairwise.wilcox.test:Соединить кадры данных по горизонтали в пределах цикла в R
pairwise.wilcox <- function(data, x, factor) {
require(reshape2)
pw <- data.frame()
for(i in x){
pw <- pairwise.wilcox.test(data[ ,i], factor, p.adj = "bonf")
#pairwise.wilcox.test output is reformatted
pw <- melt(pw[[3]])
#the first two columns of the reformatted output are merged
pw$Var1 = paste(pw$Var2, pw$Var1, sep=" - ")
#the second column is dropped
drops <- c("Var2")
pw <- pw[,!(names(pw) %in% drops)]
colnames(pw) <- c("Compared groups", "p-value")
print(pw)
}
}
pairwise.wilcox(tabel2, c(2:4), tabel2$lot)
Это отрывок выхода:
Compared groups p-value
1 1StL - 3StL 1
2 1StL - IP 1
3 1StL - ISR 1
4 1StL - ISU 1
5 1StL - StM 1
6 3StL - 3StL NA
7 3StL - IP 1
...
Compared groups p-value
1 1StL - 3StL 1,0000000
2 1StL - IP 0,1092484
3 1StL - ISR 1,0000000
4 1StL - ISU 1,0000000
5 1StL - StM 1,0000000
6 3StL - 3StL NA
7 3StL - IP 1,0000000
...
Compared groups p-value
1 1StL - 3StL 1
2 1StL - IP 1
3 1StL - ISR 1
4 1StL - ISU 1
5 1StL - StM 1
6 3StL - 3StL NA
7 3StL - IP 1
...
То, что я хочу, чтобы сцепить по горизонтали эти кадры данных по общему первому столбцу, так что в результате выводится выглядит следующим образом:
Compared groups p-value p-value p-value
1 1StL - 3StL 1 1,0000000 1
2 1StL - IP 1 0,1092484 1
3 1StL - ISR 1 1,0000000 1
4 1StL - ISU 1 1,0000000 1
5 1StL - StM 1 1,0000000 1
6 3StL - 3StL NA NA NA
7 3StL - IP 1 1,0000000 1
...
Как это сделать? Имейте в виду, что у меня больше переменных.
Мой набор данных (не все данные):
structure(list(lot = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L,
6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L,
6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L,
6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L,
6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L,
6L, 6L, 6L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 5L, 5L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L), .Label = c("1StL",
"3StL", "IP", "ISR", "ISU", "StM"), class = "factor"), RS.5_1 = c(2,
3, 3, 2, 2, 2, 2, 3, 1, 1, 2, 3, 3, 2, 2, 1, 3, 2, 3, 2, 3, 1,
3, 1, 3, 2, 3, 2, 2, 2, 2, 3, 3, 2, 2, 3, 2, 3, 2, 2, 2, 3, 3,
3, 3, 2, 3, 3, 3, 2, 2, 2, 2, 3, 2, 3, 3, 2, 2, 1, 3, 2, 2, 3,
2, 2, 3, 3, 3, 3, 3, 2, 2, 2, 2, 1, 3, 3, 3, 3, 3, 3, 1, 1, 3,
2, 1, 2, 3, 2, 2, 3, 3, 3, 2, 3, 1, 3, 2, 2, 1, 1, 3, 2, 1, 3,
1, 2, 2, 3, 2, 3, 2, 2, 2, 1, 3, 1, 2, 1, 1, 2, 2, 2, 3, 2, 2,
2, 1, 3, 3, 1, 1, 2, 2, 2, 3, 1, 2, 2, 2, 1, 2, 2, 1, 3, 3, 3,
2, 2, 2, 3, 2, 1, 2, 3, 3, 1, 1, 2, 2, 2, 3, 1, 1, 1, 2, 3, 2,
2, 3, 2, 2, 3, 3, 3, 1, 3, 3, 1, 2, 3, 3, 3, 1, 3, 3, 3, 3, 1,
3, 1, 3, 1, 2, 3, 3, 2, 3, 3, 1, 3, 2, 3, 3, 3, 2, 3, 3, 3, 3,
3, 3, 3, 3, 3, 3, 2, 3, 3, 3, 2, 2, 1, 3, 2, 3, 2, 2, 2, 1, 3,
3, 3, 1, 2, 3, 1, 1, 3, 3, 3, 2, 2, 1, 1, 2, 3, 1, 1, 2, 1, 1,
2, 3, 3, 3, 1, 2, 3, 2, 2, 2, 2, 3, 3, 1, 1, 3, 3, 2, 2, 2, 1,
3, 2, 2, 2, 2, 1, 3, 2, 2, 3, 3, 3, 1, 3, 2, 3, 2, 2, 2, 3, 2,
3, 3, 2, 2, 2, 1, 2, 2, 2, 2, 2, 3, 2, 1, 2, 2, 1, 3, 2, 2, 2,
3, 3, 2, 1, 2, 1, 2, 3, 3, 2, 3, 3, 2, 1, 3, 1, 2, 2, 2, 3, 1,
2, 3, 3, 3, 3, 2, 2, 2, 2, 3, 3, 2, 1, 2, 2, 3, 1, 3, 1, 3, 3,
3, 3, 3, 3, 2, 3, 3, 3, 2, 1, 1, 1, 3, 1, 1, 1, 3, 2, 2, 2, 3,
3, 1, 3, 1, 2, 1, 1, 3, 3, 3, 2, 1, 1, 1, 2, 3, 1, 3, 3, 3, 2,
3, 2, 2, 2, 1, 2, 2, 2, 1, 1, 1, 2, 2, 3, 3, 2, 2, 2, 1, 1, 1,
3, 3, 3, 3, 3, 2, 3, 3, 1, 1, 3, 2, 2, 2, 3, 1, 3, 2, 3, 3, 3,
2, 2, 3, 1, 2, 3, 3, 2, 3, 2, 3, 3, 3, 2, 2, 1, 3, 2, 3, 3, 2,
3, 3, 3, 2, 1, 3, 3, 1, 1, 3, 1, 3, 3, 3, 3, 1, 2, 3, 2, 3, 2,
3, 1, 3, 2, 2, 2, 2, 2, 1, 2, 1, 3, 2, 2, 2, 3, 2, 2, 2, 2, 2,
1, 2, 3, 3, 3, 3, 2, 2, 3, 1, 2, 3, 2, 1, 3, 2, 3, 3, 3, 3, 3,
3, 3, 3, 3, 1, 3, 1, 3, 3, 3, 3, 3, 2, 2, 1, 3, 3, 3, 3, 2, 3,
3, 3, 2, 2, 3, 2, 1, 3, 1, 1, 2, 2, 3, 3, 3, 3, 2, 2, 3, 2, 2,
2, 2, 2, 1, 3, 1, 2, 3, 1, 3, 3, 2, 2, 3), RS.5_2 = c(2, 2, 3,
3, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 3, 2, 1, 1, 2, 2, 2, 1,
2, 1, 2, 2, 1, 1, 3, 2, 2, 2, 2, 1, 1, 1, 1, 2, 2, 1, 1, 2, 1,
1, 3, 3, 1, 1, 2, 2, 1, 1, 3, 2, 1, 2, 1, 1, 2, 2, 1, 1, 1, 2,
2, 1, 3, 3, 2, 1, 1, 1, 2, 2, 2, 2, 3, 2, 3, 1, 1, 1, 1, 2, 1,
2, 2, 2, 2, 1, 1, 1, 1, 3, 1, 1, 2, 2, 1, 2, 1, 1, 1, 2, 2, 1,
1, 3, 2, 1, 1, 1, 1, 1, 2, 1, 3, 2, 2, 1, 1, 1, 2, 2, 1, 1, 2,
3, 1, 2, 2, 1, 1, 2, 3, 1, 2, 2, 2, 1, 2, 1, 1, 2, 3, 1, 1, 1,
1, 3, 2, 2, 3, 3, 3, 2, 1, 1, 1, 1, 2, 2, 2, 2, 3, 2, 2, 2, 1,
1, 1, 3, 1, 2, 2, 1, 3, 1, 2, 1, 2, 3, 1, 1, 3, 1, 2, 3, 1, 2,
3, 1, 3, 3, 2, 2, 1, 3, 1, 1, 1, 2, 2, 3, 1, 1, 2, 1, 1, 1, 1,
1, 1, 2, 2, 2, 2, 1, 1, 1, 3, 1, 1, 2, 3, 1, 3, 1, 1, 1, 3, 3,
1, 1, 1, 1, 1, 2, 2, 3, 2, 1, 1, 3, 2, 3, 1, 2, 2, 2, 2, 2, 3,
1, 1, 2, 3, 3, 2, 1, 1, 1, 3, 3, 2, 2, 1, 1, 1, 1, 1, 1, 1, 2,
2, 2, 2, 2, 1, 1, 1, 2, 2, 1, 2, 1, 1, 2, 2, 1, 2, 2, 1, 2, 1,
1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 2, 2, 2, 3, 2, 1, 1, 1, 2,
2, 2, 1, 2, 1, 1, 1, 2, 2, 3, 2, 2, 3, 2, 3, 2, 2, 2, 1, 2, 1,
1, 3, 1, 2, 1, 2, 3, 1, 1, 1, 2, 2, 3, 3, 1, 3, 1, 2, 3, 2, 3,
2, 1, 1, 1, 3, 1, 3, 2, 2, 2, 3, 1, 1, 1, 1, 1, 3, 3, 3, 3, 1,
2, 1, 1, 2, 2, 1, 2, 2, 1, 2, 1, 1, 1, 3, 1, 3, 3, 3, 1, 3, 2,
1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 3, 3, 1, 1, 2, 1, 1, 2, 3, 2,
3, 3, 3, 2, 2, 1, 1, 2, 1, 3, 2, 3, 3, 2, 3, 1, 3, 3, 3, 2, 2,
3, 1, 2, 3, 3, 2, 3, 1, 3, 3, 3, 1, 2, 1, 3, 2, 3, 3, 2, 3, 1,
1, 2, 1, 2, 1, 1, 1, 2, 2, 2, 3, 1, 2, 1, 3, 3, 2, 2, 1, 2, 1,
3, 3, 2, 2, 3, 2, 2, 2, 1, 3, 3, 3, 1, 1, 1, 1, 2, 1, 3, 1, 2,
3, 2, 1, 1, 2, 1, 2, 2, 1, 1, 1, 2, 1, 3, 2, 2, 3, 3, 3, 1, 1,
3, 3, 1, 3, 1, 2, 3, 2, 1, 1, 1, 1, 2, 3, 3, 3, 2, 1, 1, 2, 1,
2, 1, 1, 3, 1, 2, 2, 2, 1, 1, 3, 1, 3, 2, 1, 1, 1, 1, 1, 1, 3,
1, 2, 1, 1, 1, 3, 1, 1, 1, 1, 2, 1), RS.5_3 = c(1, 3, 3, 3, 1,
1, 3, 2, 3, 2, 1, 3, 1, 3, 3, 3, 3, 1, 1, 1, 2, 1, 2, 1, 3, 1,
2, 1, 3, 2, 1, 1, 3, 1, 2, 2, 1, 2, 3, 3, 1, 1, 2, 3, 3, 2, 3,
3, 1, 2, 1, 2, 3, 1, 3, 1, 2, 3, 3, 2, 1, 3, 1, 2, 1, 2, 3, 2,
2, 3, 2, 1, 1, 1, 3, 2, 1, 3, 3, 1, 3, 1, 2, 2, 3, 3, 2, 2, 3,
2, 2, 2, 2, 1, 2, 3, 3, 2, 2, 1, 2, 3, 1, 1, 3, 1, 1, 1, 1, 3,
2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 1, 1, 2, 1, 1, 2, 2, 2, 3, 1, 2,
2, 2, 1, 3, 1, 2, 1, 2, 2, 1, 3, 3, 2, 3, 3, 3, 1, 1, 3, 3, 1,
1, 2, 2, 3, 1, 1, 2, 1, 2, 3, 2, 2, 2, 2, 3, 2, 1, 2, 2, 2, 2,
3, 1, 3, 3, 3, 3, 2, 1, 1, 2, 3, 1, 3, 3, 3, 2, 2, 2, 2, 3, 2,
3, 1, 3, 2, 2, 3, 1, 2, 1, 1, 3, 1, 3, 1, 2, 2, 2, 1, 1, 1, 1,
3, 3, 1, 1, 3, 2, 3, 2, 1, 3, 3, 3, 2, 1, 3, 3, 3, 3, 3, 3, 1,
1, 1, 1, 3, 3, 3, 2, 2, 1, 3, 3, 3, 3, 2, 3, 3, 2, 3, 2, 2, 1,
3, 3, 3, 3, 1, 1, 1, 3, 3, 1, 2, 1, 1, 3, 1, 3, 2, 1, 3, 2, 2,
3, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 3, 2, 2, 2, 1, 1, 1, 3, 1, 3,
1, 2, 2, 1, 3, 2, 1, 2, 1, 1, 2, 3, 3, 3, 2, 1, 2, 1, 3, 2, 3,
1, 3, 3, 1, 3, 1, 1, 3, 3, 2, 2, 2, 1, 2, 1, 3, 2, 2, 1, 1, 3,
2, 1, 1, 2, 3, 2, 2, 3, 1, 2, 3, 3, 1, 3, 2, 2, 1, 3, 3, 3, 3,
1, 2, 2, 2, 3, 1, 3, 2, 3, 2, 2, 3, 1, 3, 2, 2, 3, 3, 1, 1, 1,
1, 2, 1, 2, 1, 3, 1, 2, 1, 1, 1, 3, 2, 3, 3, 3, 1, 3, 3, 1, 1,
1, 1, 2, 2, 1, 1, 1, 1, 1, 3, 3, 1, 1, 2, 1, 1, 3, 3, 3, 3, 3,
3, 2, 3, 1, 2, 2, 1, 3, 1, 3, 2, 2, 3, 1, 3, 3, 3, 2, 1, 3, 1,
3, 3, 3, 2, 2, 1, 3, 3, 3, 2, 2, 1, 3, 3, 3, 3, 1, 3, 1, 1, 3,
1, 1, 1, 2, 1, 3, 2, 2, 3, 1, 1, 2, 3, 3, 3, 1, 2, 3, 1, 3, 2,
3, 3, 3, 2, 3, 3, 2, 3, 2, 3, 2, 3, 1, 2, 3, 3, 3, 1, 2, 3, 2,
1, 3, 2, 1, 1, 3, 1, 1, 1, 1, 3, 3, 1, 3, 3, 3, 2, 3, 1, 3, 3,
3, 1, 3, 2, 3, 2, 1, 1, 2, 1, 3, 3, 3, 3, 3, 3, 3, 1, 2, 3, 1,
1, 3, 1, 1, 3, 3, 1, 1, 1, 1, 3, 2, 1, 1, 1, 1, 2, 2, 1, 2, 2,
1, 1, 3, 1, 1, 3, 3, 2, 1, 1)), .Names = c("lot", "RS.5_1", "RS.5_2",
"RS.5_3"), row.names = c(NA, -582L), class = "data.frame")
благодарственного вы, @plannapus! – Iurie