我有一个Dask数据框,我想为其计算列列表的偏斜度,如果此偏斜度超出某个阈值,则可以使用对数转换对其进行校正。我想知道是否有一种更有效的方法,通过删除以下correct_skewness()
函数中的 for循环,使correct_skewness()
函数在多列上并行工作:
import dask
import dask.array as da
from scipy import stats
# Create a dataframe
df = dask.datasets.timeseries()
df.head()
id name x y
timestamp
2000-01-01 00:00:00 1032 Oliver 0.018604 0.089191
2000-01-01 00:00:01 1032 Norbert 0.666689 -0.979374
2000-01-01 00:00:02 991 Victor 0.027691 -0.474660
2000-01-01 00:00:03 979 Kevin 0.320067 0.656949
2000-01-01 00:00:04 1087 Zelda -0.462076 0.513409
def correct_skewness(columns=None, max_skewness=2):
if columns is None:
raise ValueError(
f"columns argument is None. Please set columns argument to a list of columns"
)
for col in columns:
skewness = stats.skew(df[col])
max_val = df[col].max().compute()
min_val = df[col].min().compute()
if abs(skewness) > max_skewness and (max_val > 1 or min_val < 0):
delta = 1.0
if min_val < 0:
delta = max(1, -min_val + 1)
df[col] = da.log(delta + df[col])
return df
df = correct_skewness(columns=['x', 'y'])
答案 0 :(得分:2)
在此示例中,您可以做一些事情来改善并行度:
您可以使用dask.array.stats.skew而不是statsmodels.skew。您将必须import dask.array.stats
明确
您可以一次计算所有列的最小值/最大值
mins = [df[col].min() for col in cols]
maxes = [df[col].min() for col in cols]
skews = [da.stats.skew(df[col]) for col in cols]
mins, maxes, skews = dask.compute(mins, maxes, skews)
然后,您可以执行if-logic并适当地应用da.log
。仍然需要对数据进行两次传递,但这应该是对现有数据的一个很好的改进。