Here we'll extract the standard error of the mean for groups form the data to explore how we can define other functions in a call to aggregate() . Note that the function could equivalently be defined as a stand-alone function and called from within aggregate() - either way works.