我有一个Dask Dataframe,我想为列表列表计算偏度,如果这个偏差超过某个阈值,我会使用日志转换来纠正它。我想知道是否有更有效的方法通过删除下面的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'])
在此示例中,您可以采取一些措施来改善并行性:
您可以使用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逻辑并适当地应用da.log
。这仍然需要两次通过您的数据,但这应该是一个很好的改进,而不是你现在拥有的。