#------------------------------------------------------------------------------------------#
# Comandos:
#          fit.model = ajuste
#          attach(dados)
#          source("diag_cook_norm_dglm.txt")
#------------------------------------------------------------------------------------------#
X =  model.matrix(fit.model)
n =  nrow(X)
p =  ncol(X)
require(dglm)
fi =  fitted(fit.model$dispersion)
fi =  1/fi
for(i in 1:n){
X[i,] = fi[i]*X[i,]
}
H =  X%*%solve(t(X)%*%X)%*%t(X)
h =  diag(H)
r =  resid(fit.model)
ts =  (sqrt(fi)*r)/sqrt(1-h)
di =  (h/(1-h))*(ts^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,lty=2,lwd=2)
# Comando para ser acionado com n sendo o numero de pontos destacados
# identify(di, n=1, cex=2)
#------------------------------------------------------------------------------------------#