Predicted rank abundance should use beta from predicted values

The current approach is actually fitting the model using observed
S. This approach fits the model using the predicted values and
then switches to using observed S to make predictions for
comparable numbers of "species", i.e., points along the CDF.