Легенда с разными размерами символов в базе R

В моей диаграмме я кодирую некоторую информацию в диаметрах нарисованных кругов. У меня вопрос: как проще всего задокументировать эту информацию в легенде?

Вот то, что я пробовал до сих пор:

dat <- rnorm(100)
cex_brks <- quantile(dat, c(0.25,0.5,0.75))
cex_size <- c(1,1.4,1.8, 2.2) 
cex <- rep(NA, length(dat))
for (i in 1:3) {
    cex[is.na(cex) & dat<=cex_brks[[i]]] <- cex_size[[i]]
}
cex[is.na(cex)] <- cex_size[[4]]
plot(dat, cex=cex, pch=21)
legend(
    "bottom", 
    legend=c("very small", "small", "large", "very large"), 
    bty="n",
    pch=21,
    cex=cex_size
)

Однако, делая это таким образом, изменяется не только размер символа (pch), но и но и текст легенды тоже. Как я могу переопределить это так, чтобы только символы легенды имели разные размеры?