dat <- read.delim("clipboard"); attach(dat) #importar los datos del portapapeles
names(dat) #Chequear los nombres
str(dat) #Chequear la estructura, el tratamiento y el sexo son caracteres, hay que volverlos factores
dat$Uso=as.factor(dat$Uso)
dat$Especie=as.factor(dat$Especie)
str(dat)
dat$Presencia

model.results = glm(Presencia ~ Especie*Uso*Abundancia.plantas, binomial(link = "probit"))
anova(model.results,test="Chisq")
summary(model.results)
8289/15708

model.results1 = glm(Presencia ~ Especie*Uso+Especie*Abundancia.plantas+Uso*Abundancia.plantas, binomial(link = "probit"))
anova(model.results1,test="Chisq")
anova(model.results1,model.results,test="Chisq")

model.results2 = glm(Presencia ~ Especie*Uso+Especie*Abundancia.plantas+Uso, binomial(link = "probit"))
anova(model.results2,test="Chisq")
anova(model.results2,model.results1,test="Chisq")


model.results4 = glm(Presencia ~ Especie*Uso+Especie+Abundancia.plantas+Uso, binomial(link = "probit"))
anova(model.results4,test="Chisq")

model.results5 = glm(Presencia ~ Especie*Uso+Especie+Uso, binomial(link = "probit"))
anova(model.results5,test="Chisq")

anova(model.results5,model.results4,test="Chisq")
anova(model.results5,model.results2,test="Chisq")

summary(model.results5)

8341/15771

library(DHARMa)
residh <- simulateResiduals(model.results5)
plot(residh)

model.resultsqb = glm(Presencia ~ Especie*Uso+Especie+Uso, quasibinomial(link = "probit"))
anova(model.resultsqb,test="F")



#Poder explicativo
dev.null <- model.results5$null.deviance
dev.resid <- model.results5$deviance
dev.explained <- (dev.null-dev.resid)/dev.null
dev.explained

#Poder explicativo
dev.null1 <- model.resultsqb$null.deviance
dev.resid1 <- model.resultsqb$deviance
dev.explained1 <- (dev.null1-dev.resid1)/dev.null1
dev.explained1