#------------------------------------------------------------------------------------------#
# Comandos:
#          fit.model = ajuste
#          attach(dados)
#          source("diag_cook_norm_dglm_disp.txt")
#------------------------------------------------------------------------------------------#
X = model.matrix(fit.model)
n =  nrow(X)
p =  ncol(X)
Z=X
library(dglm)
resp =  fit.model$y
mu =  fitted(fit.model)
fi =  fitted(fit.model$dispersion)
fi =  1/fi
R =  solve(t(Z)%*%Z)
R =  Z%*%R%*%t(Z)
r =  diag(R)
t =  resp*mu - 0.5*(mu*mu + resp*resp) 
t1 = t + 1/(2*fi)
tt = t1/(sqrt(1-r)*(1/(fi*sqrt(2))))
di = (r/(1-r))*(tt^2)
#
# Reta de corte com 3 desvios padrao
plot(di,xlab="Indice", ylab="Distancia de Cook",cex=2, cex.axis=1.5, 
cex.lab=1.5,pch=16)
cut = mean(di) + 3*sd(di)
abline(cut,0,lwd=2,lty=2)
# Comando para ser acionado com n sendo o numero de pontos destacados
# identify(di, n=1,cex=2)
#------------------------------------------------------------------------------------------#