Regressão Linear Simples

##-------------------------------------------------------------##
##               INSERINDO O BANCO DE DADOS                    ##
##-------------------------------------------------------------##

banco <- read.csv2(file.choose(), sep = ";", dec = ",", header = T)
attach(banco)
banco

##     N Resistencia Peso
## 1   1         265 26.0
## 2   2         197 22.5
## 3   3         346 29.0
## 4   4         280 28.0
## 5   5         265 26.5
## 6   6         200 23.0
## 7   7         221 23.5
## 8   8         265 24.5
## 9   9         261 26.0
## 10 10         278 29.0
## 11 11         205 24.0
## 12 12         286 28.5
## 13 13         317 28.5
## 14 14         242 26.0
## 15 15         254 25.5
## 16 16         235 24.5
## 17 17         176 21.5
## 18 18         262 24.5
## 19 19         248 26.0
## 20 20         250 25.5
## 21 21         263 26.5
## 22 22         274 27.5
## 23 23         242 25.0
## 24 24         260 25.5
## 25 25         281 26.5

##-------------------------------------------------------------##
##               CRIANDO A MATRIZ DOS DADOS                    ##
##-------------------------------------------------------------##

X <- matrix(c(rep(1,25), Peso),25,2)
Y <- matrix(c(Resistencia), 25,1)

Coef <- solve(t(X)%*%X)%*%t(X)%*%Y

Yest <- X%*%Coef

##-------------------------------------------------------------##
##          GRÁFICO DE REGRESSÃO LINEAR SIMPLES                ##
##-------------------------------------------------------------##

plot(Resistencia~Peso, xlab = "Peso", ylab = "Força de Resistência à Ruptura",
     type = "p", pch = 21, bg = "yellow", las = 1)
abline(Coef, col = "red")

legend(22,350, c("Valores Observados", "Valor Estimado"), col = c("black", "red"),
       pch = c(21, NA), lty = c(1, 1), pt.bg = c("yellow", NA), bty = "n")