The `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, ...)

## Arguments

x |
A fitted model object returned by one of the
rstap modeling functions. See `stapreg-objects` . |

digits |
Number of digits to use for formatting numbers. |

include_X |
logical for whether or not to include estimated latent exposure covariate |

... |
Ignored. |

## Value

Returns `x`

, invisibly.

## Details

### Point estimates

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`

.

### Uncertainty estimates (MAD_SD)

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`

.

### Additional output

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).

## See also