具有因子名称而不是计数的列联表

Question

我有一个包含3列的数据框。 第一个包含我的类别的名称，第二个（q2）和第三个（q3）列包含其他因素。 ES：

q1 q2 q3
A  w  y
G  x  x
F  y  z

我正在使用

t1 <- xtabs(~q2+q3, data=db)

创建一个列联表并知道有多少个组合类别。

我如何在itersection单元格中获取q1列中的因子名称？

     x y z
   w 0 1 0
   x 1 0 0
   y 0 0 1

在示例中，我想要A，G或F.而不是1。

Answer 1

一个选项是使用q2和q3作为索引将q1分配给t1 。

> dn <- lapply(db[2:3], levels)
> t1 <- array(NA, lengths(dn), dn)
> t1[cbind(db$q2, db$q3)] <- levels(db$q1)[db$q1]
> t1
   q3
q2  x   y   z  
  w NA  "A" NA 
  x "G" NA  NA 
  y NA  NA  "F"