library(robustbase) adjbox(Final$Temp[Final$Loc=="UK"]~Final$Tod[Final$Loc=="UK"]) Boxplots are a visualization tool, so allow mei'll give you a visual advice. What you have is essentially functional data so you want (for visualization reasons) to use a box-plot tool that acknowledges that. Try the functional boxplot function in the fda package.