Отображение 10-го и 90-го процентилей с помощью boxplot()

Question

Я пытаюсь использовать boxplot() функция в R, чтобы показать 10-й и 95-й процентили, а не 1,5*IQR по умолчанию. Люди предоставили решения на этом форуме, но ни один из них, кажется, не работает для меня.

Для всех моих первых нескольких участков я использовал ggplot2, который работал отлично. Однако недавно я обнаружил, что stat_summary наметил среднее значение в неподходящем месте для одного из сайтов, которые я изучаю летом. Изображение, сгенерированное ggplot2, я мог решить эту проблему только с помощью Boxplot() используя следующее;

 boxplot(mydata$tgm~interaction(mydata$site, mydata$season), 
         las = 1,
         names = c("BFa", "MBa", "STa", "BFsp", "MBsp", "STsp", "BFsu", "MBsu", "STsu", "BFw", "MBw", "STw"), 
         col = c("royalblue", "red", "green", "royalblue", "red", "green", "royalblue", "red", "green", "royalblue", "red", "green"), 
         outline = FALSE, 
         xlab = "Site per season", 
         ylab = "TGM concentration (ng/m3)")

А потом:

means <- aggregate(mydata$tgm~interaction(mydata$site, mydata$season), data = mydata, mean)

points(means, col = "black", pch = 16)

Могу ли я изменить свой аргумент так, чтобы отображались 10-й и 90-й процентили?