Calculate sigma for nested lme

First stab at calculating random variance for nested mixed models in lme