Esempio n. 1
0
if len(argv) < 2:
    print Help()

excelfile = sys.argv[1]
wellcontent_file1 = sys.argv[2]
wellcontent_file2 = sys.argv[3]
outputfile = sys.argv[4]


selection_gradient = WellContent(wellcontent_file1)
inducer_gradient = WellContent(wellcontent_file2)
#print selection_gradient.conc_at('A1')
#print inducer_gradient.conc_at('A1')

expt1 = Experiment(excelfile)
expt1.extract_all_timecourses()
print expt1.all_timecourses

for keys in expt1.all_timecourses:
    print keys, expt1.all_timecourses[keys].data()
exit()
#print expt1.total_timepoints()
#print expt1.total_timecourses()
#print expt1.extract_timecourse('A1')
#A1 = TimeCourse('A1',expt1.extract_timecourse('A1'))
#A1.display()
#print A1.mean()

def sigmoid(x, x0, k):
    y = 1/(1 + np.exp(-k*(x-x0)))
    return y