Acuidade <- read.table("D:\\Dados\\Acuidade.txt", header=TRUE)

X <- Acuidade$idade

Y <- Acuidade$tempo

attach(Acuidade)

lm(Y~X)

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

Syy  # Calcule

Sxy  # Calcule

n=length(X)

B1 <- # Calcule 

B0 <- # Calcule 

par(mfrow=c(1,2)) # Gráficos em duas colunas

plot(X,Y, pch=16) 
 
abline( lm (Y ~ X) ) # Reta ajustada pelo R com a função lm: linear model

plot(X,Y, pch=16)

abline(B0, B1) # Reta com ajuste manual