alternative scoring function so that it can work with scipy optimization functions
changed rates
changed rates