more sim code for mixed linear models