## Exercícios de Regressões Lineares Simples
### Ana Carolina Luchetta

##Altura na Infância e Na Vida Adulta
jovem<-c(39,30,32,34,35,36,36,30)
adulto<-c(71,63,63,67,68,68,70,64)
relação.altura <- lm(adulto~jovem)
summary(relação.altura)                       #Como o p é inferior a 5%, podemos dizer que há uma relação significativa entre a altura do jovem e do adulto.

plot(adulto~jovem)
abline(relação.altura)
esperado <- 2*jovem
abline(lm(esperado~jovem), col=5)
confint(relação.altura)
esperado                                       #Os dados esperados superam o intervalo de confiança. Pela reta traçada no gráfico já era possível imaginar que a tese não seria corroborada.

##Seriemas e Carcarás

aves                                           #percebemos que temos muitos NAs na tabela, e vamos substituir por 0.
aves$urubu[is.na(aves$urubu)] <- 0
aves$carcara[is.na(aves$carcara)] <- 0
aves$seriema[is.na(aves$seriema)] <- 0
aves
unique (aves$fisionomia)                        #percebemos que tem um ce que precisa ser alterado para Ce
aves$fisionomia[aves$fisionomia=="ce"] <- "Ce"  #substituição
aves                                            #agora podemos começar a trabalhar
class (aves$fisionomia)
aves$fisionomia <- factor(aves$fisionomia, levels=c("Ce", "Cc", "Cl"))

aves.ce <- aves[aves$fisionomia=="Ce",]
aves.cc <- aves[aves$fisionomia=="Cc",]
aves.cl <- aves[aves$fisionomia=="Cl",]

rel.ce <- lm (aves.ce$seriema ~ aves.ce$carcara)
rel.cc <- lm (aves.cc$seriema ~ aves.cc$carcara)
rel.cl <- lm (aves.cl$seriema ~ aves.cl$carcara)
summary(rel.ce)						#não há diferença para a localidade Ce, devido ao p baixo
summary(rel.cc)						#não há diferença para a localidade Cc, pois o p é muito baixo
summary(rel.cl)						#não há diferença para a localidade Cl também, deivido ao p.

par(mfrow=c(2,2))
plot(aves.ce$seriema ~ aves.ce$carcara)
abline(lm(aves.ce$seriema ~ aves.ce$carcara))

plot(aves.cc$seriema ~ aves.cc$carcara)
abline(lm(aves.cc$seriema ~ aves.cc$carcara))

plot(aves.cl$seriema ~ aves.cl$carcara)
abline(lm(aves.cl$seriema ~ aves.cl$carcara))




##Resíduos íris

require(datasets)
head(iris)
iris_setosa<- iris[iris$Species=="setosa",]
str(iris_setosa)
setosa_sepala<-lm(iris_setosa$Sepal.Width~iris_setosa$Sepal.Length)
plot(iris_setosa$Sepal.Width~iris_setosa$Sepal.Length)
abline(setosa_sepala)
summary(setosa_sepala)

setosa_petala<-lm(iris_setosa$Petal.Width~iris_setosa$Petal.Length)
plot(iris_setosa$Petal.Width~iris_setosa$Petal.Length)
abline(setosa_petala)
residuo<- residuals(setosa_petala)
residuo
setosa_petala2<-lm(iris_setosa$Petal.Width~residuo)
plot(iris_setosa$Petal.Width~residuo)
abline(setosa_petala2)
summary(setosa_petala)
summary(setosa_petala2)