obrigado por se oferecer a ajudar. Informo que encontrei a solução e compartilho com os demais integrantes da lista r-br.
Ressalto que o gráfico exibe pontos com cores e tamanhos segundo o quartil do dado observado.
quantile(dadosT1$data)
points(dadosT1,xlab="Coordenada X (coluna)",ylab="Coordenada Y (linha)", pt.divide="quartiles")
legend(270,70, c("1º quartil", "2º quartil", "3º quartil", "4º quartil"),
col=c('black', 'black', 'black', 'black'), pt.bg=c('blue', 'green', 'yellow', 'red'), pch=c(21,21,21,21),
pt.cex=c(0.6,0.9,1.3,1.6),bty='0')