RTutor

library(ggplot2)

ggplot(data = df, aes(x = type, y = highway, color = type)) +

geom_boxplot() +

geom_jitter(width = 0.2, alpha = 0.5)

library(plotly)

plot_ly(df, x = ~dis, y = ~city, z = ~highway, color = ~type, type = "scatter3d")

# ANOVA of log-transformed highway by type and drive

# Create a new data frame with the necessary variables

df_anova <- df %>%

select(type, drive, highway)

# Log-transform the highway variable

df_anova$highway <- log(df_anova$highway)

# Conduct the ANOVA

anova <- aov(highway ~ type + drive, data = df_anova)

# Print the results

summary(anova)

Df Sum Sq Mean Sq F value Pr(>F)

type 6 11.477 1.9128 115.82 < 2e-16 ***

drive 2 0.606 0.3028 18.34 4.2e-08 ***

Residuals 225 3.716 0.0165

---

Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1