# TIC II · S27
coches <- read.csv("93cars_DATOS.csv", fileEncoding = "UTF-8-BOM")
stopifnot(nrow(coches)==93, "Mid.Price" %in% names(coches))
precio <- coches$Mid.Price
potencia <- coches$Horsepower
validos <- complete.cases(precio,potencia)
r <- cor(potencia[validos],precio[validos])
modelo <- lm(precio[validos] ~ potencia[validos])
cat("Coeficiente r:",r,"\n")
print(coef(modelo))
plot(potencia[validos],precio[validos],pch=19, xlab="Potencia (HP)", ylab="Precio (miles de dólares)",main="Potencia y precio")
abline(modelo,col="red",lwd=2)
# RETO: ¿r es positivo o negativo? ¿Qué significa?