# ===== Motivação ======

acuidade <- read.table("F:\\Dados\\acuidade.dat", header=TRUE)

X<-acuidade$idade

Y<-acuidade$tempo

lm(Y~X)

plot(X,Y, pch=16)

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


# ===== Feminino ======

sexo<-acuidade$sexo

X.fem<-acuidade$idade[which(sexo=="F")]

Y.fem<-acuidade$tempo[which(sexo=="F")]

lm(Y.fem~X.fem)

plot(X.fem,Y.fem, pch=16)

abline( lm (Y.fem ~ X.fem) ) # Reta ajustada pelo R com a função lm: linear model


# ===== Masculino ======

X.masc<-acuidade$idade[which(sexo=="M")]

Y.masc<-acuidade$tempo[which(sexo=="M")]

lm(Y.masc~X.masc)

plot(X.masc,Y.masc, pch=16)

abline( lm (Y.masc ~ X.masc) ) # Reta ajustada pelo R com a função lm: linear model


# ===== Feminino e Masculino lado a lado ======

par(mfrow=c(1,2))

plot(X.fem,Y.fem, pch=16)

abline( lm (Y.fem ~ X.fem) ) # Reta ajustada pelo R com a função lm: linear model

plot(X.masc,Y.masc, pch=16)

abline( lm (Y.masc ~ X.masc) ) # Reta ajustada pelo R com a função lm: linear model



# ===== Feminino e Masculino juntos em cores ======

par(mfrow=c(1,1))

faixaX <- range(X.fem, X.masc)

faixaY <- range(Y.fem, Y.masc)



plot(X.fem,Y.fem, pch=16, ylim=faixaY, xlim=faixaX, col="red", ylab="", xlab="")

abline( lm (Y.fem ~ X.fem), col="red") 

par(new=T)

plot(X.masc,Y.masc, pch=16, ylim=faixaY, xlim=faixaX, col="blue", ylab="tempo de reação", xlab="idade")

abline( lm (Y.masc ~ X.masc), col="blue") 

leg.names <- c("Feminino","Masculino")
 legend(locator(1), leg.names, pch=16, col=c("red", "blue"))