Closed aharach closed 9 years ago
deviance_residual <- function(act,fit){
#Residual for Gamma Log GLM model
r <- (act-fit)/fit
y <- act
mu <- fit
wt <- rep(1,length(act))
d.res <- sqrt(pmax(-2 * wt * (log(ifelse(y == 0, 1, y/mu)) - (y - mu)/mu),0))
ifelse(y > mu, d.res, -d.res)
}