As suggested by @ahunter and @ncortesi , Apply()'s arguments could be renamed to match those of base apply().
data -> X fun -> FUN margins -> MARGINS
To be reconsidered for version v3.0.0.