lmer(Z_TP~Z_Length + (1|Lake) + (1|Species), data = data, REML=TRUE)