Variance component estimation in linear models