如何检查值是否偶然出现在R中的分布中

Question

我有以下随机生成的分布:

set.seed(1)
mean=100; sd=15
x <- seq(-4,4,length=100)*sd + mean
hx <- dnorm(x,mean,sd)

plot(x, hx, type="l", lty=2, xlab="x value",
     ylab="Density", main="Some random distribution")

并且是"非随机"的价值

set.seed(1)
x <- seq(-4,4,length=100)*10 + 100
ux <- dunif(x = x, min=10, max=100)
non_random_value <- ux[1]
non_random_value
# [1] 0.01111111

我希望显示的统计数据non_random_value很重要而且不会偶然出现hx.我怎么能在R？

Answer 1

你想要的功能是pnorm(x,mean,sd).它返回正态分布中由mean和sd定义的小于x的值的比例.

您可以在双尾测试中使用pnorm来获得比x更极端的值的比例(即,从任一方向的均值更远).

p <- 2*pnorm(x, mean, sd, lower.tail=x<mean)

现在将p解释为平均p%的潜在值比x更远离均值.请记住,x的任何值都可以从正态分布中偶然出现.你不能说x不是随机的,只是看起来不太可能来自指定的分布.

这是一个很好的网站,描述了这个更详细.