3

我有这种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
4

2 回答 2

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 <-b$Var1 b$Var2.merge <-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 对象来自何处，您可以更轻松地处理原始数据。

于 2014-12-16T20:06:45.530 回答
1

如果采用 data.frame 路线，这可以避免在处理因素时可能出现的问题：

> 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
于 2014-12-16T20:19:01.670 回答