python-期货指数的单位根检验

代码如下：

             测试结果表明，以下33个品种不存在单位根.

# -*- 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
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/数据/'+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)
x = np.array(r)
result = ts.adfuller(x, 1) # maxlag is now set to 1

print(i)
if result[0]<result[4]['1%']:
print ('没有单位根')
j=j+1

else:
print('由单位跟')

print(j)


结果：
cu
没有单位根
1
al
没有单位根
2
zn
没有单位根
3
pb
没有单位根
4
sn
没有单位根
5
au
没有单位根
6
ag
没有单位根
7
rb
没有单位根
8
hc
没有单位根
9
bu
没有单位根
10
ru
没有单位根
11
m9
没有单位根
12
y9
没有单位根
13
a9
没有单位根
14
p9
没有单位根
15
c9
没有单位根
16
cs
没有单位根
17
jd
没有单位根
18
l9
没有单位根
19
v9
没有单位根
20
pp
没有单位根
21
j9
没有单位根
22
jm
没有单位根
23
i9
没有单位根
24
sr
没有单位根
25
cf
没有单位根
26
zc
没有单位根
27
fg
没有单位根
28
ta
没有单位根
29
ma
没有单位根
30
oi
没有单位根
31
rm
没有单位根
32
sm
没有单位根
33