我有这种xtabs
对象:
structure(c(1, 4, 7, 2, 5, 8, 3, 6, 9), .Dim = c(3L, 3L),
.Dimnames = structure(list(Var1 = c("A", "B", "C"), Var2 = c("A", "B", "C")),
.Names = c("Var1", "Var2")), class = c("xtabs", "table"))
给出:
Var2
Var1 A B C
A 1 2 3
B 4 5 6
C 7 8 9
我想合并两个Var1
和Var2
的值,例如" A"和" B",在一个新的值" D"保持应变矩阵属性。
所以结果将是例子:
Var2
Var1 D C
D 12 9
C 15 9
答案 0 :(得分:3)
将xtabs对象读入a
,然后可以:
b<-as.data.frame.table(a)
levels(b$Var1) <- list(D=c("A","B"),C="C")
levels(b$Var2) <- list(D=c("A","B"),C="C") #Thanks @Henrik!
xtabs(Freq ~ Var1 + Var2, b)
另一种方法,您可以使用一些合并来适应更大的因素 b $ Var1.merge&lt; -b $ Var1 b $ Var2.merge&lt; -b $ Var2
# Add the new level for the factors
levels(b$Var1.merge) <- c(levels(b$Var1.merge), "D")
levels(b$Var2.merge) <- c(levels(b$Var2.merge), "D")
#Collapse the factors
b$Var1.merge[b$Var1%in%c("A","B")] <- "D"
b$Var2.merge[b$Var2%in%c("A","B")] <- "D"
xtabs(Freq ~ Var1.merge + Var2.merge, b, drop.unused.levels=T)
当然,根据原始xtabs对象的来源,您可以更容易地对原始数据进行处理。
答案 1 :(得分:1)
如果采用data.frame route:
,则可以解决处理因素时可能出现的问题> tbl2 <- rbind(colSums(tbl[1:2, ]), tbl[3, ])
> tbl2
A B C
[1,] 5 7 9
[2,] 7 8 9
> tbl3 <- cbind(rowSums(tbl2[, 1:2]), tbl2[ ,2] )
> tbl3
[,1] [,2]
[1,] 12 7
[2,] 15 8
这显示了如何使用指定的一组字母值进行折叠:
> tbl2 <- rbind(colSums(tbl[row.names(tbl) %in% c("A","B"), ]),
tbl[!row.names(tbl) %in% c("A","B"), ])
> tbl2
A B C
[1,] 5 7 9
[2,] 7 8 9
> tbl3 <- cbind(rowSums(tbl2[, colnames(tbl) %in% c("A","B")]),
tbl2[,!colnames(tbl) %in% c("A","B")] )
> tbl3
[,1] [,2]
[1,] 12 9
[2,] 15 9