我是R的新手,正在研究qnorm(),我注意到qnorm()的值与Summary()中显示的值不同。
代码:
set.seed(123)
X<-rnorm(100,10,1)
summary(X)
qnorm(0.5,10,1)
qnorm(0.25,10,1)
qnorm(0.75,10,1)
> summary(X)
Min. 1st Qu. Median Mean 3rd Qu. Max.
7.691 9.506 10.062 10.090 10.692 12.187
> qnorm(0.5,10,1)
[1] 10
> qnorm(0.25,10,1)
[1] 9.32551
> qnorm(0.75,10,1)
[1] 10.67449
qnorm()
为0.25时是否应与第一分位数相同,因为qnorm()
返回特定cdf的值,该值也可用于查找特定分位数?
1条答案
按热度按时间2wnc66cl1#
简短的回答是:不,但随着n的增加,它会向相同的方向收敛。
考虑
x <- rnorm(1, 10, 1)
,您是否希望此单个值同时具有qnorm(0.25, 10, 1)
、qnorm(0.5, 10, 1)
和qnorm(0.75, 10, 1)
的值?在模拟值时,您是从具有某些属性的分布中随机抽取批次。在本例中,它是均值为10、标准差为1的正态分布。样本分布将向理论分布收敛。但对于100个样本,你应该预料到观察到的样本分布在平均值,标准差方面略有偏离理论分布,但样本分位数更是如此。尝试将
n
增加到一个更大的值(比如1e7),您将看到样本分位数收敛。