Распределение эмпирическое в R

Question

У меня есть вектор наблюдений, и я хотел бы получить эмпирическое значение р для каждой обервации с R. Я не знаю, какое распределение лежит в основе, и что я сейчас делаю, это просто

runif(100,0,1000)->ay
quantile(ay)

Тем не менее, это на самом деле не дает мне р-значение. Как я могу получить значение р опытным путем?

Answer 1

Я думаю, это то, что вы ищете:

rank(ay)/length(ay)

Answer 2

Я думаю, что вы хотите, это ecdf функция. Это возвращает эмпирическую интегральную функцию распределения, которую вы можете применить напрямую

ay <- runif(100)
aycdf <- ecdf(ay)

А потом

> aycdf(c(.1, .5, .7))
[1] 0.09 0.51 0.73