使用 scipy.stats 将数据与自定义分布拟合答案

【问题标题】：Fitting data with a custom distribution using scipy.stats使用 scipy.stats 将数据与自定义分布拟合
【发布时间】：2018-09-07 02:53:08
【问题描述】：

所以我注意到scipy 中没有实现Skewed generalized t distribution。将其分布到我拥有的某些数据中对我很有用。不幸的是，fit 在这种情况下似乎对我不起作用。为了进一步解释，我已经这样实现了它

import numpy as np
import pandas as pd
import scipy.stats as st
from scipy.special import beta

class sgt(st.rv_continuous):

    def _pdf(self, x, mu, sigma, lam, p, q):

        v = q ** (-1 / p) * \
            ((3 * lam ** 2 + 1) * (
                    beta(3 / p, q - 2 / p) / beta(1 / p, q)) - 4 * lam ** 2 *
             (beta(2 / p, q - 1 / p) / beta(1 / p, q)) ** 2) ** (-1 / 2)

        m = 2 * v * sigma * lam * q ** (1 / p) * beta(2 / p, q - 1 / p) / beta(
            1 / p, q)

        fx = p / (2 * v * sigma * q ** (1 / p) * beta(1 / p, q) * (
                abs(x - mu + m) ** p / (q * (v * sigma) ** p) * (
                lam * np.sign(x - mu + m) + 1) ** p + 1) ** (
                          1 / p + q))

        return fx

    def _argcheck(self, mu, sigma, lam, p, q):

        s = sigma > 0
        l = -1 < lam < 1
        p_bool = p > 0
        q_bool = q > 0

        all_bool = s & l & p_bool & q_bool

        return all_bool

这一切都很好，我可以生成具有给定参数的随机变量没有问题。 _argcheck 是必需的，因为只检查简单的正参数是不合适的。

sgt_inst = sgt(name='sgt')
vars = sgt_inst.rvs(mu=1, sigma=3, lam = -0.1, p = 2, q = 50, size = 100)

但是，当我尝试fit 这些参数时，我得到一个错误

sgt_inst.fit(vars)

RuntimeWarning：减法中遇到无效值
numpy.max(numpy.abs(fsim[0] - fsim[1:]))

它只是返回

我觉得奇怪的是，当我实现docs中所示的示例自定义高斯分布时，运行fit方法没有问题。

有什么想法吗？

【问题讨论】：

标签： python numpy scipy statistics

【解决方案1】：

正如fit docstring 所说，

拟合的初始估计由输入参数给出；对于未提供初始估计值的任何参数，调用self._fitstart(data) 来生成此类参数。

调用sgt_inst._fitstart(data)返回(1.0, 1.0, 1.0, 1.0, 1.0, 0, 1)（前五个是形状参数，后两个是loc和scale）。看起来_fitstart 不是一个复杂的过程。它选择的参数 l 不符合您的 argcheck 要求。

结论：为fit提供你自己的起始参数，例如，

sgt_inst.fit(data, 0.5, 0.5, -0.5, 2, 10)

为我的随机数据返回 (1.4587093459289049, 5.471769032259468, -0.02391466905874927, 7.07289326147152 4, 0.741434497805832, -0.07012808188413872, 0.5308181287869771)。

【讨论】：