Added script to compute variance components; unfinished, just up to model fitting.