Try reversing the order of the datasets for a and b. I believe the type of dataset will be determined from the 'a' dataset, i.e. functional bucket in this case. That type allows the statistical information. The new expression would have to be switched to "a*step(b)". It should be possible to do the 3drefit too, but that would require changing the type first.