def test1():
"""Testing if python and C++ implementation of Seigel yield the same results"""
check = True
for _ in range(100):
se = tseries.SeigelEstimator()
sse = tseries.SpanSeigelEstimator(20)
cse = cgoals.cSeigelEstimator(20)
size = random.randint(2, 20)
data = [random.uniform(-10, 10) for _ in xrange(size)]
for i, d in enumerate(data):
cse.add(d)
sse.add(d)
se.add(d)
if abs(cse.value() - se.value()) > 0.01:
check = False
print cse.value()
print sse.value()
print se.value()
print data[:(i+1)]
return check
def test1():
"""Testing if python and C++ implementation of Seigel yield the same results"""
check = True
for _ in range(100):
se = tseries.SeigelEstimator()
sse = tseries.SpanSeigelEstimator(20)
cse = cgoals.cSeigelEstimator(20)
size = random.randint(2, 20)
data = [random.uniform(-10, 10) for _ in xrange(size)]
for i, d in enumerate(data):
cse.add(d)
sse.add(d)
se.add(d)
if abs(cse.value() - se.value()) > 0.01:
check = False
print cse.value()
print sse.value()
print se.value()
print data[:(i + 1)]
return check
import testenv
import random
import time
random.seed(123)
import goals.explorer.tools.tseries as tseries
import cgoals
tse = tseries.TheilSenEstimator()
se = tseries.SeigelEstimator()
sse = tseries.SpanSeigelEstimator(20)
cse = cgoals.cSeigelEstimator(20)
# preparing data
data = []
for i in range(1000):
if i == 4:
v = random.uniform(0, 2000)
else:
v = i + random.uniform(-2, 2)
data.append(v)
t0 = time.time()
data0 = [0.0, 1.0, 2.0]
for d in data:
data0.append(d)
tseries.derivative(data0)
t0 = time.time() - t0
import testenv
import random
import time
random.seed(123)
import goals.explorer.tools.tseries as tseries
import cgoals
tse = tseries.TheilSenEstimator()
se = tseries.SeigelEstimator()
sse = tseries.SpanSeigelEstimator(20)
cse = cgoals.cSeigelEstimator(20)
# preparing data
data = []
for i in range(1000):
if i == 4:
v = random.uniform(0, 2000)
else:
v = i + random.uniform(-2, 2)
data.append(v)
t0 = time.time()
data0 = [0., 1., 2.]
for d in data:
data0.append(d)
tseries.derivative(data0)
t0 = time.time() - t0
