estou usandoecdf() função para calcular a função empírica da densidade cumulativa (ECDF) de algumas amostras aleatórias:

set.seed(0)
X = rnorm(100)
P = ecdf(X)

AgoraP fornece ECDF e podemos plotá-lo:

plot(P)
abline(h = 0.6, lty = 3)

Minha pergunta é: como posso encontrar o valor da amostrax, de tal modo queP(x) = 0.6, ou seja, o quantil 0,6 do ECDF ou a coordenada x do ponto de interseção do ECDF eh = 0.6?