Lander chp 11 aggregate function using formula and cbind