In [1]: import scipy.stats as st
...: arr = st.norm.rvs(loc=6000, scale=20, size=(300, 10))
...: st.skewtest(arr)
Out[1]:
SkewtestResult(statistic=array([ 1.21508957, -0.60569814, 0.07676428, 0.10929194, -1.78128941,
1.3046305 , 1.03649957, -0.90213437, 1.39273949, 0.87478829]), pvalue=array([0.22433193, 0.54471522, 0.93881107, 0.91297094, 0.07486518,
0.1920187 , 0.29996914, 0.36698549, 0.16369858, 0.38168911]))
In [2]: import pnumpy; pnumpy.enable()
In [3]: st.skewtest(arr)
Out[3]:
SkewtestResult(statistic=array([-4.34541497e+03, -4.34723576e+03, -4.34655330e+03, -4.34652077e+03,
-4.34841135e+03, -4.34532543e+03, -4.34559356e+03, -4.34753220e+03,
1.39273949e+00, 8.74788294e-01]), pvalue=array([0. , 0. , 0. , 0. , 0. ,
0. , 0. , 0. , 0.16369858, 0.38168911]))
In [4]: pnumpy.__version__
Out[4]: '2.0.23'
