# --------------- Dados 1 - acuidade visual ---------------------# 
# --------------- Leitura dos dados  ---------------# 

#Dados acuidade visual
dados<-read.table("http://wiki.icmc.usp.br/images/0/0f/Acuidade.txt", header=TRUE)
attach(dados)
names(dados)
Y<-acuidade
X<-idade


# ------------- (i) Ajuste do modelo ---------------# 

 modelo.ajustado <- lm(Y~X)

 modelo.ajustado

 plot(X,Y, pch=16, xlab="Idade", ylab="Tempo de reação", main="Modelo de regressão linear simples") 
 
 abline(modelo.ajustado) 
 
 summary(modelo.ajustado)


# ------------- (ii) Tabela ANOVA -----------------# 

 anova(modelo.ajustado)


## ---------- Tabela ANOVA com erro puro e falta de ajuste

 # install.packages("alr3") ## Executar somente uma vez para instalar a biblioteca

 library(alr3)

 pureErrorAnova(modelo.ajustado)


# ------------- (iii) Análise de Resíduos ---------# 

 y.ajustado <- fitted(modelo.ajustado)

 resíduos <- residuals(modelo.ajustado)

 plot(y.ajustado, resíduos)

 plot(y.ajustado, resíduos, pch=16, xlab=expression(hat(Y)), ylab=expression(e[i]), main="Gráfico de resíduos")
 
 plot(y.ajustado, resíduos, pch=16, xlab=expression(hat(Y)), ylab=expression(Y-hat(Y)), main="Gráfico de resíduos")

 n<-length(Y)

 sigma.chapéu <- sqrt(sum(modelo.ajustado$residuals^2)/(n-2))
 
 resíduos.padronizados <- resíduos/sigma.chapéu

 plot(y.ajustado, resíduos.padronizados)

 faixay <- range(c(-3,3),resíduos.padronizados)

 plot(y.ajustado, resíduos.padronizados, pch=16, ylim=faixay, xlab=expression(hat(Y)), ylab=expression(d[i]), main="Gráfico de resíduos padronizados")

  abline(qnorm(0.975),0, col=2, lty=2)

  abline(qnorm(0.025),0, col=2, lty=2)


# -------------  Gráfico de probabilidade normal ---------# 

 qqnorm(resíduos.padronizados )		


# ------------- (ii) Envelope ---------# 
# ---------- Baixe o envel_norm.txt de www.ime.usp.br/~giapaula/envel_norm 

source("http://wiki.icmc.usp.br/images/8/82/Envel_norm.txt") 


# --------------- Dados 2 - tratores ---------------------# 

#Dados tratores
dados<-read.table("http://wiki.icmc.usp.br/images/5/51/Tratores.txt", header=TRUE)
attach(dados)
names(dados)
Y<-
X<-


# --------------- Dados 3 - foguetes  ---------------------# 

#Dados foguetes
dados<-read.table("http://wiki.icmc.usp.br/images/0/06/Foguetes.txt", header=TRUE)
attach(dados)
names(dados)
Y<-
X<-