The range 0 < r < 1 is mapped on one side and the range 1 < r < Infinity is mapped on the other side of the y axis. You have to use a transformation to get something like a similar distribution of curves on both sides. I used the following, based on the x-intercept p being related to the ratio r by r = -(p-5)/(p + 5) = (5-p)/(5+p) if 0 < p < 5. The code for this adjustment is

{r,Union[
-((Range[9]/2-5)/(Range[9]/2+5)),
{1},
-((Range[9]/2+5)/(Range[9]/2-5))]}

I made some other small adjustments for the sake of performance, which you should feel no need to accept if they interfere with your purpose.