From 4f3b2d22d5d74c50d4faeac6f830cf6b0b4e8679 Mon Sep 17 00:00:00 2001
From: Fabian Joswig <fabian.joswig@ed.ac.uk>
Date: Mon, 17 Jan 2022 14:45:54 +0000
Subject: [PATCH] feat: skew and kurtosis added to Obs.plot_history

---
 pyerrors/obs.py | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/pyerrors/obs.py b/pyerrors/obs.py
index e19cc617..546a0a2e 100644
--- a/pyerrors/obs.py
+++ b/pyerrors/obs.py
@@ -4,6 +4,7 @@ import numpy as np
 import autograd.numpy as anp  # Thinly-wrapped numpy
 from autograd import jacobian
 import matplotlib.pyplot as plt
+from scipy.stats import kurtosis, skew
 import numdifftools as nd
 from itertools import groupby
 from .covobs import Covobs
@@ -564,7 +565,7 @@ class Obs:
             y = np.concatenate(tmp, axis=0)
             plt.errorbar(x, y, fmt='.', markersize=3)
             plt.xlim(-0.5, e_N - 0.5)
-            plt.title(e_name)
+            plt.title(e_name + f', skew: {skew(y):2.3f}, kurtosis: {kurtosis(y):2.3f}')
             plt.draw()
 
     def plot_piechart(self):