For the following questions you work with the built-in dataset Puromycin. It contains reaction velocity of an enzymatic reaction at different concentrations, and for two different states (treated and untreated):

 

                
               conc rate   state
1  0.02   76 treated
2  0.02   47 treated
3  0.06   97 treated
4  0.06  107 treated
5  0.11  123 treated
6  0.11  139 treated
7  0.22  159 treated
8  0.22  152 treated
9  0.56  191 treated
10 0.56  201 treated

 

Question: Use base graphics to make a boxplot of the rates for each state. Fill the boxplot with a grey color.

 

Note: The command as.grob() from the ggplotify package enables Dodona to check your result. Just put your ordinary graphical command(s) between the two comments.