In [57]:

    

# Import all libraries needed for the tutorial

# General syntax to import specific functions in a library:
##from (library) import (specific library function)
from pandas import DataFrame, read_csv

# General syntax to import a library but no functions:
##import (library) as (give the library a nickname/alias)
import matplotlib.pyplot as plt
import pandas as pd #only needed to determine version number

import glob
from scipy import randn
from pandas import *

# Enable inline plotting
%matplotlib inline


    

In [58]:

    

path =r'../expdist-1-5range'

allFiles = glob.glob(path + "/*.txt")
frame = pd.DataFrame()
list = []
for files in allFiles:
    df = pd.read_csv(files,index_col=0, header=0, delimiter="\t")
    list.append(df)


glued = pd.concat(list, axis=1)
ar = glued
#print(frame)

#print(glued)
#print(glued.groupby(level=0, axis=1).std())
stdtab = ar.groupby(level=0, axis=1).std()
meantab = glued.groupby(level=0, axis=1).mean()
#glued.groupby(level=0, axis=1).mean().to_csv("glued.txt", sep="\t")

#plt.figure(figsize=(10, 8), dpi=120);
#meantab.plot(yerr=stdtab)

#dff = pd.concat([meantab, stdtab], axis=1)

#plt.figure(figsize=(10, 8), dpi=120);
#dff.plot()
list = []
for files in allFiles:
    df = pd.read_csv(files,index_col=None, header=0, delimiter="\t")
    list.append(df)
frame = pd.concat(list)

#plt.figure(figsize=(10, 8), dpi=120);
#bp = frame.boxplot(column=['# Clock'])
#fig = plt.gcf()
#fig.savefig(path + '/clockboxplot.pdf')
#plt.figure(figsize=(10, 8), dpi=120);
#bp = frame.boxplot(column=['In Chain'])
#fig = plt.gcf()
#fig.savefig(path + '/inchainboxplot.pdf')
#plt.figure(figsize=(10, 8), dpi=120);
#bp = frame.boxplot(column=['On Target'])
#fig = plt.gcf()
#fig.savefig(path + '/ontargetboxplot.pdf')


    

In [60]:

    

#xticks=arange(0, 9000, 1000)
meantab[::100].plot(yerr=glued.groupby(level=0, axis=1).std(),figsize=(15, 12))

#fig = plt.gcf()
#fig.savefig(path + '/figure.pdf')


    

    
Out[60]:





<matplotlib.axes.AxesSubplot at 0x10c34db90>






    

In [61]:

    

path =r'../expdist2-5'


    

In [67]:

    

allFiles = glob.glob(path + "/*.txt")
frame = pd.DataFrame()
list = []
for files in allFiles:
    df = pd.read_csv(files,index_col=0, header=0, delimiter="\t")
    list.append(df)


glued = pd.concat(list, axis=1)
ar = glued
#print(frame)

#print(glued)
#print(glued.groupby(level=0, axis=1).std())
stdtab2 = ar.groupby(level=0, axis=1).std()
meantab2 = glued.groupby(level=0, axis=1).mean()
#glued.groupby(level=0, axis=1).mean().to_csv("glued.txt", sep="\t")

#plt.figure(figsize=(10, 8), dpi=120);
#meantab.plot(yerr=stdtab)

#dff = pd.concat([meantab, stdtab], axis=1)

#plt.figure(figsize=(10, 8), dpi=120);
#dff.plot()
list = []
for files in allFiles:
    df = pd.read_csv(files,index_col=None, header=0, delimiter="\t")
    list.append(df)
frame2 = pd.concat(list)

meantab.rename(columns={'In Chain':'In Chain Range 1.5'}, inplace=True)
meantab.rename(columns={'On Target':'On Target Range 1.5'}, inplace=True)
stdtab.rename(columns={'In Chain':'In Chain Range 1.5'}, inplace=True)
stdtab.rename(columns={'On Target':'On Target Range 1.5'}, inplace=True)

meantab2.rename(columns={'In Chain':'In Chain Range 2.5'}, inplace=True)
meantab2.rename(columns={'On Target':'On Target Range 2.5'}, inplace=True)
stdtab2.rename(columns={'In Chain':'In Chain Range 2.5'}, inplace=True)
stdtab2.rename(columns={'On Target':'On Target Range 2.5'}, inplace=True)


    

In [68]:

    

med = pd.concat([meantab, meantab2], axis=1)


    

In [69]:

    

error = pd.concat([stdtab, stdtab2], axis=1)


    

In [70]:

    

med.plot()
error.plot()


    

    
Out[70]:





<matplotlib.axes.AxesSubplot at 0x10ff75e90>






    













    

In [74]:

    

med[::100].plot(yerr=pd.concat([stdtab, stdtab2], axis=1),figsize=(15, 12)).set_title('Range comparison')

fig = plt.gcf()
fig.savefig(path + '/lowerversus.pdf')


    

In [ ]: