Expose ddof on covariance

lib/explorer/backend/lazy_series.ex

@@ -110,7 +110,7 @@ defmodule Explorer.Backend.LazySeries do
     nil_count: 1,
     skew: 2,
     correlation: 3,
-    covariance: 2,
+    covariance: 3,
     # Strings
     contains: 2,
     replace: 3,
@@ -504,8 +504,8 @@ defmodule Explorer.Backend.LazySeries do
   end
 
   @impl true
-  def covariance(%Series{} = left, %Series{} = right) do
-    args = [series_or_lazy_series!(left), series_or_lazy_series!(right)]
+  def covariance(%Series{} = left, %Series{} = right, ddof \\ 1) do
+    args = [series_or_lazy_series!(left), series_or_lazy_series!(right), ddof]
     data = new(:covariance, args, {:f, 64}, true)
 
     Backend.Series.new(data, {:f, 64})

lib/explorer/backend/series.ex

@@ -91,7 +91,7 @@ defmodule Explorer.Backend.Series do
   @callback skew(s, bias? :: boolean()) :: float() | non_finite() | lazy_s() | nil
   @callback correlation(s, s, ddof :: non_neg_integer()) ::
               float() | non_finite() | lazy_s() | nil
-  @callback covariance(s, s) :: float() | non_finite() | lazy_s() | nil
+  @callback covariance(s, s, ddof :: non_neg_integer()) :: float() | non_finite() | lazy_s() | nil
 
   # Cumulative
 

lib/explorer/polars_backend/expression.ex

@@ -76,7 +76,7 @@ defmodule Explorer.PolarsBackend.Expression do
     unordered_distinct: 1,
     variance: 2,
     skew: 2,
-    covariance: 2
+    covariance: 3
   ]
 
   @first_only_expressions [

lib/explorer/polars_backend/native.ex

@@ -262,7 +262,7 @@ defmodule Explorer.PolarsBackend.Native do
   def s_cumulative_product(_s, _reverse), do: err()
   def s_skew(_s, _bias), do: err()
   def s_correlation(_s1, _s2, _ddof), do: err()
-  def s_covariance(_s1, _s2), do: err()
+  def s_covariance(_s1, _s2, _ddof), do: err()
   def s_distinct(_s), do: err()
   def s_divide(_s, _other), do: err()
   def s_dtype(_s), do: err()

lib/explorer/polars_backend/series.ex

@@ -224,8 +224,8 @@ defmodule Explorer.PolarsBackend.Series do
     do: Shared.apply_series(matching_size!(left, right), :s_correlation, [right.data, ddof])
 
   @impl true
-  def covariance(left, right),
-    do: Shared.apply_series(matching_size!(left, right), :s_covariance, [right.data])
+  def covariance(left, right, ddof),
+    do: Shared.apply_series(matching_size!(left, right), :s_covariance, [right.data, ddof])
 
   # Cumulative
 

lib/explorer/series.ex

@@ -2518,10 +2518,14 @@ defmodule Explorer.Series do
       3.0
   """
   @doc type: :aggregation
-  @spec covariance(left :: Series.t() | number(), right :: Series.t() | number()) ::
+  @spec covariance(
+          left :: Series.t() | number(),
+          right :: Series.t() | number(),
+          ddof :: non_neg_integer()
+        ) ::
           float() | non_finite() | nil
-  def covariance(left, right) do
-    basic_numeric_operation(:covariance, left, right)
+  def covariance(left, right, ddof \\ 1) do
+    basic_numeric_operation(:covariance, left, right, [ddof])
   end
 
   # Cumulative

native/explorer/src/expressions.rs

@@ -501,11 +501,9 @@ pub fn expr_correlation(left: ExExpr, right: ExExpr, ddof: u8) -> ExExpr {
 }
 
 #[rustler::nif]
-pub fn expr_covariance(left: ExExpr, right: ExExpr) -> ExExpr {
+pub fn expr_covariance(left: ExExpr, right: ExExpr, ddof: u8) -> ExExpr {
     let left_expr = left.clone_inner().cast(DataType::Float64);
     let right_expr = right.clone_inner().cast(DataType::Float64);
-    // TODO: make this a parameter.
-    let ddof: u8 = 1;
     ExExpr::new(cov(left_expr, right_expr, ddof))
 }
 

native/explorer/src/series.rs

@@ -1006,11 +1006,9 @@ pub fn s_correlation(
 }
 
 #[rustler::nif(schedule = "DirtyCpu")]
-pub fn s_covariance(env: Env, s1: ExSeries, s2: ExSeries) -> Result<Term, ExplorerError> {
+pub fn s_covariance(env: Env, s1: ExSeries, s2: ExSeries, ddof: u8) -> Result<Term, ExplorerError> {
     let s1 = s1.clone_inner().cast(&DataType::Float64)?;
     let s2 = s2.clone_inner().cast(&DataType::Float64)?;
-    // TODO: make ddof a parameter.
-    let ddof: u8 = 1;
     let cov = cov(s1.f64()?, s2.f64()?, ddof);
     Ok(term_from_optional_float(cov, env))
 }