fab.dat=scan("e:/UEPB1/fabrica.dat",what=list(falha=0,mes=0))
attach(fab.dat)
y=fab.dat$falha
x=fab.dat$mes
fit=glm(y~x,family=poisson())
summary(fit)
fit0=glm(y~1,family=poisson())
anova(fit0,fit,test="Chi")
yh=fitted(fit)
X=model.matrix(fit)
w=fit$weights
W=diag(w)
H=sqrt(W)%*%X%*%solve(t(X)%*%W%*%X)%*%t(X)%*%sqrt(W)
h=diag(H)
plot(h,xlab="Índice")
identify(h)
ro=residuals(fit,type="response")
fi=1
rd=residuals(fit,type="deviance")
td=rd*sqrt(fi/(1-h))
rp=residuals(fit,type="pearson")
rp=sqrt(fi)*rp
ts=rp/sqrt(1-h)
LD=h*(ts^2)/(1-h)
envelope.poisson(form=y~x,Fam=poisson(),k=100,alfa=0.05)
plot(td,ylab="Componente do desvio",xlab="Índice")
identify(td)
plot(LD,ylab="Distância de Cook",xlab="Índice")
identify(LD)