##Exercício 7##

#7.2#
library(MASS)
data(Animals)
anim.m2<-lm(log(brain)~log(body),data=Animals,subset=!(log(Animals$body)>8&log(Animals$brain)<6))
#Ajuste do modelo
anim.m0<-lm(log(brain)~1, data=Animals,subset=!(log(Animals$body)>8&log(Animals$brain)<6))
#Compare dos 2:
anova(anim.m0,anim.m2)

#1. Qual a relação do comando 'anova' acima com:
anova(anim.m2)
# O primeiro faz a ANOVA entre os dois modelos informados, enquanto o segundo compara o modelo informado (anim.2) com o modelo nulo. Ambos geram as mesmas informações.

#2. Qual a relação entre os valores obtidos por estes comandos:
summary(anim.m0)
mean(log(Animals$brain[!(log(Animals$body)>8&log(Animals$brain)<6)]))
sd(log(Animals$brain[!(log(Animals$body)>8&log(Animals$brain)<6)]))
# o comando "mean" fornece a média, mesmo valor da estimativa do desvio padrão fornecida pelo comando summary
# o comando "sd" fornece o desvio padrão, que é o mesmo valor fornecido pelo resíduo do desvio padrão fornecido pelo summary.


#7.3#



aves<-read.table("aves_cerrado.csv",header=TRUE,sep=";")
str(aves_cerrado)
aves$carcara[is.na(aves$carcara)]<-0
aves$seriema[is.na(aves$seriema)]<-0
aves$urubu[is.na(aves$urubu)] <- 0
aves

#modelo linear
mod.ce<-lm(seriema~carcara,data=aves,subset=aves$fisionomia=="Ce")
mod.cc<-lm(seriema~carcara,data=aves,subset=aves$fisionomia=="CC")
mod.cl<-lm(seriema~carcara,data=aves,subset=aves$fisionomia=="CL")

#gráficos
plot(seriema~carcara,data=aves)
abline(mod.cc, col="blue")
abline(mod.ce,col="red")
abline(mod.cl,col="green")
par(mfrow=c(3,1))
plot(seriema~carcara,data=aves,subset=aves$fisionomia=="Ce")
abline(mod.ce,col="red")
plot(seriema~carcara,data=aves,subset=aves$fisionomia=="CC")
abline(mod.cc, col="blue")
plot(seriema~carcara,data=aves,subset=aves$fisionomia=="CL")
abline(mod.cl,col="green")
par(mfrow=(1,1))

#Coeficientes
coef.ce<-coef(mod.ce) # 7.3921076  -0.1458764 
coef.cc <-coef(mod.cc) # 5.67584746 -0.08262712 
coef.cl<-coef(mod.cl) #7.5833333  -0.4507576 

#anova
anova(mod.ce)
anova(mod.cc)
anova(mod.cl)

p.ce<-0.2401
p.cc<-0.6454
p.cl<-0.1045

# 1.Há relação entre o número de avistamentos de seriemas e carcarás em alguma das fisionomias amostradas?
 Não há relação significativa entre o avistamento de carcarás e de seriemas em nenhuma das regiões do cerrado.

# 2.Há diferenças na relação entre o número de avistamentos de seriemas e carcarás entre as fisionomias do Cerrado?
 Há uma relação maior na área cl do que nas áreas cc e ce, pois possui um maior valor de P.