feat: p-values added to skew and kurtosis in Obs.plot_history

This commit is contained in:
Fabian Joswig 2022-01-17 14:55:14 +00:00
parent 4f3b2d22d5
commit 0918330f13

View file

@ -4,7 +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
from scipy.stats import skew, skewtest, kurtosis, kurtosistest
import numdifftools as nd
from itertools import groupby
from .covobs import Covobs
@ -565,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 + f', skew: {skew(y):2.3f}, kurtosis: {kurtosis(y):2.3f}')
plt.title(e_name + f', skew: {skew(y):.3f} (p={skewtest(y).pvalue:.2f}), kurtosis: {kurtosis(y):.3f} (p={kurtosistest(y).pvalue:.2f})')
plt.draw()
def plot_piechart(self):