Copy of the dMcast function from the Matrix.utils package, since this is off CRAN and does not seem to be maintained anymore internally

Usage

dMcast(
  data,
  formula,
  fun.aggregate = "sum",
  value.var = NULL,
  as.factors = FALSE,
  factor.nas = TRUE,
  drop.unused.levels = TRUE
)