COVAR

semver fix

Source/Statistics/COVAR_Series.cs

@@ -0,0 +1,67 @@
+namespace QuanTAlib;
+using System;
+
+/* <summary>
+COVAR: Covariance
+   Covariance is defined as the expected value (or mean) of the product
+   of their deviations from their individual expected values.
+
+Sources:
+  https://en.wikipedia.org/wiki/Covariance
+
+</summary> */
+
+public class COVAR_Series : Pair_TSeries_Indicator
+{
+  public COVAR_Series(TSeries d1, TSeries d2, int period, bool useNaN = false) : base(d1, d2, period, useNaN)
+  {
+    if (base._d1.Count > 0 && base._d2.Count > 0) { for (int i = 0; i < base._d1.Count; i++) { this.Add(base._d1[i], base._d2[i], false); } }
+  }
+
+  private readonly System.Collections.Generic.List<double> _x = new();
+  private readonly System.Collections.Generic.List<double> _xx = new();
+  private readonly System.Collections.Generic.List<double> _y = new();
+  private readonly System.Collections.Generic.List<double> _yy = new();
+  private readonly System.Collections.Generic.List<double> _xy = new();
+
+  public override void Add((System.DateTime t, double v) TValue1, (System.DateTime t, double v) TValue2, bool update)
+  {
+    if (update)
+    {
+      _x[_x.Count - 1] = TValue1.v;
+      _xx[_xx.Count - 1] = TValue1.v * TValue1.v;
+      _y[_y.Count - 1] = TValue2.v;
+      _y[_yy.Count - 1] = TValue2.v * TValue2.v;
+      _xy[_xy.Count - 1] = TValue1.v * TValue2.v;
+    }
+    else
+    {
+      _x.Add(TValue1.v);
+      _xx.Add(TValue1.v * TValue1.v);
+      _y.Add(TValue2.v);
+      _yy.Add(TValue2.v * TValue2.v);
+      _xy.Add(TValue1.v * TValue2.v);
+    }
+    if (_x.Count > this._p) { _x.RemoveAt(0); }
+    if (_xx.Count > this._p) { _xx.RemoveAt(0); }
+    if (_y.Count > this._p) { _y.RemoveAt(0); }
+    if (_yy.Count > this._p) { _yy.RemoveAt(0); }
+    if (_xy.Count > this._p) { _xy.RemoveAt(0); }
+
+    double _sumx = 0;
+    for (int i = 0; i < _x.Count; i++) { _sumx += _x[i]; }
+    double _sumxx = 0;
+    for (int i = 0; i < _xx.Count; i++) { _sumxx += _xx[i]; }
+    double _sumy = 0;
+    for (int i = 0; i < _y.Count; i++) { _sumy += _y[i]; }
+    double _sumyy = 0;
+    for (int i = 0; i < _yy.Count; i++) { _sumyy += _yy[i]; }
+    double _sumxy = 0;
+    for (int i = 0; i < _xy.Count; i++) { _sumxy += _xy[i]; }
+
+    double _covar = (_sumxy / _p) - ((_sumx / _p) * (_sumy / _p));
+
+    var result = (TValue1.t, (this.Count < this._p - 1 && this._NaN) ? double.NaN : _covar);
+    if (update) { base[base.Count - 1] = result; } else { base.Add(result); }
+  }
+}