Recommended way to compute difference statistics

[phinate]

Hi! When using the create dataset functionality, I'm interested in the most efficient way to create statistics based on the per-variable timestep differences.

Should I be doing some kind of chunk-wise call to np.diff and/or the mean/variance of those values? Or is there a way to create a dataset of differences, and then use the built-in writing of statistics to get these natively?