print.stapreg.RdThe print method for stapreg objects displays a compact summary of the
fitted model. See the Details section below for descriptions of the
different components of the printed output. For additional summary statistics
and diagnostics use the summary method.
# S3 method for stapreg print(x, digits = 1, include_X = FALSE, ...)
| x | A fitted model object returned by one of the
rstap modeling functions. See |
|---|---|
| digits | Number of digits to use for formatting numbers. |
| include_X | logical for whether or not to include estimated latent exposure covariate |
| ... | Ignored. |
Returns x, invisibly.
Point estimates are medians computed from simulations.
For models fit using MCMC ("sampling") the posterior
sample is used. The point estimates reported are the same as the values
returned by coef.
The standard deviations reported (labeled MAD_SD in the print output)
are computed from the same set of draws described above and are proportional
to the median absolute deviation (mad) from the median.
Compared to the raw posterior standard deviation, the MAD_SD will be
more robust for long-tailed distributions. These are the same as the values
returned by se.
The median and MAD_SD are also reported for mean_PPD, the sample
average posterior predictive distribution of the outcome. This is useful as a
quick diagnostic. A useful heuristic is to check if mean_PPD is
plausible when compared to mean(y). If it is plausible then this does
not mean that the model is good in general (only that it can reproduce
the sample mean), however if mean_PPD is implausible then it is a sign
that something is wrong (severe model misspecification, problems with the
data, computational issues, etc.).
For GLMs with group-specific terms (see stap_glmer) the printed
output also shows point estimates of the standard deviations of the group
effects (and correlations if there are both intercept and slopes that vary by
group).