如何在 Python 中使用零偏度对数变换

如何在Python中进行零偏度对数变换？

例如在 Stata 中它的实现是lnskew0 (see https://www.stata.com/manuals13/rlnskew0.pdf https://www.stata.com/manuals13/rlnskew0.pdf).

我没有找到Python 的实现。有人知道实施吗？

否则，第一次尝试将是：

from scipy.stats import skew
import numpy as np
from scipy.optimize import root_scalar


def lnskew0(x):
    def skew_ln(k):
        return skew(np.log(x - k))
    res = root_scalar(
        skew_ln,
        bracket=[-x.min(), x.max()*0.99999],
        method='bisect'
    )
    return np.log(x - res.root)

运行良好numpy仅包含正数的数组。 Stata怎么样lnskew0是否实现了它也适用于负数？

我又试了一次，使它也适用于负数：

from scipy.stats import skew
import numpy as np
from scipy.optimize import root_scalar


def lnskew0(x):

    x0 = x + 1

    def skew_ln_pos(k):
        with np.errstate(all='ignore'):
            return skew(np.log(x0 - k))
    res_pos = root_scalar(
        skew_ln_pos,
        bracket=[-150, 150],
        method='bisect'
    )

    def skew_ln_neg(k):
        with np.errstate(all='ignore'):
            return skew(np.log(-x0 - k))
    res_neg = root_scalar(
        skew_ln_neg,
        bracket=[-150, 150],
        method='bisect'
    )

    res = (res_pos.root - 1, res_neg.root + 1)

    lnskew0_res = (
        np.log(x - res[0]),
        np.log(-x - res[1])
    )

    whichmin = np.nanargmin([abs(skew(x)) for x in lnskew0_res])
    return lnskew0_res[whichmin]

注意：它仍然有一个问题。的括号root_scalar需要手动选择。