First cut at changes discussed in issue #3423

This version also needs cleaning (the two cholfact methods defined in the loop should be condensed to one) but I foolishly went ahead and used these facilities in the latest version of the MixedModels package.