Update jointplot/pairplot to use new distribution functions

[mwaskom]

c.f. #2157

[mwaskom]

Forgot that this was mostly taken care of in #2125, but there is an issue where code will crash when rug= is defined in the marginal_kws dict. (Maybe other distplot-specific kws too?) This appeared in an API example.