Se tivesse a posteriori de um parâmetro dava para fazer um histograma/densityplot e pegar intervalos com a função quantile. Como o caso é bidimensional, tem-se que fazer uma densidade kernel em duas dimensões. Para isso veja o help da função::kde2d() (http://stat.ethz.ch/R-manual/R-devel/library/MASS/html/kde2d.html). Outra opção é usar KernSmooth::bkde2D(). Inclusive, dá para fazer os gráficos com contornos partindo desse código parcial da mensagem.
========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 skype: walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================