新聞中心
這里有您想知道的互聯(lián)網(wǎng)營銷解決方案
python實現(xiàn)檢驗33品種數(shù)據(jù)是否是正態(tài)分布-創(chuàng)新互聯(lián)
我就廢話不多說了,直接上代碼吧!
# -*- coding: utf-8 -*- """ Created on Thu Jun 22 17:03:16 2017 @author: yunjinqi E-mail:yunjinqi@qq.com Differentiate yourself in the world from anyone else. """ import pandas as pd import numpy as np import matplotlib.pyplot as plt import statsmodels.tsa.stattools as ts import statsmodels.api as sm from statsmodels.graphics.api import qqplot from statsmodels.sandbox.stats.runs import runstest_1samp import scipy.stats as sts namelist=['cu','al','zn','pb','sn','au','ag','rb','hc','bu','ru','m9','y9','a9', 'p9','c9','cs','jd','l9','v9','pp','j9','jm','i9','sr','cf', 'zc','fg','ta','ma','oi','rm','sm'] j=0 for i in namelist: filename='C:/Users/HXWD/Desktop/數(shù)據(jù)/'+i+'.csv' data=pd.read_csv(filename,encoding='gbk') data.columns=['date','open','high','low','close','amt','opi'] data.head() data=np.log(data['close']) r=data-data.shift(1) r=r.dropna() #print(r) rate = np.array(list(r)) print('品種{}數(shù)據(jù)長度{}均值{}標(biāo)準(zhǔn)差{}方差{}偏度{}峰度{}'.format(i,len(rate), rate.mean(),rate.std(),rate.var(),sts.skew(rate), sts.kurtosis(rate)))
當(dāng)前題目:python實現(xiàn)檢驗33品種數(shù)據(jù)是否是正態(tài)分布-創(chuàng)新互聯(lián)
分享URL:http://www.ef60e0e.cn/article/dggoss.html