residualsVGAM: Extract Residuals from models of VGAM
In PAsso: Assessing the Partial Association Between Ordinal Variables

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