我有一个观察载体，想用R获得每个obervation的经验p值。我不知道潜在的分布，我现在做的只是
runif(100,0,1000)->ay
quantile(ay)
然而，这并没有给我一个p值。如何凭经验获得p值？
答案 0 :(得分：1)
我认为这是你正在寻找的东西：
rank(ay)/length(ay)
答案 1 :(得分：1)
我认为你想要的是ecdf功能。这将返回经验累积分布函数，您可以直接应用
ay <- runif(100)
aycdf <- ecdf(ay)
然后
> aycdf(c(.1, .5, .7))
[1] 0.09 0.51 0.73