#Exemplo 1.10.3 do livro Modelos de regressao com apoio computacional. Paula, G.A. (2004)
num<-c(175,108,95,82,71,50,49,31,28,17,16,11)
tempo<-1:12
win.graph()

#Modelo Linear-Y
fit<-lm(num~tempo)
summary(fit)
form<-num~tempo
k<-100
alfa<-0.05
envelope.normal(form=form,k=k,alfa=alfa)
#Modelo quadratico-Y
tempo2<-tempo^2
fit2<-lm(num~tempo+tempo2)
summary(fit2)
form<-num~tempo
k<-100
alfa<-0.05
envelope.normal(form=form,k=k,alfa=alfa)
#Modelo Linear-sqrt(Y)
sqnum<-sqrt(num)
fit3<-lm(sqnum~tempo)
summary(fit3)
form<-sqnum~tempo
k<-100
alfa=0.05
envelope.normal(form=form,k=k,alfa=alfa)
#Modelo quadratico-sqrt(Y)
fit4<-lm(sqnum~tempo+tempo2)
summary(fit4)
form<-sqnum~tempo
k<-100
alfa<-0.05
envelope.normal(form=form,k=k,alfa=alfa)

#Modelo Poisson
fit5<-glm(num~tempo,family=Poisson())
summary(fit)
form<-num~tempo
k<-100
alfa<-0.05
Fam<-poisson()
envelope.poisson(form = form, Fam = Fam, k = k, alfa = alfa)