Se os contornos de confiança são para algo que se assume ter distribuição normal, pode-se usar a função ellipse() do pacote ellipse. Você vai precisar do vetor de médias e matriz de covariâncias. Veja a documentação para maiores detalhes.

À disposição.

Walmes.