Scipy.stats的fit和pdf函数

时间:2018-09-07 13:12:35

标签: python scipy

我想知道scipy.stats如何使用其方法fitpdf。根据文档,fit(data, a, loc = 0, scale = 1)估计datapdf(x, a, loc=0, scale=1) computes probability density function的参数。但是我找不到统计和数学上fitpdf的实际执行情况。

我正在使用sm.datasets.elnino数据,正在使用tmthydvnprt的code

import warnings
import numpy as np
import pandas as pd
import scipy.stats as st
import statsmodels as sm
import matplotlib
import matplotlib.pyplot as plt

data = pd.Series(sm.datasets.elnino.load_pandas().data.set_index('YEAR').values.ravel())


y, x = np.histogram(data, bins = 50, density = True)
x = (x + np.roll(x, -1))[:-1] / 2.0
distribution = st.gennorm

params = distribution.fit(data)
arg = params[:-2]
loc = params[-2]
scale = params[-1]
pdf = distribution.pdf(x, loc = loc, scale = scale, *arg)
sse = np.sum(np.power(y - pdf, 2.0))

使用dataarg = 4.3836,loc = 23.2991,scale = 3.8499。

我想知道arglocscale的含义以及它们的计算方式。

谢谢。

0 个答案:

没有答案