Support factor predictors with continuous moderators in sim_slopes/johnson_neyman

[jacob-long]

This would involve finding a j_n interval at each level of the factor predictor. For now, would involve too much of a rework of the underlying functions for me to do it. It's not too much unlike what is done for second moderators, though.