In [4]:
# import some modules using standard naming conventions
from pylab import *
import pandas as pd
# find the .csv files and read them into pandas data sets
import glob
flist = glob.glob('/usgs/data1/csherwood/poking_eyeball/*.csv')
In [5]:
size(flist)
ma = zeros(size(flist))
sa = zeros(size(flist))
mg = zeros(size(flist))
i=0
for f in flist:
exec( 'd'+f[-6:-4]+' = pd.read_csv(f)' )
exec( 'ma[i]=mean(d'+f[-6:-4]+'.arith_mean)' )
exec( 'sa[i]=std(d'+f[-6:-4]+'.arith_mean)' )
exec( 'mg[i]=mean(d'+f[-6:-4]+'.geom_mean)' )
#print f, ma[i], sa[i]
i=i+1
In [6]:
def Phi(mm):
return -1.44296504*log(mm)
phi = array([0, .5, 1., 1.5, 2., 2.5, 3., 3.5, 4]);
phima = Phi(ma)
phisa = Phi(sa)
phimg = Phi(mg)
figure()
plot(phi,phi,'-k')
plot(phi-0.5,phi,'--k') # one sieve larger
plot(phi[0]*ones_like(d00.geom_mean),Phi(d00.geom_mean),'o',color='0.8')
plot(phi[1]*ones_like(d05.geom_mean),Phi(d05.geom_mean),'o',color='0.8')
plot(phi[2]*ones_like(d10.geom_mean),Phi(d10.geom_mean),'o',color='0.8')
plot(phi[3]*ones_like(d15.geom_mean),Phi(d15.geom_mean),'o',color='0.8')
plot(phi[4]*ones_like(d20.geom_mean),Phi(d20.geom_mean),'o',color='0.8')
plot(phi[5]*ones_like(d25.geom_mean),Phi(d25.geom_mean),'o',color='0.8')
plot(phi[6]*ones_like(d30.geom_mean),Phi(d30.geom_mean),'o',color='0.8')
plot(phi[7]*ones_like(d35.geom_mean),Phi(d35.geom_mean),'o',color='0.8')
plot(phi[8]*ones_like(d40.geom_mean),Phi(d40.geom_mean),'o',color='0.8')
plot(phi,phimg,'or')
xlim([-.5,4.5])
xlabel('Sieve size (phi)')
ylabel('DGS geometric mean (phi)')
title('MOF Poking Eyeball Tests, May 3, 2013')
#plt.savefig('D:\crs\proj\poking_eyeball_tests\MOF_test_050313\compare.png')
Out[6]:
In [7]:
flist
Out[7]:
In [7]: