# Ajuste de modelo de regressão linear simples via Estimadores de Mínimos Quadrados

Acuidade = read.table("F:\\Pastas\\Dropbox\\USP\\SME0122 - Int. a Inf. Estatistica\\Dados - tempo de reacao.txt", header=TRUE)

attach(Acuidade)  # Anexa o conjunto de dados e pode-se acessar as variáveis direto pelo nome

Y <- tempo

X <- idade

n <- length(X)

Sxx <- sum(X^2) - n * mean(X)^2

Sxy <- sum(X*Y) - n * mean(X) * mean(Y)

b <- Sxy/Sxx

a <- mean(Y) - b*mean(X)

plot(X,Y, pch=16)

abline(a, b, col=2, lwd=2) # Reta com ajuste manual 

legend(locator(1), paste("Reta ajustada: Y=", a, "+", b, "X")) # Adicionando legenda



# Alternativamente, usando direto a função lm (linear model)

plot(X,Y, pch=16)

abline( lm (Y ~ X), col=4, lwd=2 ) # Reta ajustada pelo R com a função lm: linear model

legend(locator(1), paste("Reta ajustada: Y=", a, "+", b, "X")) # Adicionando legenda