aves<-read.table("aves_cerrado.csv", header=T, sep=";")
head(aves)
str(aves)
aves$urubu[is.na(aves$urubu)] <- 0
aves$carcara[is.na(aves$carcara)] <- 0
aves$seriema[is.na(aves$seriema)] <- 0
aves$fisionomia[aves$fisionomia=="ce"] <- "Ce"
aves$fisionomia <- factor(aves$fisionomia, levels=c("CL","CC","Ce"))
str(aves)
head(aves)
#1.Ha relacao entre o numero de avistamentos de seriemas e carcaras em alguma das fisionomias amostradas?
modelce<-lm(seriema~carcara, data=aves[aves$fisionomia=="Ce",])
modelcc<-lm(seriema~carcara, data=aves[aves$fisionomia=="CC",])
modelcl<-lm(seriema~carcara, data=aves[aves$fisionomia=="CL",])
summary(modelce)
summary(modelcc)
summary(modelcl)
ce<-summary(modelce)$coefficients[2,4]
cc<-summary(modelcc)$coefficients[2,4]
cl<-summary(modelcl)$coefficients[2,4]
coef(modelce)
coef(modelcc)
coef(modelcl)
#Ha diferencas na relacao entre o numero de avistamentos de seriemas e carcaras entre as fisionomias do cerrado?
aves.dif<-lm(seriema~carcara*fisionomia, data=aves)
summary(aves.dif)
#Resposta: Nao ha diferencas significativas na relacao entre o numero de avistamentos de seriemas e carcaras entre as fisionomias do cerrado