compute_gaulss_scores: Compute per-observation scores for gaulss family
Description
gaulss uses `tau = 1/sigma` with logb link and defines `family$sandwich`,
which prevents the generic GLM score path. Scores are computed analytically
from the log-likelihood and transformed to linear predictor space.
Usage
compute_gaulss_scores(b, X)
Value
Score matrix (n_obs x p).
Arguments
- b
Fitted GAM object.
- X
Model matrix with `lpi` attribute.