Outlier(d.pizza$temperature, na.rm=TRUE)
# find the corresponding rows
d.pizza[which(d.pizza$temperature %in% Outlier(d.pizza$temperature, na.rm=TRUE)),]
# outliers for the drivers
tapply(d.pizza$temperature, d.pizza$driver, Outlier, na.rm=TRUE)
# see also
boxplot(temperature ~ driver, d.pizza)$out
Run the code above in your browser using DataLab