Scipy.stats的fit和pdf函数

Question

我想知道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的含义以及它们的计算方式。

谢谢。