Maturation time update

Maturation time 1/km + tauM

rcode/distributedDelayvsdelayequation.R

@@ -19,7 +19,8 @@ idx = which.min(tmpCsv$norm)
 concatDF = tmpCsv[idx,]
 concatDF['dtauM'] = 0
 concatDF['idx'] = 0
-inData = rbind(inData,  concatDF[,c('idx', 'norm',  'tauM', 'dtauM')])
+concatDF['TM'] = 1/concatDF['kM']+ concatDF['tauM']
+inData = rbind(inData,  concatDF[,c('idx', 'norm',  'tauM', 'dtauM', 'TM')])
 
 
 
@@ -32,7 +33,8 @@ for (idx in seq(1,length(delay),1)) {
   concatDF = tmpCsv[idxC,]
   concatDF['idx'] = idx
   concatDF['dtauM'] = concatDF['dtauM']/2
-  inData = rbind(inData, concatDF[,c('idx','norm',  'tauM', 'dtauM')])
+  concatDF['TM'] = 1/concatDF['kM']+ concatDF['tauM']
+  inData = rbind(inData, concatDF[,c('idx','norm',  'tauM', 'dtauM', 'TM')])
 }
 
 
@@ -58,7 +60,7 @@ ggsave(filename="distributedDelay_chi2.pdf", plot=pChi2, width = 71, height = 30
 
 # ---- 
 
-ptaum<- ggplot(inData, aes(x = dtauM, y = tauM, fill = idx)) + geom_bar(stat = "identity", position = 'dodge')  +
+ptaum<- ggplot(inData, aes(x = dtauM, y = TM, fill = idx)) + geom_bar(stat = "identity", position = 'dodge')  +
   theme_classic(base_size = 8) + xlab('Half-width distribution (min)')  +   scale_fill_manual(values = cols[seq(6,1,-1)]) + ylab('Avg. Mat. time (min)\n')  + 
   coord_cartesian(ylim= c(0,50), xlim = c(0.4,6.6), expand = c(0)) + 
   theme(legend.position = 'none', axis.line.x = element_line(color="black", size = 0.5),