```
Outlier(d.pizza$temperature, na.rm=TRUE)
# it's the same as the result from boxplot
sort(d.pizza$temperature[Outlier(d.pizza$temperature, value=FALSE, na.rm=TRUE)])
b <- boxplot(d.pizza$temperature, plot=FALSE)
sort(b$out)
# nice to find the corresponding rows
d.pizza[Outlier(d.pizza$temperature, value=FALSE, na.rm=TRUE), ]
# compare to Hampel's rule
Outlier(d.pizza$temperature, method="hampel", na.rm=TRUE)
# outliers for the each driver
tapply(d.pizza$temperature, d.pizza$driver, Outlier, na.rm=TRUE)
# the same as:
boxplot(temperature ~ driver, d.pizza)$out
```

Run the code above in your browser using DataLab