Write a program in R that fits and interprets a regression model to the smiles.csv data. Your program should do the following: 1. Produce a boxplot of leniency score for the various types of smiles. 2. Fit the regression model that treats "neutral" as the reference level of type. Display the model summary. 3. Report the ANOVA table with the information for testing whether mean leniency differs by smile type. 4. Estimate the population mean leniency scores by smile type, along with their standard errors. Display the results. 5. Estimate all pairwise differences in mean leniency score among smile type, along with the standard errors of the differences. Display the results. 6. Produce an effects plot of the effect of the levels of smile type on leniency score