###Exercícios - Aula 7

getwd()
setwd("C:/Users/Edmundo/Documents/Meus Documentos/BIANCA/MESTRADO/Linguagem R")
getwd()

##Altura na Infância e na Vida Adulta

crianca <- c(39,30,32,34,35,36,36,30)
adulto <- c(71,63,63,67,68,68,70,64)
plot(adulto~crianca)
crescimento1 <- lm(adulto~crianca)
abline(crescimento1,col="blue")
anova(crescimento1)
#P=0.000386. P<0.05 = há relação significativa entre a altura com dois anos e a altura do adulto.
coef(crescimento1)
confint(crescimento1)
#Valores esperados:
crescimento2 <- lm(2*crianca~crianca)
abline(crescimento2,col="red")
anova(crescimento2)
#
coef(crescimento2)
confint(crescimento2)
#
#Os dados não corroboram a hipótese do senso comum porque os limites de confiança do primeiro modelo não abrangem os coeficientes do segundo modelo.

##Seriemas e Carcarás

aves <- read.table("aves_cerrado.csv",header=TRUE,sep=";",as.is=TRUE)
aves
names(aves)
list(aves$fisionomia)
aves$fisionomia[18]<-"Ce"
aves$fisionomia
is.na(aves)
aves$urubu[is.na(aves$urubu)] <- 0
aves$seriema[is.na(aves$seriema)] <- 0
aves$carcara[is.na(aves$carcara)] <- 0
is.na(aves)
aves
#
avesCC <- aves[aves$fisionomia=="CC",]
avesCC
avesCL <- aves[aves$fisionomia=="CL",]
avesCL
avesCe <- aves[aves$fisionomia=="Ce",]
avesCe

#1.Há relação entre o número de avistamentos das aves em alguma das fisionomias amostradas?
plot(aves$seriema~aves$carcara)
avistamentos1 <- lm(seriema~carcara,data=avesCC)
avistamentos1
abline(avistamentos1,col="red")
avistamentos2 <- lm(seriema~carcara,data=avesCL)
avistamentos2
abline(avistamentos2,col="blue")
avistamentos3 <- lm(seriema~carcara,data=avesCe)
avistamentos3
abline(avistamentos3)
anova(avistamentos1)
#P=0.6454. P>0.05 = há relação significativa.
anova(avistamentos2)
#P=0.1045. P>0.05 = há relação significativa.
anova(avistamentos3)
#P=0.2446. P>0.05 = há relação significativa.

#2.Há diferenças na relação entre o número de avistamentos de seriemas e carcarás entre as fisionomias de Cerrado?
difseriema <- aov(seriema~factor(fisionomia),data=aves)
summary(difseriema)
#P=0.518. P>0.05 = não significativo = não existe diferença para avistamento de seriema.
difcarcara <- aov(carcara~factor(fisionomia),data=aves)
summary(difcarcara)
#P=3.127e-14. P<0.05 = significativo = existe diferença para avistamento de carcará.

##Resíduos de Iris

iris
#1.
plot(iris$Sepal.Length~iris$Sepal.Width)
questao1 <- lm(Sepal.Length~Sepal.Width,data=iris)
abline(questao1)
anova(questao1)
#P=0.1519. P>0.05 = não signicativo = não existe relação entre a largura da sépala e o comprimento da mesma.
#2.
#Não entendi o objetivo da questão.