data(therapy)
plot(therapy ~ perstest, data=therapy, pch=16)
abline(lm(therapy ~ perstest, data=therapy), col="red")
plot(therapy ~ perstest, data=therapy, cex=1.5, pch=16,
col=ifelse(sex=="M", "red","blue"))
Run the code above in your browser using DataLab