Properly scale standard deviations depending on the error distribution
scaleFam(y, model, newdata, standardize = "scale", standardize.type = "latent.linear")