En R, puedo devolver los resultados del recuento utilizando los nombres de columna específicos que me interesan como una matriz como se muestra a continuación.

require("plyr")
bevs <- data.frame(cbind(name = c("Bill", "Llib"), drink = c("coffee", "tea", "cocoa", "water"), cost = seq(1:8)))
count(bevs, c("name", "drink"))

# produces
  name  drink freq
1 Bill  cocoa    2
2 Bill coffee    2
3 Llib    tea    2
4 Llib  water    2

¿Cómo puedo obtener el resultado del conteo de dos nombres de columna específicos en una matriz que tiene columnas: todas las bebidas únicas, filas: todos los nombres únicos y celdas: freqs (como a continuación)?

     cocoa  coffee tea water
Bill   2      2     0   0
Llib   0      0     2   2

P.S: Obviamente, la solución no necesita usar plyr.