 = cond)) +
  geom_boxplot()

# Replot without the outliers
mariokart |&gt;
  filter(total_pr &lt; 80) |&gt;
  ggplot(aes(x = total_pr, y = cond)) +
  geom_boxplot()

# Fit a multiple regression models
mariokart_no &lt;- mariokart |&gt; filter(total_pr &lt; 80)
m1 &lt;- lm(total_pr ~ cond + stock_photo + duration + wheels, data = mariokart_no)
tidy(m1)
m2 &lt;- lm(total_pr ~ cond + stock_photo + wheels, data = mariokart_no)
tidy(m2)
m3 &lt;- lm(total_pr ~ cond + wheels, data = mariokart_no)
tidy(m3)

# Fit 