Fits candidate probability distributions to residuals, optionally with parameters that depend on state variables.
fit_residual_distribution(
residuals,
candidates = c("normal", "t", "skew-normal"),
conditional_on = NULL,
data = NULL
)List with best fitting distribution and parameters
Numeric vector of residuals
Character vector of distribution families to try
Formula for conditional parameters (optional)
Data frame (required if conditional_on specified)