dplyr

Bee Trait Aggregate

This is a follow-up to Bee Trait Data Basics.

For data.frame df this dplyr code calculates the average day of emergence by Family:

Famdf <- group_by(df[!is.na(Pheno_mean)], Family)
summarize(Famdf, emerge_Julian = mean(Pheno_mean))

Modify the code to calculate and print the average emergence of a bee species in each Subfamily.
Use max() to determine the latest bee emergence in each Family.

[click here for output]