Question

请帮助重新排列行和cloumns 。我有两张桌子（垫子和垫子1）表1是垫

mat <- matrix("", 5, 7)
mat[c(1),] <- sample(c("AA", "AG", "GG", "--"), 7, TRUE)
mat[c(2),] <- sample(c("AA", "AT", "TT", "--"), 7, TRUE)
mat[c(3),] <- sample(c("CC", "CG", "GG", "--"), 7, TRUE)
mat[c(4,5),] <- sample(c("CC", "CG", "GG", "--"), 7, TRUE)
colnames(mat)= c("Line_1", "Line_2", "Line_3", "Line_4", "Line_5", "Line_6", "Line_7")
row.names(mat)= c("Marker_1", "Marker_2", "Marker_3", "Marker_4", "Marker_5")

表2是mat1

mat1<-matrix(c("1", "1","3","8","10", 2.3, 4.3, 45.5, 102.8, 74.8), nrow = 5)
row.names(mat1)= c("Marker_1", "Marker_2", "Marker_3", "Marker_4", "Marker_5")
colnames(mat1)= c("Linkage_Group", "Map_Position")

我用cbind组合了两个表

mat2 <- cbind(mat, mat1)

然而，我需要通过比较每一行来重新分配数据，并通过以下格式的比较列获得结果。如果标记数据匹配则为“Mono”，如果其中一个或两个都丢失则为“Mis”，如果不匹配则为“Poly＆＃39;

Comparison           Marker  Linkage_Group Map_Position   LineX1  LineX2   Mon/Poly
Line_1/Line_2       "Marker_1"      "1"     "2.3"       "AA"        "--"      Mis
Line_1/Line_3       "Marker_1"      "1"     "2.3"       "AA"        "TT"      Poly
Line_1/Line_4       "Marker_1"      "1"     "2.3"       "AA"        "TT"      Poly
Line_1/Line_5       "Marker_1"      "1"     "2.3"       "AA"        "AT"      Poly
Line_1/Line_6       "Marker_1"      "1"     "2.3"       "AA"        "--"      Mis
Line_1/Line_7       "Marker_1"      "1"     "2.3"       "AA"        "AA"      Mono

.
.
.

以上示例将有21个线组合。

这里有很多R大师......请帮助我！

Answer 1

检查这是否有帮助：

列名Line的组合。使用grep查找列名称中的模式

indx <- combn(colnames(mat2)[grep("Line", colnames(mat2))], 2)

将列名indx按indx列拆分，以便每个列表元素是两列的组合

  res <- do.call(rbind, lapply(split(indx, col(indx)), function(x) {
 m1 <- mat2[, c("Linkage_Group", "Map_Position", x)] ##subset mat2 
#compare the concerned columns rowwise  
m2 <- cbind(Comparison = paste(x, collapse = "/"), Marker = row.names(m1), m1, 
    `Mon/Poly` = apply(m1[, 3:4], 1, function(x) ifelse(any(grepl('--', x)), "Mis",
 ifelse(x[1] == x[2], "Mono", "Poly"))))
}))

如果需要，请订购

res1 <- res[order(res[,"Marker"], res[,"Linkage_Group"], res[,"Map_Position"]),]

tail(res1,3)
#       Comparison      Marker     Linkage_Group Map_Position Line_1 Line_2
# Marker_5 "Line_5/Line_6" "Marker_5" "10"          "74.8"       "--"   "--"  
# Marker_5 "Line_5/Line_7" "Marker_5" "10"          "74.8"       "--"   "--"  
# Marker_5 "Line_6/Line_7" "Marker_5" "10"          "74.8"       "--"   "--"  
#         Mon/Poly
# Marker_5 "Mis"   
# Marker_5 "Mis"   
# Marker_5 "Mis"

BTW：最好将它们存储在data.frame

中

如何重新排列R中的行和列

1 个答案: