Using Custom Aggregate Functions with cast() in R reshape2: A Practical Guide to Resolving the Limitation of vapply and fill=0