Untitled

# Create a boxplot for each UPF tertile group
boxplot(scores ~ group, data = data,
        col = c("lightblue", "lightgreen", "lightpink"),
        main = "Verbal Scores by UPF Consumption Tertile",
        xlab = "UPF Consumption Tertile",
        ylab = "Verbal Score (McCarthy Scale)")