substitua sua penultima linha por: x = dens[['x']][dens[['x']] > 1] y = dens[['y']][dens[['x']] > 1] lines(x, y, col=2) (mas vc precisa ver os detalhes sobre estimacao de densidade por kernel para entender pq o seu grafico original nao comeca em 1) b