r - 计算经验概率

Question

感觉应该问这个问题，但搜索 stats.stackexchange 我找不到答案。

我在 R 中有这个数字向量:

x <- c(0.01, 0.50, 0.00, 0.28, 0.00, 0.15, 0.67, 0.00, 0.23, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.23, 0.00, 0.15, 0.00)

使用 R 代码，如何计算一系列值的概率？例如，某个值在 0 到 0.01 之间的概率是多少？

此外，使用 R 代码，某个值等于或小于 0.3 的概率是多少？

Answer 1

当您说“之间”时，您应该说明包含哪些端点(如果有)

有很多方法可以解决这些问题；我没有以有效的方式执行这些操作，只是举例说明了各种方法。

解决 P(a < X <= b) 类型问题的一种方法:

diff(ecdf(x)(c(0,0.01)))

对于 P(X <= b):

ecdf(x)(0.3)

另一种获取 P(X=x) 的方法:

table(x)/length(x)

并获取 P(X <= x):

cumsum(table(x))/length(x)

因此，如果您想包含两个端点，您可以这样做:

对于 P(a <= X <= b):

p <- table(x)/length(x)
diff(ecdf(x)(c(0,0.01)))+ p[names(p)=="0"]

或者:

p <- table(x)/length(x)
sum(p[names(p) %in% c("0","0.01")])

或者只是

length(x[0<=x & x<=0.01])/length(x)