stats::sigma() is not defined for multivariate gllvmTMB() fits: a
single residual scale is not meaningful across mixed response families.
Use family-specific dispersion extractors instead (for example
extract_Sigma() on Gaussian traits).
Usage
# S3 method for class 'gllvmTMB_multi'
sigma(object, ...)