Three-way crosstabulation using R

Here is one way to do it.