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

 

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

Banco de Dados para Desenvolver Algoritmo de MLS para Previsão

DBO	ICobV	ICArb	Bcont	Dis_Pl	IBD_A
1,604	89	60	11	9	90
0,385	90	61	10	8,9	91
0,216	91	62	9	9,1	92
0,303	90	59	10	8,8	89
1,961	20	12	81	0,2	20
0,782	21	14	79	0,3	22
0,57	22	15	78	0,25	23
2,187	22	12	77	0,2	24
0,764	59	35	41	6	60
0,273	60	32	40	6,5	61
1,883	64	33	38	5,8	63
0,581	62	32	37	5,6	62
0,18	79	50	21	8,2	80
0,007	80	49	20	7,8	79
2,028	80	48	18	8,2	81
2,431	79	47	21	7,7	78

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

---------------

# 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