FormulaSums function

Sums (aggregates) and/or sparse model matrix with possible cross table