BUG, TST, DOC: fix core.missing._akima_interpolate will raise AttributeError

pandas/core/missing.py

@@ -445,8 +445,9 @@ def _akima_interpolate(xi, yi, x, der=0, axis=0):
         A 1-D array of real values.  `yi`'s length along the interpolation
         axis must be equal to the length of `xi`. If N-D array, use axis
         parameter to select correct axis.
-    x : scalar or array_like of length M.
-    der : int or list, optional
+    x : scalar or array_like
+        Of length M.
+    der : int, optional
         How many derivatives to extract; None for all potentially
         nonzero derivatives (that is a number equal to the number
         of points), or a list of derivatives to extract. This number
@@ -468,12 +469,7 @@ def _akima_interpolate(xi, yi, x, der=0, axis=0):
 
     P = interpolate.Akima1DInterpolator(xi, yi, axis=axis)
 
-    if der == 0:
-        return P(x)
-    elif interpolate._isscalar(der):
-        return P(x, der=der)
-    else:
-        return [P(x, nu) for nu in der]
+    return P(x, nu=der)
 
 
 def _cubicspline_interpolate(xi, yi, x, axis=0, bc_type="not-a-knot", extrapolate=None):

pandas/tests/series/methods/test_interpolate.py

@@ -133,17 +133,28 @@ def test_interpolate_akima(self):
 
         ser = Series([10, 11, 12, 13])
 
+        # interpolate at new_index where `der` is zero
         expected = Series(
             [11.00, 11.25, 11.50, 11.75, 12.00, 12.25, 12.50, 12.75, 13.00],
             index=Index([1.0, 1.25, 1.5, 1.75, 2.0, 2.25, 2.5, 2.75, 3.0]),
         )
-        # interpolate at new_index
         new_index = ser.index.union(Index([1.25, 1.5, 1.75, 2.25, 2.5, 2.75])).astype(
             float
         )
         interp_s = ser.reindex(new_index).interpolate(method="akima")
         tm.assert_series_equal(interp_s[1:3], expected)
 
+        # interpolate at new_index where `der` is a non-zero int
+        expected = Series(
+            [11.0, 1.0, 1.0, 1.0, 12.0, 1.0, 1.0, 1.0, 13.0],
+            index=Index([1.0, 1.25, 1.5, 1.75, 2.0, 2.25, 2.5, 2.75, 3.0]),
+        )
+        new_index = ser.index.union(Index([1.25, 1.5, 1.75, 2.25, 2.5, 2.75])).astype(
+            float
+        )
+        interp_s = ser.reindex(new_index).interpolate(method="akima", der=1)
+        tm.assert_series_equal(interp_s[1:3], expected)
+
     @td.skip_if_no_scipy
     def test_interpolate_piecewise_polynomial(self):
         ser = Series([10, 11, 12, 13])