def dataReader2(dataFile): generation = 'single colonies' colonyCounts = {} dilutionFactors = {} with open(dataFile, 'r') as f: for line in f: vector = line.split(',') if vector[0] == 'final': if vector[1] == 'dilution factor': a = int(vector[2].split('in')[0]) b = int(vector[2].split('in')[1]) c = int(vector[3].split('in')[0]) d = int(vector[3].split('in')[1]) dilutionFactors[generation] = [((a + b) / a)**3, ((c + d) / c)**3] elif 'plate' in vector[2]: if vector[1] != '': colonyName = vector[1] if 'colony 30' not in colonyName: if (colonyName == 'colony 8.1' and vector[2] == 'plate # 3') is False: a = int(vector[3]) b = int(vector[4]) c = int(vector[5]) d = int(vector[6]) if colonyName not in colonyCounts.keys(): colonyCounts[colonyName] = [[a], [b], [c], [d]] else: colonyCounts[colonyName][0].append(a) colonyCounts[colonyName][1].append(b) colonyCounts[colonyName][2].append(c) colonyCounts[colonyName][3].append(d) # converting colony counts into cf validNames = list(colonyCounts.keys()) validNames.sort() converted = {} for name in validNames: print(name) treatment = [colonyCounts[name][2], colonyCounts[name][3]] nt = [colonyCounts[name][0], colonyCounts[name][1]] cf_mu, cf_sd, pvalue = calculateStatistics.main(treatment, nt) converted[name] = [cf_mu, cf_sd, pvalue] return converted
import matplotlib,numpy,sys,scipy import matplotlib.pyplot sys.path.append('/Users/alomana/gDrive2/projects/centers/ap/src/assessmentGraphs/publicationFigures/lib') import calculateStatistics ### MAIN matplotlib.rcParams.update({'font.size':36,'font.family':'Times New Roman','xtick.labelsize':28,'ytick.labelsize':28}) thePointSize=12 # clonal 2 xSignal=numpy.array([[187, 144, 166, 171, 158],[77, 78, 60, 61, 49]]) xNoSignal=numpy.array([[170, 143, 163, 141, 162],[56, 35, 45, 73, 42]]) cf_mu_0, cf_sd_0, pvalue_0 = calculateStatistics.main(xSignal, xNoSignal) xSignal=numpy.array([[200, 171, 173, 176, 199],[174, 176, 147, 186, 197]]) xNoSignal=numpy.array([[203, 218, 231, 210, 210],[124, 125, 130, 161, 145]]) cf_mu_50, cf_sd_50, pvalue_50 = calculateStatistics.main(xSignal, xNoSignal) xSignal=numpy.array([[153, 154, 124, 167, 157],[126, 161, 115, 114, 110]]) xNoSignal=numpy.array([[187, 179, 175, 160, 150],[127, 151, 115, 133, 89]]) cf_mu_100, cf_sd_100, pvalue_100 = calculateStatistics.main(xSignal, xNoSignal) xSignal=numpy.array([[188, 207, 186, 183, 193],[155, 162, 161, 160, 145]]) xNoSignal=numpy.array([[202, 185, 237, 217, 186],[136, 162, 139, 151, 146]]) cf_mu_150, cf_sd_150, pvalue_150 = calculateStatistics.main(xSignal, xNoSignal) xSignal=numpy.array([[183, 202, 205, 214, 213],[200, 144, 185, 159, 168]]) xNoSignal=numpy.array([[255, 221, 218, 213, 235],[221, 201, 185, 176, 166]]) cf_mu_200, cf_sd_200, pvalue_200 = calculateStatistics.main(xSignal, xNoSignal)
sys.path.append("/Users/alomana/gDrive2/projects/centers/ap/src/assessmentGraphs/publicationFigures/lib") import calculateStatistics ### MAIN matplotlib.rcParams.update( {"font.size": 36, "font.family": "Times New Roman", "xtick.labelsize": 28, "ytick.labelsize": 28} ) thePointSize = 12 # clonal 2 xSignal = numpy.array([[175, 153, 186, 189, 157], [37, 59, 46, 67, 70]]) xNoSignal = numpy.array([[200, 202, 224, 194, 193], [71, 66, 71, 87, 60]]) cf_mu_0, cf_sd_0, pvalue_0 = calculateStatistics.main(xSignal, xNoSignal) xSignal = numpy.array([[25, 28, 19, 18, 16], [0, 9, 4, 9, 1]]) xNoSignal = numpy.array([[24, 16, 29, 17, 23], [4, 7, 5, 3, 4]]) cf_mu_50, cf_sd_50, pvalue_50 = calculateStatistics.main(xSignal, xNoSignal) xSignal = numpy.array([[96, 97, 94, 127, 80], [32, 36, 36, 42, 36]]) xNoSignal = numpy.array([[104, 137, 110, 128, 113], [52, 36, 32, 50, 41]]) cf_mu_100, cf_sd_100, pvalue_100 = calculateStatistics.main(xSignal, xNoSignal) xSignal = numpy.array([[204, 223, 199, 249, 193], [141, 131, 125, 154, 139]]) xNoSignal = numpy.array([[171, 217, 240, 200, 168], [166, 192, 163, 196, 170]]) cf_mu_150, cf_sd_150, pvalue_150 = calculateStatistics.main(xSignal, xNoSignal) xSignal = numpy.array([[197, 216, 224, 219, 208], [181, 182, 186, 179, 116]]) xNoSignal = numpy.array([[261, 227, 229, 188, 236], [179, 169, 174, 183, 164]])