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]])
