text: Title: jointplot with kind="hex" fails with datetime64[ns]
  Body: Minimal example:
  ```python
  import seaborn as sns
  import numpy as np
  dates = np.array(["2023-01-01", "2023-01-02", "2023-01-03"], dtype="datetime64[ns]")
  sns.jointplot(x=dates, y=[1, 2, 3], kind="hex")
  ```
  Error:
  ```
  Traceback (most recent call last):
    File "/.../seaborn_bug.py", line 21, in <module>
      sns.jointplot(x=dates, y=[1, 2, 3], kind="hex")
    File "/Library/Frameworks/Python.framework/Versions/3.10/lib/python3.10/site-packages/seaborn/axisgrid.py", line 2307, in jointplot
      x_bins = min(_freedman_diaconis_bins(grid.x), 50)
    File "/Library/Frameworks/Python.framework/Versions/3.10/lib/python3.10/site-packages/seaborn/distributions.py", line 2381, in _freedman_diaconis_bins
      iqr = np.subtract.reduce(np.nanpercentile(a, [75, 25]))
  TypeError: the resolved dtypes are not compatible with subtract.reduce. Resolved (dtype('<M8[ns]'), dtype('<M8[ns]'), dtype('<m8[ns]'))
  ```
  I think this should work, as datetime64[ns] is the default type for datetimes in pandas. It works when I omit `kind="hex"` or use `kind="kde"`. The error was in version 0.13.2. In 0.11.2 I got a error in the same cases but it was a integer overflow with numpy during conversions.