Чтобы добавить легенду к ответу joran в базе R:
legend("topleft",title="Decile",legend=c(1:10),col =rbPal(10),pch=20)
В этом примере добавляется «, cex = 0.8» только для красивости: