Scipy.stats公司的fit和pdf函数

2024-04-20 04:37:04 发布

您现在位置:Python中文网/ 问答频道 /正文

我想知道scipy.stats如何使用其方法fit和{}。根据文档,fit(data, a, loc = 0, scale = 1)估计data的参数,而pdf(x, a, loc=0, scale=1) computes probability density function。但我无法找到fit和{}是如何在统计和数学上执行的。在

我使用的是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。在

我想知道argloc和{}代表什么以及它们是如何计算的。

谢谢。在


Tags: importdatapdfasstatsnpargscipy