如何合并意外矩阵的两个输入值?

时间:2014-12-16 19:42:43

标签: r merge summary contingency

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

我想合并两个Var1Var2的值,例如" A"和" B",在一个新的值" D"保持应变矩阵属性。

所以结果将是例子:

     Var2
Var1  D C
   D 12 9 
   C 15 9 

2 个答案:

答案 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