将直方图拟合到自由度未知的 chi2

问题描述

有人可以验证某些数据的 scipy.stats.chi2.fit 中的第一个参数实际上是自由度数 df 吗?这是我正在检查的示例

x=np.random.rand(1000)
y=stats.chi2.pdf(x,df=1)
z=np.log10(y)
# assume z is the distribution of a rv in logspace

# is this line correct for the mles?
df,loc,scale = stats.chi2.fit(z)

仅当我检查例如类似的东西

x=np.random.rand(1000)
y=stats.chi2.pdf(x,df=1.5)
z=np.log10(y)
# assume z is the distribution in logspace
df,scale = stats.chi2.fit(z)

它回来了

(2.686925741184449,-0.5319434585863583,0.08411242959515701)

我不确定合身是否正确?

谢谢!

解决方法

暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!

如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@)