kernel = [kindex, prior, sprior]
fnames = ["../cache/fith/PI{}.p".format(i) for i in xrange(5)]
statesPI = search.multiPESIPS(ojfa, lb, ub, kernel, 10, fnames)
x = []
y = []
for stateEI in statesEI:
a[2].plot(stateEI[11], 'r')
#a[3].plot([sum(stateEI[5][:i]) for i in xrange(len(stateEI[5]))],stateEI[11],'r')
x.append([sum(stateEI[5][:i]) for i in xrange(len(stateEI[5]))])
y.append(stateEI[11].flatten())
print 'reccomended under EI: {} : {}'.format(
[10**i for i in stateEI[4][-1]],
ojf(stateEI[4][-1], None, None)[0])
X_, Y_, lb_, ub_ = OPTutils.mergelines(x, y)
a[3].fill_between(X_,
lb_,
ub_,
facecolor='lightcoral',
edgecolor='lightcoral',
alpha=0.5)
a[3].plot(X_, Y_, 'r')
x = []
y = []
for statePE in statesPE:
a[2].plot(statePE[11], 'b')
#a[3].plot([sum(statePE[5][:i]) for i in xrange(len(statePE[5]))],statePE[11],'b')
x.append([sum(statePE[5][:i]) for i in xrange(len(statePE[5]))])
y.append(statePE[11].flatten())
lsu = 3
lsl = -7
names = ["../cache/braninnoise/PESVS_"+str(int(100*sp.log10(s)))+"_"+str(pwr)+"_"+str(i)+"_"+str(lsl)+"_"+str(lsu)+".p" for i in xrange(nreps)]
results = search.multiPESVS(ojf,lb,ub,kernel,s,bd,lambda x,s:cfn(s),lsl,lsu,names)
Rz = [sp.log10(sp.array(r[11])-braninmin).flatten() for r in results]
Cz = [[sum(r[5][:j]) for j in xrange(len(r[5]))] for r in results]
#for i in xrange(nreps):
# a[2].plot(Cz[i],Rz[i].flatten(),'rx-')
#[a[1].semilogy([1e-6/(c**(1./pwr)) for c in r[5]],'cyan') for r in results]
x=[]
y=[]
for r in results:
y.append(sp.array(sp.log10([1e-6/(c**(1./pwr)) for c in r[5]])))
print y
x.append(range(len(r[5])))
X_,Y_,lb_,ub_ = OPTutils.mergelines(x,y)
a[1].fill_between(X_,lb_,ub_,facecolor='lightcoral',edgecolor='lightcoral',alpha=0.5)
a[1].plot(X_,Y_,'r')
#a[2].plot(sp.array([sum(C2[:j]) for j in xrange(len(C2))]).flatten(),(sp.log(Yreg)).flatten(),'ro-')
def linint(i,c):
if c>Cz[i][-1]:
return sp.NaN
else:
j=0
while Cz[i][j]<c:
j+=1
return Rz[i][j-1]+(c-Cz[i][j-1])*(Rz[i][j]-Rz[i][j-1])/(Cz[i][j]-Cz[i][j-1])
#for i in xrange(nreps):
#!/usr/bin/env python2
#encoding: UTF-8
# To change this license header, choose License Headers in Project Properties.
# To change this template file, choose Tools | Templates
# and open the template in the editor.
import scipy as sp
from matplotlib import pyplot as plt
n=12
def f(x):
return 0.1*x**2-0.5*x+1.+0.1*(0.2*abs(x)+1)*sp.random.normal()
r=[]
y=[]
for i in xrange(n):
r.append(sorted(sp.random.uniform(0,10,40)))
y.append(map(f,r[-1]))
plt.plot(r[-1],y[-1],'c')
import OPTutils
X,Y,lb,ub = OPTutils.mergelines(r,y)
plt.fill_between(X,lb,ub,facecolor='lightgreen')
plt.plot(X,Y)
plt.show()
#encoding: UTF-8
# To change this license header, choose License Headers in Project Properties.
# To change this template file, choose Tools | Templates
# and open the template in the editor.
import scipy as sp
from matplotlib import pyplot as plt
n = 12
def f(x):
return 0.1 * x**2 - 0.5 * x + 1. + 0.1 * (0.2 * abs(x) +
1) * sp.random.normal()
r = []
y = []
for i in xrange(n):
r.append(sorted(sp.random.uniform(0, 10, 40)))
y.append(map(f, r[-1]))
plt.plot(r[-1], y[-1], 'c')
import OPTutils
X, Y, lb, ub = OPTutils.mergelines(r, y)
plt.fill_between(X, lb, ub, facecolor='lightgreen')
plt.plot(X, Y)
plt.show()
