A internal function to simulate surrogate residuals for models fitted by
vglm and vglm. Now, this one support "vglm" and
"vgam", and adjacent categories regression model by "vglm". This one may
need to update to support more models from VGAM.
residualsVGAM(
object,
type = c("surrogate", "sign", "general", "deviance"),
jitter = c("latent", "uniform"),
jitter.uniform.scale = c("probability", "response"),
nsim = 1L,
...
)