Variance analysis and components

# Variance components
library(nlme)
varcomp (lme(dep_trait ~ 1, random = ~ 1 | grouping_factor, data=data_set))

# Variance analysis
summary(aov(dep_trait ~ grouping_factor, data=set_of_data))