如何在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如何实现它也可以与负数一起使用?

解决方法

暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!

如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@)

相关问答

依赖报错 idea导入项目后依赖报错,解决方案:https://blog....
错误1:代码生成器依赖和mybatis依赖冲突 启动项目时报错如下...
错误1:gradle项目控制台输出为乱码 # 解决方案:https://bl...
错误还原:在查询的过程中,传入的workType为0时,该条件不起...
报错如下,gcc版本太低 ^ server.c:5346:31: 错误:‘struct...