Variance components model