Based on the combination of common factors for the best-fitting model, maximum likelihood genetic factor scores were estimated by computing the conditional likelihood of the twin pairs' item responses, weighted by the joint likelihood of the factor score estimates. This factor scoring model was iteratively fitted, separately for each twin zygosity group, to each twin pair's raw data to obtain estimates of scores on the factors of the best-fitting model for each individual.
