more r code on log-normal mixed linear models