Added function to compute average ru/mi over threshold values and plot them: RUMIplot