如何在Python中使用零偏度日志转换

如何解决如何在Python中使用零偏度日志转换

如何在Python中进行零偏度日志转换？

例如，在Stata中，它是在lnskew0中实现的（请参见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如何实现它也可以与负数一起使用？