If you have 2 cross classifying variables you can use `rowSums` and `colSums` to produce margin totals on an `xtabs` output. But how can it be done if you have 3 classifying variables (ie margin totals in each sub table)?