In [ ]:
# carga de los datos
insurance <- read.csv("insurance.csv", stringsAsFactors = TRUE)
str(insurance)
In [ ]:
summary(insurance$expenses)
In [ ]:
hist(insurance$expenses)
In [ ]:
table(insurance$region)
In [ ]:
cor(insurance[c("age", "bmi", "children", "expenses")])
In [ ]:
pairs(insurance[c("age", "bmi", "children", "expenses")])
In [ ]:
pairs.panels(insurance[c("age", "bmi", "children", "expenses")])
In [ ]:
ins_model <- lm(expenses ~ age + children + bmi + sex + smoker + region,
data = insurance)
In [ ]:
ins_model <- lm(expenses ~ ., data = insurance)
In [ ]:
ins_model
In [ ]:
In [ ]:
In [ ]:
summary(ins_model)
In [ ]:
insurance$age2 <- insurance$age^2
In [ ]:
insurance$bmi30 <- ifelse(insurance$bmi >= 30, 1, 0)
In [ ]:
ins_model2 <- lm(expenses ~ age + age2 + children + bmi + sex + bmi30*smoker + region,
data = insurance)
summary(ins_model2)
In [ ]:
In [ ]: