random factors are now directly defined from the model by corr=corCompSymm(form = model$call$random).
blocking is added
create a function which compais different random factor structures