Generalise scipy.stats.binned_statistic_dd to handle multi dimensional values -- 2

Currently, the module is general with regards to binning, but not data. For reference, see: [login to view URL] I would need a module, that build on binned_statistics_dd, so as it can handle "values" being a sequence of E arrays of length N, or as an (N,E) array. The purpose if this is that I want to be able to pass a function that uses more than one list as input, allowing me to solve linear least square or similar problems.