生成一个样本
import numpy as np
from scipy import stats
'''
在这里,我们将t分布的所需形状参数(在统计中对应于自由度)设置为10。使用size = 1000表示我们的样本包含1000个独立绘制的(伪)随机数。由于未指定关键字参数loc和scale,因此将它们设置为默认值零和一。
'''
np.random.seed(2020)
x=stats.t.rvs(10,size=1000)
d:\Anaconda3\lib\importlib\_bootstrap.py:219: RuntimeWarning: numpy.ufunc size changed, may indicate binary incompatibility. Expected 192 from C header, got 216 from PyObject
return f(*args, **kwds)
d:\Anaconda3\lib\importlib\_bootstrap.py:219: RuntimeWarning: numpy.ufunc size changed, may indicate binary incompatibility. Expected 192 from C header, got 216 from PyObject
return f(*args, **kwds)
x
array([-1.79937219e+00, -9.16947599e-01, -6.23672658e-02, 6.93167535e-01,
1.80118163e+00, -1.02592060e+00, 4.01066025e-01, -1.91015338e-01,
-1.19563987e-01, -2.26134986e-01, 3.23050886e-01, -9.22061883e-01,
1.35064128e+00, 5.42949841e-01, -4.60522342e-01, 1.01424367e+00,
-2.20308350e-01, -2.59316115e-01, 1.52749609e-01, 4.63332848e-01,
-7.63435123e-01, 3.51077958e-01, -1.20694886e+00, -9.74257986e-02,
-5.82930066e-02, 2.90404773e-01, -1.61086154e-01, -1.20916124e-01,
7.52397419e-01, 3.99534824e-01, 8.25653983e-01, 6.69627593e-01,
-2.08560499e+00, -2.33150645e+00, 1.26080788e+00, -1.06018145e+00,
-2.64190038e-01, -6.06938922e-01, -2.76862279e-01, 2.19626508e+00,
-5.08492116e-01, 1.03688480e+00, 7.64800971e-01, 6.85051528e-01,
1.44908752e-01, 7.84647158e-02, -1.33507204e+00, -2.91500267e-03,
-3.41229510e-01, -1.71745063e+00, -1.05830288e+00, 3.64122554e-01,
-3.08349226e-01, -5.16484855e-01, -8.98755695e-01, 1.19968196e-01,
7.75893767e-01, 7.48960897e-01, -5.22244842e-01, -5.01603404e-01,
2.42525394e-01, 1.11303062e+00, -2.40953584e-01, 1.08495827e-01,
-3.48788202e-01, 1.90611916e-01, 4.64715652e-01, 1.44927602e+00,
4.62647827e-01, 1.29131093e+00, -1.91232214e+00, 7.63557003e-01,
4.39071347e-01, 2.70055877e-01, -4.90372934e-01, 8.64151858e-01,
3.86246914e-01, 1.36265801e+00, -1.88792648e+00, 1.36666145e+00,
-5.16043304e-01, -2.22099152e+00, -7.18979232e-01, -6.20482628e-01,
-7.15841369e-01, -4.54799087e-01, -9.55602901e-02, -3.45630716e-01,
1.24773888e+00, 7.81183082e-01, -3.70699115e-01, 2.83458583e+00,
4.89026489e-02, -4.22891142e-01, -5.58655211e-01, -4.02121958e-01,
-9.72855346e-01, -2.64769519e-01, 3.06353094e-01, 7.88867173e-01,
3.34265051e-02, 1.51818327e-01, -1.01860792e+00, -2.00735348e+00,
1.53877595e+00, -1.83679404e-01, 2.62407339e-01, -2.36135884e+00,
1.35516891e-01, -1.22892795e+00, -4.58815730e-01, -1.43076291e+00,
-1.02552654e+00, -4.14160111e-01, 8.66666545e-01, 7.60102259e-01,
5.07226554e-02, 8.97690959e-01, -6.21611401e-01, -1.37674039e+00,
3.14506644e+00, -1.05763693e-01, -2.17657105e+00, 6.60361160e-01,
-9.98428328e-01, -1.09877083e+00, 1.18607328e+00, 1.40722016e+00,
2.57477060e+00, 1.31960769e-01, -3.69164901e-01, -9.64919664e-01,
-8.09420250e-01, 2.21967948e-01, -3.02601823e-01, -2.43457161e+00,
-1.75800766e+00, 1.94960559e-01, -4.70040682e-01, -1.54720971e+00,
2.83967200e-01, -1.31312282e+00, -1.21888345e+00, 2.17262865e+00,
-4.67741401e-01, 3.65184673e-01, -5.17542541e-01, -1.38431605e-01,
1.20425425e+00, -3.59548580e-01, 6.81515393e-01, 8.10493589e-01,
-1.90650953e+00, 3.47785640e-01, -7.93668736e-01, -2.12985117e+00,
3.28987538e+00, -1.92844296e+00, -2.93629540e-01, 2.89118286e-01,
-5.59469458e-01, -9.53687384e-01, -1.49961443e+00, 4.18682537e-01,
1.62982834e-01, -7.95447423e-01, -1.67759534e-01, 1.29297515e+00,
9.91589670e-01, -2.91405755e-01, 3.39254829e-01, -7.45802187e-01,
1.28416736e-01, -1.30075076e-02, 6.15615192e-01, -9.60827114e-01,
5.85180189e-01, -6.54609677e-03, -1.70012744e+00, -2.01086610e+00,
-1.55337488e+00, -1.48792330e-01, -4.87362644e-01, -1.46177538e-02,
-7.61429996e-01, 7.26066729e-03, 7.61362964e-01, -1.52205730e-01,
-1.21967690e+00, 1.80757986e-01, 4.83092697e-01, -5.73581980e-01,
-2.75049063e+00, -4.73260606e-01, 1.43346909e+00, 3.40165863e-01,
9.62642309e-01, 1.26195277e+00, -1.01484047e+00, 5.74035611e-01,
3.06959938e+00, -1.06769702e+00, -1.11481136e+00, -6.54523649e-01,
-3.64085138e-02, -1.60232903e+00, -4.73992998e-01, -2.53760021e-01,
-1.66595060e+00, 2.57650511e+00, -1.62012373e+00, 1.69824712e+00,
9.91285915e-01, -1.20208981e+00, -1.17894856e+00, 2.46615552e-01,
7.14680934e-01, 5.60965530e-01, 2.89920106e-01, -4.17970148e-02,
1.07576669e+00, 7.97204795e-01, 1.46495756e+00, -8.74192281e-01,
-4.13581633e-01, 7.52948119e-01, 6.42370467e-01, 8.37751364e-01,
3.49279637e-01, 3.20054985e-01, -1.22765020e+00, 1.15293577e+00,
1.19246800e+00, -1.45251058e+00, -8.57420389e-01, 7.01822972e-01,
-6.31196323e-01, 4.74819915e-01, -3.97675576e-01, 3.99861943e-01,
-3.22211572e-01, 4.69449430e-01, 1.90734110e-01, -5.77826928e-01,
1.39639925e-01, -6.77690141e-01, -3.24901738e-01, 7.26996665e-01,
-9.65225515e-01, -5.87541183e-01, -2.45436379e+00, -1.46933565e+00,
-3.50510148e-01, 6.21482343e-01, 1.07609277e-01, 7.10654443e-02,
-7.38847949e-01, 1.23334104e+00, 2.25958988e-01, -6.59597265e-01,
-3.44780451e-01, 3.94303881e-01, 2.89530992e-01, 1.78517003e+00,
-6.02216889e-01, -1.36605026e+00, 4.65872583e-02, 1.08180913e+00,
1.23751595e+00, -5.36831313e-01, -3.58100319e-01, 5.01469481e-01,
-1.32375572e+00, 9.04311818e-02, -1.02552670e+00, 1.33671711e+00,
-4.24495154e-01, 1.25094432e-01, -5.84271712e-01, -8.98355284e-01,
-6.31428046e-01, -1.28628534e+00, -1.38725639e-01, 6.39099571e-02,
-6.27458927e-02, 3.49503912e-01, -1.49697063e+00, 1.46353500e+00,
-1.12010503e-01, -8.32609087e-01, -3.93624330e-01, -7.95758265e-01,
-2.56427774e-01, 6.88291393e-01, 4.07696874e-01, 2.20811481e-01,
-1.77355407e+00, 1.53306673e+00, -8.77953903e-01, 1.75035291e-01,
1.41491266e-04, 5.14717682e-01, 1.07605822e+00, -7.72818413e-02,
-8.68092400e-01, 5.35060965e-01, 4.00531067e-01, -6.38972284e-01,
1.35654431e-01, -6.85395862e-01, -5.64921163e-01, -4.57333264e-01,
1.16398213e+00, 1.61745391e+00, -9.35842375e-01, -7.67413907e-01,
3.81203063e-01, 2.83199925e-01, -6.58560058e-01, -5.75123773e-01,
1.77037607e+00, 1.15374591e+00, 1.65424362e+00, -2.91742701e-01,
-8.67794157e-01, -1.13712304e+00, 2.22675326e+00, -7.18912178e-01,
-1.99020525e+00, -6.78266021e-01, 8.18983944e-01, -5.63827377e-01,
1.65298380e-01, 1.61059028e+00, -4.75257953e-01, 9.28472771e-01,
-2.14187239e+00, -2.09125662e-01, 1.83804388e-01, -1.10644847e+00,
3.81194542e-01, -1.16031785e+00, -7.93882599e-01, 1.39007476e+00,
-1.05070145e-01, -3.58467504e-01, -7.57766660e-02, 1.66509589e+00,
-7.64445677e-01, -7.78358147e-01, -3.76169489e-01, 2.29208116e+00,
-3.55574465e+00, 1.93148247e-01, 3.69624880e-01, 1.04401457e+00,
-4.82090136e-01, 3.02605257e-01, 1.27945570e+00, 1.42889518e+00,
-1.25080706e-01, 1.15118961e-01, 1.06259696e+00, -4.93663778e-01,
-4.82725126e-01, -9.33504356e-01, 1.06817947e+00, 1.35031138e-01,
-4.72107707e-01, -9.01418963e-01, 1.02287328e+00, -1.09974515e+00,
-6.02084789e-01, -1.62122913e+00, 1.12320555e+00, -3.10711556e+00,
3.40643711e-01, -1.31555565e+00, 1.08834375e+00, 2.94294385e-01,
-1.04399975e+00, 6.56700165e-02, -4.04297787e-01, -2.50936648e+00,
-3.40409573e+00, -1.17415783e+00, 1.55202882e+00, -2.43021346e+00,
-2.06654428e-01, 2.77203185e-01, -1.92454320e+00, -3.21138712e+00,
7.81846978e-01, -4.94343116e-01, -1.42285605e+00, 1.01165779e+00,
-1.05621289e+00, -4.70858405e-01, -1.84855538e+00, 1.83544233e+00,
8.54216877e-02, 1.20768462e+00, -5.00878167e-02, -3.01017318e-01,
6.74436614e-02, 3.57643587e+00, 9.73112849e-03, -1.41660006e+00,
9.55491387e-01, 1.55697817e+00, -1.76747086e+00, -8.36885348e-02,
1.32484606e-01, 6.59713537e-01, -4.39207927e-01, 7.29622464e-01,
-8.83102516e-01, -5.68212399e-01, -1.54701409e+00, -1.83484749e+00,
-1.11207542e+00, -6.30133038e-01, 3.59731712e+00, -1.22433911e+00,
4.81430610e-01, 1.69157704e+00, 1.24325615e+00, 1.89304306e+00,
-2.52977962e-01, 4.13964819e-01, -1.29019087e-02, 2.32770647e+00,
-2.18338566e+00, 9.56747474e-01, -7.86999019e-01, 2.72178666e-01,
7.83949548e-01, -6.15045341e-01, 1.81493448e-01, 1.84227027e-01,
6.23226125e-01, -3.13434625e-01, 2.03515829e+00, -3.26083617e-01,
1.33384415e+00, 1.53257334e+00, 2.17861109e-01, -3.25970270e-01,
-1.71850221e-01, 8.62358831e-01, 1.14420483e+00, 1.74114274e+00,
2.76831195e-01, 5.44802610e-01, -7.49641276e-01, -8.43549430e-01,
9.40281558e-01, 2.22689248e+00, -8.61136444e-01, 2.04034863e-01,
-2.88333582e-01, 1.29941573e-01, -9.44808543e-01, -1.07027467e+00,
-1.24412738e+00, 1.20884480e-01, -8.87989794e-01, -5.05801499e-01,
5.71754111e-01, -2.20892129e+00, -1.03869720e+00, -1.19586039e-02,
-5.91241725e-01, -1.11146963e+00, 4.75419608e-01, -2.17156133e-01,
-3.24144476e-01, -8.88137773e-02, 2.89674600e+00, 1.35049104e+00,
9.91113967e-01, 1.20965203e+00, -3.29231573e-01, 4.85379819e-01,
1.29285142e+00, 8.44040495e-03, -9.01242040e-01, 8.34947862e-01,
1.70705398e+00, 9.32748859e-01, 1.52478129e-02, -1.08514648e+00,
7.43950996e-01, 5.20194904e-02, 1.07873251e-01, 6.72791526e-01,
1.03483294e+00, -1.76558135e-01, 6.93114569e-01, 2.91510294e-02,
7.26748970e-01, 4.94061553e-01, 2.16663518e+00, 8.04155766e-01,
8.88196452e-01, 9.53864778e-01, 1.92699086e+00, -1.38556787e-01,
1.51858059e-01, -1.69771866e+00, -1.08174290e+00, 2.49518963e-01,
-2.36819129e-01, 3.91158717e-01, -6.51798125e-01, -1.99902160e+00,
-6.80987615e-01, -1.87441803e-01, 1.44935725e+00, -6.41260947e-01,
-4.99065166e-01, -6.45009215e-01, -9.17364747e-01, -8.55997669e-01,
2.19335212e+00, -5.66664302e-01, -1.97867758e+00, 1.94087914e-01,
1.68614966e-01, -6.60388994e-01, 1.34923590e+00, -1.37132068e-01,
-9.89474077e-01, 2.39157298e+00, 2.04987816e-01, 3.05089529e-01,
1.13792090e+00, -1.41417687e+00, 3.72546926e-01, 8.90006152e-01,
-2.37727905e-01, 7.15436698e-01, 8.74792313e-01, -1.69586708e+00,
1.71054562e+00, -7.88000038e-01, -3.11766786e-01, -1.23965323e+00,
4.68960542e-01, -1.18410977e+00, -1.18461208e+00, 2.04842309e-01,
1.32590028e-01, -8.84822591e-01, -2.34845349e+00, 6.28220250e-01,
-1.30874939e-02, 1.13634500e+00, -2.17638776e+00, -1.22577437e+00,
7.88243135e-01, -1.20009066e+00, 1.91180242e+00, -1.29433083e+00,
2.28446800e-01, -4.50400877e-01, -7.30477576e-01, 3.55069314e-01,
-3.06574408e-01, -1.23074720e+00, 5.51622900e-01, -7.32176553e-01,
9.40304844e-01, 2.12630177e+00, -7.00436261e-01, -1.07417304e+00,
3.65847962e-01, -1.44356346e+00, 4.74709400e-01, -1.06116166e-01,
-6.66407243e-01, 1.90500484e+00, 2.75017564e-01, 4.57761931e-02,
-1.28097188e+00, -9.71959081e-01, -1.07396235e-01, 1.06247882e+00,
3.58182102e-01, 4.22178850e-01, -4.31319965e-01, -4.90864508e-01,
-6.47302201e-01, 2.37596386e+00, -6.96921942e-02, 3.09228602e-01,
1.21870025e+00, -2.70134220e-01, 3.08017056e-01, -5.84184856e-01,
-1.27852639e+00, 4.18588129e-01, -7.45331945e-01, -1.59068294e-01,
2.80776902e-01, 3.46668564e-01, 6.76588944e-01, -1.11051783e+00,
-1.08971194e+00, 1.41879573e-01, -7.63544990e-01, -6.65690087e-01,
-1.76808103e+00, 4.22807571e-01, -4.39634215e-01, 1.24194445e-01,
7.61460961e-01, -3.99892508e-01, 3.93130817e-01, -3.43281082e-02,
-9.86332136e-01, 2.66332391e+00, 9.72884761e-01, -1.18004766e+00,
-1.09838424e+00, 6.40122114e-01, 8.55557360e-01, 3.96390559e-01,
1.94218446e+00, -5.52968777e-01, -1.62050947e-01, -2.65894704e-01,
2.33393748e+00, 6.57008336e-01, 2.58477191e-01, 2.43219079e+00,
-2.32021713e+00, 3.97121111e-01, 9.77122138e-02, -3.13478838e-01,
1.09516389e+00, 5.53267244e-01, 2.70629023e-01, -7.58102582e-01,
-1.66491499e+00, 1.14093117e+00, -7.91930409e-01, -4.91190473e-01,
1.27886657e+00, -5.69314357e-01, 9.06733565e-01, 3.40679223e-01,
6.41147048e-02, -5.72075184e-01, -1.05409977e-02, -1.64904015e+00,
-3.50495148e+00, -1.07449345e+00, -1.07177849e+00, 3.74000256e-01,
6.71714697e-01, 4.29318166e-01, -8.40398760e-01, 9.80742549e-01,
1.14754252e-02, -3.04449258e-01, -6.98813734e-01, 6.64483704e-02,
6.99714715e-01, -7.40769598e-01, 5.16941065e-01, -1.52159277e+00,
7.25639074e-01, 3.01346301e-01, 3.86813048e-01, -3.06866377e-01,
-1.13479221e+00, 1.03880390e+00, -2.71805294e+00, -8.21892545e-01,
-1.26687403e+00, 1.47390986e+00, 3.95502192e-02, -6.74448802e-01,
3.63877950e-01, 3.61230522e-01, 2.41602128e-01, 2.15287848e+00,
1.40127340e+00, -1.51486354e-02, -1.11754724e+00, -1.47355465e-02,
2.09592683e-01, -1.95353546e+00, 5.72710722e-01, 9.78364994e-01,
1.21984164e+00, 6.90469416e-02, -1.35689390e+00, 3.22861905e+00,
3.61823853e+00, -9.87534567e-01, 5.11165110e-01, 5.32794477e-01,
2.08582449e-01, 3.72183451e-01, 2.81463979e-01, 1.07231679e+00,
-7.10901741e-01, 1.78485769e-01, 6.72799179e-02, 2.20681631e-01,
1.22115922e+00, -4.27986412e-01, -7.96295559e-01, -4.14046066e-01,
-3.64119089e-02, -3.10820027e-01, -6.50293560e-01, -6.77101477e-01,
-7.78321108e-02, -7.95070770e-01, 5.94612609e-02, -1.82026772e+00,
-2.37428376e-01, -1.39288448e+00, 6.96976667e-01, 9.02302818e-01,
-1.25895725e+00, 1.28672905e+00, 1.78920416e+00, -2.50457967e-01,
1.99555644e-02, 1.31965662e+00, 3.84139171e-01, -3.28156125e-02,
2.04359887e-01, -2.26492308e+00, -5.52222954e-01, 9.77392353e-01,
-1.36763824e+00, -1.05187782e+00, 2.62200011e-02, 1.04246432e+00,
7.41673361e-02, 3.31257408e-01, -1.59514855e-01, -1.03112795e+00,
-5.70931036e-01, -1.64417057e-01, -5.92448957e-02, 1.99017808e+00,
7.99713912e-01, 5.22006783e-01, 2.39040174e+00, 9.01094668e-02,
-1.89556983e-01, -3.78931260e+00, -6.90240883e-01, 1.36924651e-01,
-1.43466386e+00, -1.03118382e+00, -2.17888811e-01, -7.71372386e-01,
1.34727901e+00, -3.49114017e-01, -1.25712954e+00, -1.09320595e+00,
1.17191204e+00, -5.60330678e-01, -2.02828820e+00, -1.15174936e+00,
-9.09022130e-01, 3.53276860e-01, 4.03861482e-02, 2.07662361e-01,
-2.01047245e+00, -7.84043337e-01, 5.96431590e-01, -2.71363318e+00,
6.48175176e-02, 1.07548132e-01, -8.40236372e-02, -5.88029636e-01,
-1.01843183e+00, -7.75717652e-01, -9.97061999e-01, -2.09938703e-01,
1.36240452e-01, 1.12394611e+00, 3.30243432e-01, -6.80752603e-01,
-2.02732860e+00, -1.07136590e+00, -5.39416801e-01, -9.22231451e-02,
-2.75038929e-01, 1.13240300e+00, 1.80717194e-01, -6.96441382e-01,
-1.34698859e+00, -2.83031805e+00, 3.57892652e-01, -3.12758418e+00,
5.31434015e-01, -6.72184025e-01, 1.40090147e+00, 3.23006925e-01,
2.74605145e-01, 6.28495840e-01, -2.43919087e+00, -3.73581222e-01,
-7.78091987e-01, 3.40000949e-01, -1.64136220e-01, -1.74345397e+00,
-1.01895973e+00, -5.09155078e-02, 1.23765182e+00, -1.01082482e+00,
1.12175979e-01, -2.54944459e+00, -5.99798217e-01, 1.28582406e-01,
3.46721204e-01, 4.15834370e-01, -1.58833310e+00, -1.22441551e+00,
1.87723603e+00, -2.28946794e+00, -1.55878215e+00, 1.01398465e+00,
-1.77897731e+00, -4.51670894e+00, 1.20557656e+00, 3.94199177e-01,
5.74476546e-01, -6.38153861e-01, -6.60106900e-01, 9.75272601e-01,
-6.16124750e-01, 1.03447295e+00, 2.52370756e-01, -3.20622788e-01,
9.01685902e-01, 8.96673400e-01, 1.70066988e-01, 6.68706515e-01,
-1.75531656e+00, 1.02578301e+00, -8.44902227e-01, 9.04132173e-01,
6.85724088e-01, -1.84000820e+00, -4.07477191e-01, -1.82119522e+00,
-1.01164182e-01, -8.15360799e-01, -4.53228821e-01, 2.19557738e-01,
-5.99549119e-01, -2.22593096e-01, 7.15476674e-01, -5.54389695e-01,
-2.76667105e-01, 6.81784255e-01, -9.87076936e-01, 2.51663653e-01,
-8.60889766e-01, -1.62639294e-02, 1.25996146e-01, -6.02397003e-01,
-3.79509914e-01, 1.15651248e+00, -3.70510248e-01, 4.51471255e-01,
8.49687922e-01, 1.52503857e-01, -3.81245200e-01, -3.97668144e-01,
-1.38012519e+00, -3.42556608e-01, -9.65798827e-01, -4.60301848e-02,
4.96586428e-03, 1.49271501e+00, 3.04949487e-01, -5.78192510e-01,
-9.47462912e-02, -2.60891093e-01, -1.62449811e-01, 9.52706107e-01,
-8.53091255e-01, -7.28591504e-01, -1.13119730e+00, 6.30989700e-01,
7.99143599e-01, 2.01677473e+00, -1.14975420e-01, 2.87907952e-01,
-3.77346381e-01, -1.84760885e+00, 1.07076299e+00, -9.91416876e-01,
-9.64751836e-01, 1.60627994e+00, 2.18403904e-01, -1.65083839e+00,
-3.72828414e-02, -9.80715232e-01, -9.13371110e-01, -1.08877786e-01,
-1.22094387e-01, -1.38673390e-01, -5.93839831e-01, 4.54203420e-01,
-1.53949310e+00, -3.91304425e-01, -4.98499095e-01, -9.65375597e-01,
-6.45086720e-01, 2.21712572e-01, 1.47852358e+00, -7.38363353e-01,
3.79611546e-01, -5.03879825e-01, -1.47343248e-01, 5.99098313e-02,
-2.55774579e-01, -6.57828973e-01, -2.48280864e+00, -7.59752834e-01,
4.01349800e-01, 2.23696398e-01, 3.96695750e-01, 6.71896290e-01,
9.34436908e-01, 1.58745228e-01, -8.47427795e-01, -1.79698788e+00,
7.70740713e-01, 9.93491448e-01, 4.86574804e-01, 1.85259905e+00,
6.88065845e-02, 5.97567123e-01, -7.16941778e-01, 1.76719572e+00,
-6.77591123e-01, 1.95433831e+00, 6.42688371e-01, -1.66799497e-02,
4.57572640e-01, 3.95743410e-02, 1.97152412e-01, 1.68384677e+00,
7.17119962e-01, -1.09555242e+00, -1.98426712e-01, -4.98962833e-01,
2.08847976e+00, 9.55048249e-02, 5.34907227e-02, -1.07698637e+00,
2.12435189e-01, 2.42523782e+00, 4.82199454e-01, -1.72664098e+00,
4.40842285e-01, 1.06823313e+00, -4.71414760e-01, -9.89152854e-01,
1.56733198e+00, 7.01238847e-02, 1.64552555e+00, -5.05172246e-01,
-1.53124141e+00, -7.26814216e-01, 2.96880989e+00, -4.57751355e-01,
9.23176296e-01, -1.57182499e+00, -6.10276616e-01, 6.33556483e-01,
-8.29672373e-01, -9.51414725e-02, -1.01612493e+00, -1.65914088e+00,
1.50240690e-02, -4.04605525e-01, 4.54091957e-01, 1.96104806e-01,
1.92194868e+00, -1.56156827e+00, 1.32702829e+00, 3.66661037e-02,
1.26848221e-01, -1.07561704e-02, 9.33018317e-01, -1.52379256e-01])
import pandas as pd
df_x=pd.DataFrame(x)
df_x.describe().T
count | mean | std | min | 25% | 50% | 75% | max | |
---|---|---|---|---|---|---|---|---|
0 | 1000.0 | -0.072891 | 1.081381 | -4.516709 | -0.739328 | -0.050502 | 0.563663 | 3.618239 |
描述性统计
x.mean()
-0.07289058666672985
x.min()
-4.5167089414384085
x.max()
3.618238525954051
x.var()
1.168214980664953
m, v, s, k = stats.t.stats(10, moments='mvsk')
m
array(0.)
v
array(1.25)
s
array(0.)
k
array(1.)
stats.t.stats(10, moments='mvsk')
(array(0.), array(1.25), array(0.), array(1.))
sstr = '%-14s mean = %6.4f, variance = %6.4f, skew = %6.4f, kurtosis = %6.4f'
print(sstr % ('distribution:', m, v, s ,k))
distribution: mean = 0.0000, variance = 1.2500, skew = 0.0000, kurtosis = 1.0000
'''
注意:stats.describe使用无偏估计量作为方差,而np.var是有偏估计量。
对于我们的样本,样本统计量与理论值相差很小。
'''
n, (smin, smax), sm, sv, ss, sk = stats.describe(x)
print(sstr % ('sample:', sm, sv, ss, sk))
sample: mean = -0.0729, variance = 1.1694, skew = 0.0162, kurtosis = 0.8654
T检验和KS检验
print('t-statistic = %6.3f pvalue = %6.4f' % stats.ttest_1samp(x, m))
t-statistic = -2.132 pvalue = 0.0333
tt = (sm-m)/np.sqrt(sv/float(n)) # t-statistic for mean
pval = stats.t.sf(np.abs(tt), n-1)*2 # two-sided pvalue = Prob(abs(t)>tt)
print('t-statistic = %6.3f pvalue = %6.4f' % (tt, pval))
t-statistic = -2.132 pvalue = 0.0333
Kolmogorov-Smirnov检验可用于检验样本来自标准t分布的假设
print('KS-statistic D = %6.3f pvalue = %6.4f' % stats.kstest(x, 't', (10,)))
KS-statistic D = 0.057 pvalue = 0.0027
print('KS-statistic D = %6.3f pvalue = %6.4f' % stats.kstest(x, 'norm'))
KS-statistic D = 0.053 pvalue = 0.0069
d, pval = stats.kstest((x-x.mean())/x.std(), 'norm')
print('KS-statistic D = %6.3f pvalue = %6.4f' % (d, pval))
KS-statistic D = 0.039 pvalue = 0.0871
注意:Kolmogorov-Smirnov检验假设我们针对给定参数的分布进行检验,因为在最后一种情况下,我们估计均值和方差,因此违反了该假设,并且检验统计量的分布(p值是根据,是不正确的。