Fitting Multilevel Factor Models