R/qresid_binom.R

In glmxdiag: A Collection of Graphic Tools for GLM Diagnostics and some Extensions

qresid_binom = function(model) {
  p.obs <- model$y
  n <- model$prior.weights
  y <- n * p.obs
  pi.hat <- fitted(model)
  u <- runif(n)
  binom.cdf <- pbinom(round(y)-1, size = n, prob = pi.hat) + u * dbinom(round(y), size = n, prob = pi.hat)
  res <- qnorm(binom.cdf)
  res
  
}