我想知道scipy.stats
如何使用其方法fit
和pdf
。根据文档,fit(data, a, loc = 0, scale = 1)
估计data
时pdf(x, a, loc=0, scale=1) computes probability density function
的参数。但是我找不到统计和数学上fit
和pdf
的实际执行情况。
我正在使用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))
使用data
,arg
= 4.3836,loc
= 23.2991,scale
= 3.8499。
我想知道arg
,loc
和scale
的含义以及它们的计算方式。
谢谢。