我想知道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))
使用data
,arg
=4.3836,loc
=23.2991,scale
=3.8499。在
我想知道arg
、loc
和{
谢谢。在
目前没有回答
相关问题 更多 >
编程相关推荐