jarque bera test.py

from scipy.stats import jarque_bera
import numpy as np

# generate a sample dataset
data = np.random.normal(0, 1, 1000)

# calculate the test statistic and p-value
jb_stat, jb_p = jarque_bera(data)

# print the results
print("Jarque-Bera statistic:", jb_stat)
print("p-value:", jb_p)

# interpret the results
if jb_p < 0.05:
    print("The dataset is not normally distributed.")
else:
    print("The dataset is normally distributed.")