By default, this function will provide means at -1SD, the mean, and +1SD for
continuous variables, and at each levele of categorical variables. This can
be overridden with the levels
parameter.
If there are additional covariates in the model other than what are selected
in the function call, these variables will be set to their respective means.
In the case of a categorical covariate, the results will be averaged across
all its levels.