找出偏态分布的概率

Question

我有一个变量X，下面给出了摘要。看看分布情况似乎偏向正态分布。我想知道如何在R？

中找到与变量X的每个值相关的p值

very hardworking

Answer 1

根据您的情节判断，我假设您有一些名为averageRank的数字序列，您试图使用density.default()拟合内核密度估算器。

基本上，您对概率密度函数（density.default()的输出）进行了估计。你的问题陈述：

  

我想知道如何在R？

中找到与变量X的每个值相关联的p值

density.default()间接给你的是什么，除了512个等距点，而不是连续点。对于特定值x，您可以找到density.default()和相应y输出的最接近值：

# Create some random data:
set.seed(123)
averageRank <- rnorm(100, mean = 300, sd = 75)
x           <- 250

kde       <- density.default(averageRank)
closeness <- abs(kde$x - x)
kde$y[which(closeness == min(closeness))]
# [1] 0.004370022    # is P(X = 250) 

然而，从您提出问题的方式来看，我感觉您正在寻找累积分布函数：

sum(kde$y[1:which(closeness == min(closeness))]/sum(kde$y)
# [1] 0.2138626     # is P(X < 250)