Segunda aula do R dada pelo Professor Cristian Villegas do LCE/ESALQ

 

#------------------------------------------------

# DBO Demanda Bioquímica de Oxigênio

# ICobV Índice de Cobertura Vegetal

# ICArb Índice de Cobertura Arbórea

# IBCont Bioindicador de Contaminação (agrotóxicos)

# Dis_Pl Distancia do plantio de grãos.

#------------------------------------------------

library(readxl)

dados <- read_xlsx("dados_regressao_multipla.xlsx")

str(dados)

#------------------------------------------------

# Ajuste modelo 1

#------------------------------------------------

modelo1<- lm(IBD_A ~ DBO+ICobV+ICArb+Bcont+Dis_Pl,

             data = dados)

summary(modelo1)

resumo_modelo1<- summary(modelo1)

names(resumo_modelo1)

resumo_modelo1$coefficients

resumo_modelo1$sigma

resumo_modelo1$r.squared

resumo_modelo1$adj.r.squared

resumo_modelo1$fstatistic

vcov(modelo1)

sigma(modelo1)

par(mfrow=c(2,2))

plot(modelo1, pch = 20)

#------------------------------------------------

# Ajuste modelo 2

#------------------------------------------------

modelo2<- lm(IBD_A ~ ICArb+Bcont+Dis_Pl,

             data = dados)

summary(modelo2)

resumo_modelo2<- summary(modelo2)

names(resumo_modelo2)

resumo_modelo2$coefficients

resumo_modelo2$sigma

resumo_modelo2$r.squared

resumo_modelo2$adj.r.squared

resumo_modelo2$fstatistic

vcov(modelo2)

sigma(modelo2)

par(mfrow=c(2,2))

plot(modelo2, pch = 20)

#------------------------------------------------

# Alternativas para estudar depois

#------------------------------------------------

modelo1_opcao<- lm(IBD_A ~ ., data = dados)

summary(modelo1_opcao)

summary(modelo1)

modelo2_opcao<- update(modelo1_opcao, 

                       .~ICArb+Bcont+Dis_Pl,

                       data = dados)

modelo2_opcao2<- update(modelo1_opcao, 

                        .~.-DBO-ICobV,

                        data = dados)

summary(modelo2)

summary(modelo2_opcao)

summary(modelo2_opcao2)

Atenciosamente

Cristian Villegas

http://lattes.cnpq.br/7868115035277497
https://orcid.org/0000-0003-3176-5236

University of São Paulo, Department of Exact Science

Avenida Pádua Dias, 11 - Piracicaba/SP - CEP 13418-900
+551934478912