我想用qq图来检验一个分布是否为对数正态分布。因此,为了方便起见,我使用stats创建了一个对数正态分布,并使用dist函数在探测图中检查它。但是我仍然没有得到我在45度线上的分数,即qq图显示我的z不是对数正态的。为什么?我是新来的,任何帮助都将不胜感激。下面是我的代码
z=stats.lognorm.rvs(2.5,size=500)
stats.probplot(z,dist=stats.lognorm, sparams=(2.5,),plot=pylab)
pylab.show()这是我的结果:output
发布于 2020-01-11 19:47:05
这似乎是对数正态分布指数性质的统计伪影。注意图的最右边的点之间有很大的间隙。有趣的是,即使样本量很大,它似乎也不是很好的视觉拟合。
但是如果你首先取对数,你会发现随机采样点确实是正态分布的(在对数尺度上)。
from scipy import stats
import pylab
import numpy as np
z=np.log(stats.lognorm.rvs(2.5,size=500))
stats.probplot(z,dist=stats.norm, sparams=(0,2.5),plot=pylab)
pylab.show()因此,从编程的角度来看,您的SciPy代码似乎工作正常。出于数学原因,最好在https://stats.stackexchange.com/上询问
https://stackoverflow.com/questions/59692757
复制相似问题