##rm(list=ls())
##
## Exemplo 1
##  - estimar a média geral
##  - estimar a média por grupo
df1 <- read.csv("http://www.leg.ufpr.br/ce092/dados/ModMix1.csv", header = TRUE)
head(df1)
df1

with(df1, plot(y ~grupo, pch = 20))

##
## Exemplo 2
##  - Descrever a relação entre x e y
rm(list=ls())
df2 <- read.csv("http://www.leg.ufpr.br/ce092/dados/ModMix2.csv", header = TRUE)

with(df2, plot(x, y))
with(df2, cor(x, y))

## alguns ajustes "manuais"
abline(c(52.72, 1.287), col = "blue", lty = 2)
abline(c(81.14, -0.879), col = "red", lty = 2)
## Qual o melhor ajuste? ... será?