The current statistical_inefficiency function only returns the subsampled data, but I also want the value of statistical inefficiency. And I have to use the statistical_inefficiency in 'pymbar', and do a similar calculation again.

It'd help if statistical_inefficiency returns the subsampled data as well as the value of statistical inefficiency or has an option to control whether to return the value.