diff --git a/results/fitnup107nup358/.RData b/results/fitnup107nup358/.RData new file mode 100644 index 0000000000000000000000000000000000000000..b28f31e878e247bd8312ca952493a18f60184a60 Binary files /dev/null and b/results/fitnup107nup358/.RData differ diff --git a/results/fitnup107nup358/.Rhistory b/results/fitnup107nup358/.Rhistory new file mode 100644 index 0000000000000000000000000000000000000000..2a3e3a68785c34a0d255501b7531da5345143eb9 --- /dev/null +++ b/results/fitnup107nup358/.Rhistory @@ -0,0 +1,512 @@ +'lower' = datain$NUP358_C - datain$NUP358_C_std) +plotDataAll(nc107, c107, simulNUP107, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP107.pdf') +plotDataAll(nc358, c358, simulNUP358, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP358.pdf') +plotDataF(nc107, simulNUP107, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray40'), 0.4, 'NUP107_NC.pdf') +plotDataF(c107, simulNUP107, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray40'), 0.4, 'NUP107_C.pdf') +plotDataF(nc358, simulNUP358, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray40' ), 0.4, 'NUP358_NC.pdf') +plotDataF(c358, simulNUP358, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray40'), 0.4, 'NUP358_C.pdf') +plotData_IP(simulNUP107, simulNUP358, c('time','C_ip_density'), c('#A2C964', '#C1FF00'), 'NUP107_358_C_den.pdf' ) +# ----- +plotData_IP <-function(simul1, simul2, names, cols, fname) { +pdf(fname, width = 2.8, height=2, useDingbats = F,pointsize = 8, paper="special", fonts="Helvetica") +lw<- 2 +par(mar=c(4, 4, 1,1), mgp=c(2.5,1,0)) +plot(simul1[,names[1]], simul1[,names[2]], xlab = "", ylab = "Normalized interphase density", +type = "n", yaxt = 'n', xaxt = 'n', xlim = c(0, 120), ylim = c(0, 1.1), bty='l') +axis(1, tck = -0.05, at=c(0,20,40,60,80,100,120)) +axis(2, las=2, tck =-0.05, at=c(0,0.2,0.4,0.6,0.8,1)) +mtext(side=1, text="Time after AO (min)", line=2) +lines(simul1[,names[1]], simul1[,names[2]], lwd = lw, lty = 1, col = cols[1]) +lines(simul2[,names[1]], simul2[,names[2]], lwd = lw, lty = 1, col = cols[2]) +par(xpd=TRUE) +legend(0,1.1, c("Nup107", "Nup358"), +lwd = c(1.5,1.5), +col=c(cols[1], cols[2]), horiz = F, bty='n', +y.intersp = 0.8, cex = 0.8) +dev.off() +} +# ---- +nc107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_NC, 'upper'=datain$NUP107_NC + datain$NUP107_NC_std, +'lower' = datain$NUP107_NC - datain$NUP107_NC_std) +c107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_C, 'upper'=datain$NUP107_C + datain$NUP107_C_std, +'lower' = datain$NUP107_C - datain$NUP107_C_std) +nc358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_NC, 'upper'=datain$NUP358_NC + datain$NUP358_NC_std, +'lower' = datain$NUP358_NC - datain$NUP358_NC_std) +c358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_C, 'upper'=datain$NUP358_C + datain$NUP358_C_std, +'lower' = datain$NUP358_C - datain$NUP358_C_std) +plotDataAll(nc107, c107, simulNUP107, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP107.pdf') +plotDataAll(nc358, c358, simulNUP358, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP358.pdf') +plotDataF(nc107, simulNUP107, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray40'), 0.4, 'NUP107_NC.pdf') +plotDataF(c107, simulNUP107, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray40'), 0.4, 'NUP107_C.pdf') +plotDataF(nc358, simulNUP358, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray40' ), 0.4, 'NUP358_NC.pdf') +plotDataF(c358, simulNUP358, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray40'), 0.4, 'NUP358_C.pdf') +plotData_IP(simulNUP107, simulNUP358, c('time','C_ip_density'), c('green', 'yellow'), 'NUP107_358_C_den.pdf' ) +#c('#A2C964', '#C1FF00') +# ---- +nc107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_NC, 'upper'=datain$NUP107_NC + datain$NUP107_NC_std, +'lower' = datain$NUP107_NC - datain$NUP107_NC_std) +c107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_C, 'upper'=datain$NUP107_C + datain$NUP107_C_std, +'lower' = datain$NUP107_C - datain$NUP107_C_std) +nc358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_NC, 'upper'=datain$NUP358_NC + datain$NUP358_NC_std, +'lower' = datain$NUP358_NC - datain$NUP358_NC_std) +c358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_C, 'upper'=datain$NUP358_C + datain$NUP358_C_std, +'lower' = datain$NUP358_C - datain$NUP358_C_std) +plotDataAll(nc107, c107, simulNUP107, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP107.pdf') +plotDataAll(nc358, c358, simulNUP358, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP358.pdf') +plotDataF(nc107, simulNUP107, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray40'), 0.4, 'NUP107_NC.pdf') +plotDataF(c107, simulNUP107, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray40'), 0.4, 'NUP107_C.pdf') +plotDataF(nc358, simulNUP358, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray40' ), 0.4, 'NUP358_NC.pdf') +plotDataF(c358, simulNUP358, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray40'), 0.4, 'NUP358_C.pdf') +plotData_IP(simulNUP107, simulNUP358, c('time','C_ip_density'), c('#A2C964', 'yellow'), 'NUP107_358_C_den.pdf' ) +#c('#A2C964', '#C1FF00') +# ---- +nc107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_NC, 'upper'=datain$NUP107_NC + datain$NUP107_NC_std, +'lower' = datain$NUP107_NC - datain$NUP107_NC_std) +c107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_C, 'upper'=datain$NUP107_C + datain$NUP107_C_std, +'lower' = datain$NUP107_C - datain$NUP107_C_std) +nc358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_NC, 'upper'=datain$NUP358_NC + datain$NUP358_NC_std, +'lower' = datain$NUP358_NC - datain$NUP358_NC_std) +c358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_C, 'upper'=datain$NUP358_C + datain$NUP358_C_std, +'lower' = datain$NUP358_C - datain$NUP358_C_std) +plotDataAll(nc107, c107, simulNUP107, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP107.pdf') +plotDataAll(nc358, c358, simulNUP358, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP358.pdf') +plotDataF(nc107, simulNUP107, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray40'), 0.4, 'NUP107_NC.pdf') +plotDataF(c107, simulNUP107, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray40'), 0.4, 'NUP107_C.pdf') +plotDataF(nc358, simulNUP358, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray40' ), 0.4, 'NUP358_NC.pdf') +plotDataF(c358, simulNUP358, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray40'), 0.4, 'NUP358_C.pdf') +plotData_IP(simulNUP107, simulNUP358, c('time','C_ip_density'), c('#A2C964', '#C1FF00'), 'NUP107_358_C_den.pdf' ) +#c('#A2C964', '#C1FF00') +# ---- +nc107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_NC, 'upper'=datain$NUP107_NC + datain$NUP107_NC_std, +'lower' = datain$NUP107_NC - datain$NUP107_NC_std) +c107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_C, 'upper'=datain$NUP107_C + datain$NUP107_C_std, +'lower' = datain$NUP107_C - datain$NUP107_C_std) +nc358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_NC, 'upper'=datain$NUP358_NC + datain$NUP358_NC_std, +'lower' = datain$NUP358_NC - datain$NUP358_NC_std) +c358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_C, 'upper'=datain$NUP358_C + datain$NUP358_C_std, +'lower' = datain$NUP358_C - datain$NUP358_C_std) +plotDataAll(nc107, c107, simulNUP107, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP107.pdf') +plotDataAll(nc358, c358, simulNUP358, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP358.pdf') +plotDataF(nc107, simulNUP107, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray40'), 0.4, 'NUP107_NC.pdf') +plotDataF(c107, simulNUP107, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray40'), 0.4, 'NUP107_C.pdf') +plotDataF(nc358, simulNUP358, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray40' ), 0.4, 'NUP358_NC.pdf') +plotDataF(c358, simulNUP358, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray40'), 0.4, 'NUP358_C.pdf') +plotData_IP(simulNUP107, simulNUP358, c('time','C_ip_density'), c('#A2C964', 'yellow'), 'NUP107_358_C_den.pdf' ) +#c('#A2C964', '#C1FF00') +# ----- +plotDataF <-function(d, simul, names, cols, alpha, fname) { +# function to plot data no separation between interphase and postmitotic +pdf(fname, width = 2.8, height=2, useDingbats = F,pointsize = 8, paper="special", fonts="Helvetica") +lw<- 2 +par(mar=c(4, 4, 1,1), mgp=c(2.5,1,0)) +plot(d$x, d$y, xlab = "", ylab = "Normalized total intensity", +type = "n", yaxt = 'n', xaxt = 'n', xlim = c(0, 120), ylim = c(0, 1.4), bty='l') +axis(1, tck = -0.05, at=c(0,20,40,60,80,100,120)) +axis(2, las=2, tck =-0.05, at=c(0,0.2,0.4,0.6,0.8,1,1.2)) +mtext(side=1, text="Time after AO (min)", line=2) +cols1_t <-adjustcolor(cols[1], alpha.f = alpha) +polygon(c(d$x, rev(d$x)), c(d$upper, rev(d$lower)), col = cols1_t, border = NA) +points(d$x, d$y, pch = 17, col = cols[1], cex=0.6) +colfun <- colorRampPalette(c(cols[1], cols[2])) +lines(simul[,names[1]], simul[,names[2]], lwd = lw, col = colfun(50)[30]) +lines(simul[,names[1]], simul[,names[3]], lwd = lw, lty = 2, col = cols[2]) +lines(simul[,names[1]], simul[,names[4]], lwd = lw, col=cols[2]) +par(xpd=TRUE) +legend(0,1.65, c("data", "post-mitotic+interphase", "post-mitotic", "interphase"), +lwd = c(0,1.5,1.5,1.5),pch = c(17,-1,-1,-1), lty = c(0,1,2,1), +col=c(cols[1], colfun(50)[30], cols[2], cols[2]), horiz = F, bty='n', +y.intersp = 0.8, cex = 0.8) +dev.off() +} +# ---- +nc107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_NC, 'upper'=datain$NUP107_NC + datain$NUP107_NC_std, +'lower' = datain$NUP107_NC - datain$NUP107_NC_std) +c107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_C, 'upper'=datain$NUP107_C + datain$NUP107_C_std, +'lower' = datain$NUP107_C - datain$NUP107_C_std) +nc358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_NC, 'upper'=datain$NUP358_NC + datain$NUP358_NC_std, +'lower' = datain$NUP358_NC - datain$NUP358_NC_std) +c358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_C, 'upper'=datain$NUP358_C + datain$NUP358_C_std, +'lower' = datain$NUP358_C - datain$NUP358_C_std) +plotDataAll(nc107, c107, simulNUP107, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP107.pdf') +plotDataAll(nc358, c358, simulNUP358, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP358.pdf') +plotDataF(nc107, simulNUP107, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray40'), 0.4, 'NUP107_NC.pdf') +plotDataF(c107, simulNUP107, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray40'), 0.4, 'NUP107_C.pdf') +plotDataF(nc358, simulNUP358, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray40' ), 0.4, 'NUP358_NC.pdf') +plotDataF(c358, simulNUP358, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray40'), 0.4, 'NUP358_C.pdf') +plotData_IP(simulNUP107, simulNUP358, c('time','C_ip_density'), c('#A2C964', 'yellow'), 'NUP107_358_C_den.pdf' ) +#c('#A2C964', '#C1FF00') +# ---- +nc107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_NC, 'upper'=datain$NUP107_NC + datain$NUP107_NC_std, +'lower' = datain$NUP107_NC - datain$NUP107_NC_std) +c107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_C, 'upper'=datain$NUP107_C + datain$NUP107_C_std, +'lower' = datain$NUP107_C - datain$NUP107_C_std) +nc358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_NC, 'upper'=datain$NUP358_NC + datain$NUP358_NC_std, +'lower' = datain$NUP358_NC - datain$NUP358_NC_std) +c358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_C, 'upper'=datain$NUP358_C + datain$NUP358_C_std, +'lower' = datain$NUP358_C - datain$NUP358_C_std) +plotDataAll(nc107, c107, simulNUP107, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP107.pdf') +plotDataAll(nc358, c358, simulNUP358, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP358.pdf') +plotDataF(nc107, simulNUP107, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray80'), 0.4, 'NUP107_NC.pdf') +plotDataF(c107, simulNUP107, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray80'), 0.4, 'NUP107_C.pdf') +plotDataF(nc358, simulNUP358, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray80' ), 0.4, 'NUP358_NC.pdf') +plotDataF(c358, simulNUP358, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray80'), 0.4, 'NUP358_C.pdf') +plotData_IP(simulNUP107, simulNUP358, c('time','C_ip_density'), c('#A2C964', 'yellow'), 'NUP107_358_C_den.pdf' ) +#c('#A2C964', '#C1FF00') +# ---- +nc107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_NC, 'upper'=datain$NUP107_NC + datain$NUP107_NC_std, +'lower' = datain$NUP107_NC - datain$NUP107_NC_std) +c107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_C, 'upper'=datain$NUP107_C + datain$NUP107_C_std, +'lower' = datain$NUP107_C - datain$NUP107_C_std) +nc358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_NC, 'upper'=datain$NUP358_NC + datain$NUP358_NC_std, +'lower' = datain$NUP358_NC - datain$NUP358_NC_std) +c358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_C, 'upper'=datain$NUP358_C + datain$NUP358_C_std, +'lower' = datain$NUP358_C - datain$NUP358_C_std) +plotDataAll(nc107, c107, simulNUP107, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP107.pdf') +plotDataAll(nc358, c358, simulNUP358, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP358.pdf') +plotDataF(nc107, simulNUP107, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray20'), 0.4, 'NUP107_NC.pdf') +plotDataF(c107, simulNUP107, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray20'), 0.4, 'NUP107_C.pdf') +plotDataF(nc358, simulNUP358, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray20' ), 0.4, 'NUP358_NC.pdf') +plotDataF(c358, simulNUP358, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray20'), 0.4, 'NUP358_C.pdf') +plotData_IP(simulNUP107, simulNUP358, c('time','C_ip_density'), c('#A2C964', 'yellow'), 'NUP107_358_C_den.pdf' ) +#c('#A2C964', '#C1FF00') +# ----- +plotDataF <-function(d, simul, names, cols, alpha, fname) { +# function to plot data no separation between interphase and postmitotic +pdf(fname, width = 2.8, height=2, useDingbats = F,pointsize = 8, paper="special", fonts="Helvetica") +lw<- 2 +par(mar=c(4, 4, 1,1), mgp=c(2.5,1,0)) +plot(d$x, d$y, xlab = "", ylab = "Normalized total intensity", +type = "n", yaxt = 'n', xaxt = 'n', xlim = c(0, 120), ylim = c(0, 1.4), bty='l') +axis(1, tck = -0.05, at=c(0,20,40,60,80,100,120)) +axis(2, las=2, tck =-0.05, at=c(0,0.2,0.4,0.6,0.8,1,1.2)) +mtext(side=1, text="Time after AO (min)", line=2) +cols1_t <-adjustcolor(cols[1], alpha.f = alpha) +polygon(c(d$x, rev(d$x)), c(d$upper, rev(d$lower)), col = cols1_t, border = NA) +points(d$x, d$y, pch = 17, col = cols[1], cex=0.6) +colfun <- colorRampPalette(c(cols[1], cols[2])) +lines(simul[,names[1]], simul[,names[2]], lwd = lw, col = colfun(50)[20]) +lines(simul[,names[1]], simul[,names[3]], lwd = lw, lty = 2, col = cols[2]) +lines(simul[,names[1]], simul[,names[4]], lwd = lw, col=cols[2]) +par(xpd=TRUE) +legend(0,1.65, c("data", "post-mitotic+interphase", "post-mitotic", "interphase"), +lwd = c(0,1.5,1.5,1.5),pch = c(17,-1,-1,-1), lty = c(0,1,2,1), +col=c(cols[1], colfun(50)[20], cols[2], cols[2]), horiz = F, bty='n', +y.intersp = 0.8, cex = 0.8) +dev.off() +} +# ---- +nc107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_NC, 'upper'=datain$NUP107_NC + datain$NUP107_NC_std, +'lower' = datain$NUP107_NC - datain$NUP107_NC_std) +c107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_C, 'upper'=datain$NUP107_C + datain$NUP107_C_std, +'lower' = datain$NUP107_C - datain$NUP107_C_std) +nc358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_NC, 'upper'=datain$NUP358_NC + datain$NUP358_NC_std, +'lower' = datain$NUP358_NC - datain$NUP358_NC_std) +c358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_C, 'upper'=datain$NUP358_C + datain$NUP358_C_std, +'lower' = datain$NUP358_C - datain$NUP358_C_std) +plotDataAll(nc107, c107, simulNUP107, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP107.pdf') +plotDataAll(nc358, c358, simulNUP358, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP358.pdf') +plotDataF(nc107, simulNUP107, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray20'), 0.4, 'NUP107_NC.pdf') +plotDataF(c107, simulNUP107, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray20'), 0.4, 'NUP107_C.pdf') +plotDataF(nc358, simulNUP358, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray20' ), 0.4, 'NUP358_NC.pdf') +plotDataF(c358, simulNUP358, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray20'), 0.4, 'NUP358_C.pdf') +plotData_IP(simulNUP107, simulNUP358, c('time','C_ip_density'), c('#A2C964', 'yellow'), 'NUP107_358_C_den.pdf' ) +#c('#A2C964', '#C1FF00') +# ----- +plotDataF <-function(d, simul, names, cols, alpha, fname) { +# function to plot data no separation between interphase and postmitotic +pdf(fname, width = 2.8, height=2, useDingbats = F,pointsize = 8, paper="special", fonts="Helvetica") +lw<- 2 +par(mar=c(4, 4, 1,1), mgp=c(2.5,1,0)) +plot(d$x, d$y, xlab = "", ylab = "Normalized total intensity", +type = "n", yaxt = 'n', xaxt = 'n', xlim = c(0, 120), ylim = c(0, 1.4), bty='l') +axis(1, tck = -0.05, at=c(0,20,40,60,80,100,120)) +axis(2, las=2, tck =-0.05, at=c(0,0.2,0.4,0.6,0.8,1,1.2)) +mtext(side=1, text="Time after AO (min)", line=2) +cols1_t <-adjustcolor(cols[1], alpha.f = alpha) +polygon(c(d$x, rev(d$x)), c(d$upper, rev(d$lower)), col = cols1_t, border = NA) +points(d$x, d$y, pch = 17, col = cols[1], cex=0.6) +colfun <- colorRampPalette(c(cols[1], gray)) +lines(simul[,names[1]], simul[,names[2]], lwd = lw, col = colfun(50)[20]) +lines(simul[,names[1]], simul[,names[3]], lwd = lw, lty = 2, col = cols[2]) +lines(simul[,names[1]], simul[,names[4]], lwd = lw, col=cols[2]) +par(xpd=TRUE) +legend(0,1.65, c("data", "post-mitotic+interphase", "post-mitotic", "interphase"), +lwd = c(0,1.5,1.5,1.5),pch = c(17,-1,-1,-1), lty = c(0,1,2,1), +col=c(cols[1], colfun(50)[20], cols[2], cols[2]), horiz = F, bty='n', +y.intersp = 0.8, cex = 0.8) +dev.off() +} +# ---- +nc107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_NC, 'upper'=datain$NUP107_NC + datain$NUP107_NC_std, +'lower' = datain$NUP107_NC - datain$NUP107_NC_std) +c107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_C, 'upper'=datain$NUP107_C + datain$NUP107_C_std, +'lower' = datain$NUP107_C - datain$NUP107_C_std) +nc358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_NC, 'upper'=datain$NUP358_NC + datain$NUP358_NC_std, +'lower' = datain$NUP358_NC - datain$NUP358_NC_std) +c358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_C, 'upper'=datain$NUP358_C + datain$NUP358_C_std, +'lower' = datain$NUP358_C - datain$NUP358_C_std) +plotDataAll(nc107, c107, simulNUP107, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP107.pdf') +plotDataAll(nc358, c358, simulNUP358, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP358.pdf') +plotDataF(nc107, simulNUP107, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray20'), 0.4, 'NUP107_NC.pdf') +plotDataF(c107, simulNUP107, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray20'), 0.4, 'NUP107_C.pdf') +plotDataF(nc358, simulNUP358, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray20' ), 0.4, 'NUP358_NC.pdf') +plotDataF(c358, simulNUP358, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray20'), 0.4, 'NUP358_C.pdf') +plotData_IP(simulNUP107, simulNUP358, c('time','C_ip_density'), c('#A2C964', 'yellow'), 'NUP107_358_C_den.pdf' ) +#c('#A2C964', '#C1FF00') +# ----- +plotDataF <-function(d, simul, names, cols, alpha, fname) { +# function to plot data no separation between interphase and postmitotic +pdf(fname, width = 2.8, height=2, useDingbats = F,pointsize = 8, paper="special", fonts="Helvetica") +lw<- 2 +par(mar=c(4, 4, 1,1), mgp=c(2.5,1,0)) +plot(d$x, d$y, xlab = "", ylab = "Normalized total intensity", +type = "n", yaxt = 'n', xaxt = 'n', xlim = c(0, 120), ylim = c(0, 1.4), bty='l') +axis(1, tck = -0.05, at=c(0,20,40,60,80,100,120)) +axis(2, las=2, tck =-0.05, at=c(0,0.2,0.4,0.6,0.8,1,1.2)) +mtext(side=1, text="Time after AO (min)", line=2) +cols1_t <-adjustcolor(cols[1], alpha.f = alpha) +polygon(c(d$x, rev(d$x)), c(d$upper, rev(d$lower)), col = cols1_t, border = NA) +points(d$x, d$y, pch = 17, col = cols[1], cex=0.6) +colfun <- colorRampPalette(c(cols[1], 'gray')) +lines(simul[,names[1]], simul[,names[2]], lwd = lw, col = colfun(50)[20]) +lines(simul[,names[1]], simul[,names[3]], lwd = lw, lty = 2, col = cols[2]) +lines(simul[,names[1]], simul[,names[4]], lwd = lw, col=cols[2]) +par(xpd=TRUE) +legend(0,1.65, c("data", "post-mitotic+interphase", "post-mitotic", "interphase"), +lwd = c(0,1.5,1.5,1.5),pch = c(17,-1,-1,-1), lty = c(0,1,2,1), +col=c(cols[1], colfun(50)[20], cols[2], cols[2]), horiz = F, bty='n', +y.intersp = 0.8, cex = 0.8) +dev.off() +} +# ---- +nc107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_NC, 'upper'=datain$NUP107_NC + datain$NUP107_NC_std, +'lower' = datain$NUP107_NC - datain$NUP107_NC_std) +c107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_C, 'upper'=datain$NUP107_C + datain$NUP107_C_std, +'lower' = datain$NUP107_C - datain$NUP107_C_std) +nc358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_NC, 'upper'=datain$NUP358_NC + datain$NUP358_NC_std, +'lower' = datain$NUP358_NC - datain$NUP358_NC_std) +c358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_C, 'upper'=datain$NUP358_C + datain$NUP358_C_std, +'lower' = datain$NUP358_C - datain$NUP358_C_std) +plotDataAll(nc107, c107, simulNUP107, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP107.pdf') +plotDataAll(nc358, c358, simulNUP358, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP358.pdf') +plotDataF(nc107, simulNUP107, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray20'), 0.4, 'NUP107_NC.pdf') +plotDataF(c107, simulNUP107, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray20'), 0.4, 'NUP107_C.pdf') +plotDataF(nc358, simulNUP358, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray20' ), 0.4, 'NUP358_NC.pdf') +plotDataF(c358, simulNUP358, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray20'), 0.4, 'NUP358_C.pdf') +plotData_IP(simulNUP107, simulNUP358, c('time','C_ip_density'), c('#A2C964', 'yellow'), 'NUP107_358_C_den.pdf' ) +#c('#A2C964', '#C1FF00') +# ----- +plotDataF <-function(d, simul, names, cols, alpha, fname) { +# function to plot data no separation between interphase and postmitotic +pdf(fname, width = 2.8, height=2, useDingbats = F,pointsize = 8, paper="special", fonts="Helvetica") +lw<- 2 +par(mar=c(4, 4, 1,1), mgp=c(2.5,1,0)) +plot(d$x, d$y, xlab = "", ylab = "Normalized total intensity", +type = "n", yaxt = 'n', xaxt = 'n', xlim = c(0, 120), ylim = c(0, 1.4), bty='l') +axis(1, tck = -0.05, at=c(0,20,40,60,80,100,120)) +axis(2, las=2, tck =-0.05, at=c(0,0.2,0.4,0.6,0.8,1,1.2)) +mtext(side=1, text="Time after AO (min)", line=2) +cols1_t <-adjustcolor(cols[1], alpha.f = alpha) +polygon(c(d$x, rev(d$x)), c(d$upper, rev(d$lower)), col = cols1_t, border = NA) +points(d$x, d$y, pch = 17, col = cols[1], cex=0.6) +colfun <- colorRampPalette(c(cols[1], 'gray')) +lines(simul[,names[1]], simul[,names[2]], lwd = lw, col = colfun(50)[40]) +lines(simul[,names[1]], simul[,names[3]], lwd = lw, lty = 2, col = cols[2]) +lines(simul[,names[1]], simul[,names[4]], lwd = lw, col=cols[2]) +par(xpd=TRUE) +legend(0,1.65, c("data", "post-mitotic+interphase", "post-mitotic", "interphase"), +lwd = c(0,1.5,1.5,1.5),pch = c(17,-1,-1,-1), lty = c(0,1,2,1), +col=c(cols[1], colfun(50)[40], cols[2], cols[2]), horiz = F, bty='n', +y.intersp = 0.8, cex = 0.8) +dev.off() +} +# ---- +nc107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_NC, 'upper'=datain$NUP107_NC + datain$NUP107_NC_std, +'lower' = datain$NUP107_NC - datain$NUP107_NC_std) +c107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_C, 'upper'=datain$NUP107_C + datain$NUP107_C_std, +'lower' = datain$NUP107_C - datain$NUP107_C_std) +nc358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_NC, 'upper'=datain$NUP358_NC + datain$NUP358_NC_std, +'lower' = datain$NUP358_NC - datain$NUP358_NC_std) +c358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_C, 'upper'=datain$NUP358_C + datain$NUP358_C_std, +'lower' = datain$NUP358_C - datain$NUP358_C_std) +plotDataAll(nc107, c107, simulNUP107, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP107.pdf') +plotDataAll(nc358, c358, simulNUP358, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP358.pdf') +plotDataF(nc107, simulNUP107, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray20'), 0.4, 'NUP107_NC.pdf') +plotDataF(c107, simulNUP107, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray20'), 0.4, 'NUP107_C.pdf') +plotDataF(nc358, simulNUP358, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray20' ), 0.4, 'NUP358_NC.pdf') +plotDataF(c358, simulNUP358, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray20'), 0.4, 'NUP358_C.pdf') +plotData_IP(simulNUP107, simulNUP358, c('time','C_ip_density'), c('#A2C964', 'yellow'), 'NUP107_358_C_den.pdf' ) +#c('#A2C964', '#C1FF00') +# ----- +plotDataF <-function(d, simul, names, cols, alpha, fname) { +# function to plot data no separation between interphase and postmitotic +pdf(fname, width = 2.8, height=2, useDingbats = F,pointsize = 8, paper="special", fonts="Helvetica") +lw<- 2 +par(mar=c(4, 4, 1,1), mgp=c(2.5,1,0)) +plot(d$x, d$y, xlab = "", ylab = "Normalized total intensity", +type = "n", yaxt = 'n', xaxt = 'n', xlim = c(0, 120), ylim = c(0, 1.4), bty='l') +axis(1, tck = -0.05, at=c(0,20,40,60,80,100,120)) +axis(2, las=2, tck =-0.05, at=c(0,0.2,0.4,0.6,0.8,1,1.2)) +mtext(side=1, text="Time after AO (min)", line=2) +cols1_t <-adjustcolor(cols[1], alpha.f = alpha) +polygon(c(d$x, rev(d$x)), c(d$upper, rev(d$lower)), col = cols1_t, border = NA) +points(d$x, d$y, pch = 17, col = cols[1], cex=0.6) +colfun <- colorRampPalette(c(cols[1], 'gray60')) +lines(simul[,names[1]], simul[,names[2]], lwd = lw, col = colfun(50)[40]) +lines(simul[,names[1]], simul[,names[3]], lwd = lw, lty = 2, col = cols[2]) +lines(simul[,names[1]], simul[,names[4]], lwd = lw, col=cols[2]) +par(xpd=TRUE) +legend(0,1.65, c("data", "post-mitotic+interphase", "post-mitotic", "interphase"), +lwd = c(0,1.5,1.5,1.5),pch = c(17,-1,-1,-1), lty = c(0,1,2,1), +col=c(cols[1], colfun(50)[40], cols[2], cols[2]), horiz = F, bty='n', +y.intersp = 0.8, cex = 0.8) +dev.off() +} +# ---- +nc107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_NC, 'upper'=datain$NUP107_NC + datain$NUP107_NC_std, +'lower' = datain$NUP107_NC - datain$NUP107_NC_std) +c107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_C, 'upper'=datain$NUP107_C + datain$NUP107_C_std, +'lower' = datain$NUP107_C - datain$NUP107_C_std) +nc358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_NC, 'upper'=datain$NUP358_NC + datain$NUP358_NC_std, +'lower' = datain$NUP358_NC - datain$NUP358_NC_std) +c358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_C, 'upper'=datain$NUP358_C + datain$NUP358_C_std, +'lower' = datain$NUP358_C - datain$NUP358_C_std) +plotDataAll(nc107, c107, simulNUP107, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP107.pdf') +plotDataAll(nc358, c358, simulNUP358, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP358.pdf') +plotDataF(nc107, simulNUP107, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray20'), 0.4, 'NUP107_NC.pdf') +plotDataF(c107, simulNUP107, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray20'), 0.4, 'NUP107_C.pdf') +plotDataF(nc358, simulNUP358, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray20' ), 0.4, 'NUP358_NC.pdf') +plotDataF(c358, simulNUP358, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray20'), 0.4, 'NUP358_C.pdf') +plotData_IP(simulNUP107, simulNUP358, c('time','C_ip_density'), c('#A2C964', 'yellow'), 'NUP107_358_C_den.pdf' ) +#c('#A2C964', '#C1FF00') +8*0.8 +8*0.9 +8*0.85 +8*0.86 +8*0.87 +8*0.88 +# ----- +plotDataF <-function(d, simul, names, cols, alpha, fname) { +# function to plot data no separation between interphase and postmitotic +pdf(fname, width = 2.8, height=2, useDingbats = F,pointsize = 8, paper="special", fonts="Helvetica") +lw<- 2 +par(mar=c(4, 4, 1,1), mgp=c(2.5,1,0)) +plot(d$x, d$y, xlab = "", ylab = "Normalized total intensity", +type = "n", yaxt = 'n', xaxt = 'n', xlim = c(0, 120), ylim = c(0, 1.4), bty='l') +axis(1, tck = -0.05, at=c(0,20,40,60,80,100,120)) +axis(2, las=2, tck =-0.05, at=c(0,0.2,0.4,0.6,0.8,1,1.2)) +mtext(side=1, text="Time after AO (min)", line=2) +cols1_t <-adjustcolor(cols[1], alpha.f = alpha) +polygon(c(d$x, rev(d$x)), c(d$upper, rev(d$lower)), col = cols1_t, border = NA) +points(d$x, d$y, pch = 17, col = cols[1], cex=0.6) +colfun <- colorRampPalette(c(cols[1], 'gray60')) +lines(simul[,names[1]], simul[,names[2]], lwd = lw, col = colfun(50)[40]) +lines(simul[,names[1]], simul[,names[3]], lwd = lw, lty = 2, col = cols[2]) +lines(simul[,names[1]], simul[,names[4]], lwd = lw, col=cols[2]) +par(xpd=TRUE) +legend(0,1.65, c("data", "post-mitotic+interphase", "post-mitotic", "interphase"), +lwd = c(0,1.5,1.5,1.5),pch = c(17,-1,-1,-1), lty = c(0,1,2,1), +col=c(cols[1], colfun(50)[40], cols[2], cols[2]), horiz = F, bty='n', +y.intersp = 0.88, cex = 0.88) +dev.off() +} +# ----- +plotDataAll <-function(d1,d2, simul, names, cols, alpha, fname) { +# function to plot data no separation between interphase and postmitotic +pdf(fname, width = 2.8, height=2, useDingbats = F,pointsize = 8, paper="special", fonts="Helvetica") +lw<- 2 +par(mar=c(4, 4, 1,1), mgp=c(2.5,1,0)) +plot(d1$x, d1$y, xlab = "", ylab = "Normalized total intensity", +type = "n", yaxt = 'n', xaxt = 'n', xlim = c(0, 120), ylim = c(0, 1.4), bty='l') +#title('NUP107 inner-core region') +axis(1, tck = -0.05, at=c(0,20,40,60,80,100,120)) +axis(2, las=2, tck =-0.05, at=c(0,0.2,0.4,0.6,0.8,1,1.2)) +mtext(side=1, text="Time after AO (min)", line=2) +#transparent colors +cols1_t <-adjustcolor(cols[1], alpha.f = alpha) +cols2_t <-adjustcolor(cols[2], alpha.f = alpha) +#plot data for Non-core and core region +polygon(c(d1$x, rev(d1$x)), c(d1$upper, rev(d1$lower)), col = cols1_t, border = NA) +polygon(c(d1$x, rev(d1$x)), c(d2$upper, rev(d2$lower)), col = cols2_t, border = NA) +points(d1$x, d1$y, pch = 17, col = cols[1], cex=0.6) +points(d1$x, d2$y, pch = 17, col = cols[2], cex=0.6) +#plot simulations +lines(simul[,names[1]], simul[,names[2]], lwd = lw, col = cols[1]) +lines(simul[,names[1]], simul[,names[3]], lwd = lw, col = cols[2]) +par(xpd=TRUE) +legend(0,1.65, c("Non-core", "Inner-core","Non-core model", "Inner-core model"), +lwd = c(0,0, 1.5,1.5), pch = c(17,17, -1,-1), lty = c(0,0, 1,1), +col= c(cols[1], cols[2], cols[1], cols[2]), horiz = F, bty='n', +y.intersp = 0.88, cex = 0.88) +dev.off() +} +# ----- +plotData_IP <-function(simul1, simul2, names, cols, fname) { +pdf(fname, width = 2.8, height=2, useDingbats = F,pointsize = 8, paper="special", fonts="Helvetica") +lw<- 2 +par(mar=c(4, 4, 1,1), mgp=c(2.5,1,0)) +plot(simul1[,names[1]], simul1[,names[2]], xlab = "", ylab = "Normalized interphase density", +type = "n", yaxt = 'n', xaxt = 'n', xlim = c(0, 120), ylim = c(0, 1.1), bty='l') +axis(1, tck = -0.05, at=c(0,20,40,60,80,100,120)) +axis(2, las=2, tck =-0.05, at=c(0,0.2,0.4,0.6,0.8,1)) +mtext(side=1, text="Time after AO (min)", line=2) +lines(simul1[,names[1]], simul1[,names[2]], lwd = lw, lty = 1, col = cols[1]) +lines(simul2[,names[1]], simul2[,names[2]], lwd = lw, lty = 1, col = cols[2]) +par(xpd=TRUE) +legend(0,1.1, c("Nup107", "Nup358"), +lwd = c(1.5,1.5), +col=c(cols[1], cols[2]), horiz = F, bty='n', +y.intersp = 0.88, cex = 0.88) +dev.off() +} +# ---- +nc107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_NC, 'upper'=datain$NUP107_NC + datain$NUP107_NC_std, +'lower' = datain$NUP107_NC - datain$NUP107_NC_std) +c107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_C, 'upper'=datain$NUP107_C + datain$NUP107_C_std, +'lower' = datain$NUP107_C - datain$NUP107_C_std) +nc358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_NC, 'upper'=datain$NUP358_NC + datain$NUP358_NC_std, +'lower' = datain$NUP358_NC - datain$NUP358_NC_std) +c358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_C, 'upper'=datain$NUP358_C + datain$NUP358_C_std, +'lower' = datain$NUP358_C - datain$NUP358_C_std) +plotDataAll(nc107, c107, simulNUP107, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP107.pdf') +plotDataAll(nc358, c358, simulNUP358, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP358.pdf') +plotDataF(nc107, simulNUP107, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray20'), 0.4, 'NUP107_NC.pdf') +plotDataF(c107, simulNUP107, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray20'), 0.4, 'NUP107_C.pdf') +plotDataF(nc358, simulNUP358, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray20' ), 0.4, 'NUP358_NC.pdf') +plotDataF(c358, simulNUP358, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray20'), 0.4, 'NUP358_C.pdf') +plotData_IP(simulNUP107, simulNUP358, c('time','C_ip_density'), c('#A2C964', 'yellow'), 'NUP107_358_C_den.pdf' ) +#c('#A2C964', '#C1FF00') +# ---- +nc107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_NC, 'upper'=datain$NUP107_NC + datain$NUP107_NC_std, +'lower' = datain$NUP107_NC - datain$NUP107_NC_std) +c107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_C, 'upper'=datain$NUP107_C + datain$NUP107_C_std, +'lower' = datain$NUP107_C - datain$NUP107_C_std) +nc358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_NC, 'upper'=datain$NUP358_NC + datain$NUP358_NC_std, +'lower' = datain$NUP358_NC - datain$NUP358_NC_std) +c358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_C, 'upper'=datain$NUP358_C + datain$NUP358_C_std, +'lower' = datain$NUP358_C - datain$NUP358_C_std) +plotDataAll(nc107, c107, simulNUP107, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP107.pdf') +plotDataAll(nc358, c358, simulNUP358, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP358.pdf') +plotDataF(nc107, simulNUP107, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray20'), 0.4, 'NUP107_NC.pdf') +plotDataF(c107, simulNUP107, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray20'), 0.4, 'NUP107_C.pdf') +plotDataF(nc358, simulNUP358, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray20' ), 0.4, 'NUP358_NC.pdf') +plotDataF(c358, simulNUP358, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray20'), 0.4, 'NUP358_C.pdf') +plotData_IP(simulNUP107, simulNUP358, c('time','C_ip_density'), c('red', 'green'), 'NUP107_358_C_den.pdf' ) +#c('#A2C964', '#C1FF00') +# ---- +nc107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_NC, 'upper'=datain$NUP107_NC + datain$NUP107_NC_std, +'lower' = datain$NUP107_NC - datain$NUP107_NC_std) +c107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_C, 'upper'=datain$NUP107_C + datain$NUP107_C_std, +'lower' = datain$NUP107_C - datain$NUP107_C_std) +nc358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_NC, 'upper'=datain$NUP358_NC + datain$NUP358_NC_std, +'lower' = datain$NUP358_NC - datain$NUP358_NC_std) +c358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_C, 'upper'=datain$NUP358_C + datain$NUP358_C_std, +'lower' = datain$NUP358_C - datain$NUP358_C_std) +plotDataAll(nc107, c107, simulNUP107, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP107.pdf') +plotDataAll(nc358, c358, simulNUP358, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP358.pdf') +plotDataF(nc107, simulNUP107, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray20'), 0.4, 'NUP107_NC.pdf') +plotDataF(c107, simulNUP107, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray20'), 0.4, 'NUP107_C.pdf') +plotDataF(nc358, simulNUP358, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray20' ), 0.4, 'NUP358_NC.pdf') +plotDataF(c358, simulNUP358, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray20'), 0.4, 'NUP358_C.pdf') +plotData_IP(simulNUP107, simulNUP358, c('time','C_ip_density'), c('#A2C964' , '#d01c8b'), 'NUP107_358_C_den.pdf' ) +#c('#A2C964', '#C1FF00') diff --git a/results/fitnup107nup358/NUP107.pdf b/results/fitnup107nup358/NUP107.pdf new file mode 100644 index 0000000000000000000000000000000000000000..89a0f9378b2f944985d8a8cf8546d42477ef78d3 Binary files /dev/null and b/results/fitnup107nup358/NUP107.pdf differ diff --git a/results/fitnup107nup358/NUP107_358_C_den.pdf b/results/fitnup107nup358/NUP107_358_C_den.pdf new file mode 100644 index 0000000000000000000000000000000000000000..04cd9f74d76630400ff4900a1dcdeef7cb5b1afe Binary files /dev/null and b/results/fitnup107nup358/NUP107_358_C_den.pdf differ diff --git a/results/fitnup107nup358/NUP107_C.pdf b/results/fitnup107nup358/NUP107_C.pdf new file mode 100644 index 0000000000000000000000000000000000000000..2483394d7f5a0cfb8265d03252161c92c5e95ec5 Binary files /dev/null and b/results/fitnup107nup358/NUP107_C.pdf differ diff --git a/results/fitnup107nup358/NUP107_NC.pdf b/results/fitnup107nup358/NUP107_NC.pdf new file mode 100644 index 0000000000000000000000000000000000000000..4395629aa7bb4c6f6b2057a616112cad8a29ffa0 Binary files /dev/null and b/results/fitnup107nup358/NUP107_NC.pdf differ diff --git a/results/fitnup107nup358/NUP358.pdf b/results/fitnup107nup358/NUP358.pdf new file mode 100644 index 0000000000000000000000000000000000000000..650e49b289864f0c091b6ed9df7f96e82111ae8e Binary files /dev/null and b/results/fitnup107nup358/NUP358.pdf differ diff --git a/results/fitnup107nup358/NUP358_C.pdf b/results/fitnup107nup358/NUP358_C.pdf new file mode 100644 index 0000000000000000000000000000000000000000..a0e8c76269d27df67e19539512a3d4e928724df5 Binary files /dev/null and b/results/fitnup107nup358/NUP358_C.pdf differ diff --git a/results/fitnup107nup358/NUP358_NC.pdf b/results/fitnup107nup358/NUP358_NC.pdf new file mode 100644 index 0000000000000000000000000000000000000000..97053e34dda762ba724bdddff78c0401e0d92605 Binary files /dev/null and b/results/fitnup107nup358/NUP358_NC.pdf differ diff --git a/results/fitnup107nup358/evaluation_of_fits b/results/fitnup107nup358/evaluation_of_fits new file mode 100644 index 0000000000000000000000000000000000000000..55a38c0591b9c4a1e8a3e7e95f0438541024ce21 Binary files /dev/null and b/results/fitnup107nup358/evaluation_of_fits differ diff --git a/results/fitnup107nup358/k1_ip_CI.txt b/results/fitnup107nup358/k1_ip_CI.txt new file mode 100644 index 0000000000000000000000000000000000000000..0804f3b53da0bc2984c0f81dc9c7568dd8864bb3 --- /dev/null +++ b/results/fitnup107nup358/k1_ip_CI.txt @@ -0,0 +1,101 @@ +%k1_pm k2_pm k1_ip k2_ip norm normBest +3.528e-01 4.546e-02 2.797e-02 8.925e-02 1.301e+02 9.449e+01 +3.529e-01 4.545e-02 2.823e-02 8.760e-02 1.284e+02 9.449e+01 +3.530e-01 4.544e-02 2.849e-02 8.612e-02 1.267e+02 9.449e+01 +3.531e-01 4.543e-02 2.876e-02 8.469e-02 1.251e+02 9.449e+01 +3.533e-01 4.541e-02 2.902e-02 8.354e-02 1.236e+02 9.449e+01 +3.533e-01 4.540e-02 2.929e-02 8.205e-02 1.221e+02 9.449e+01 +3.534e-01 4.540e-02 2.956e-02 8.057e-02 1.207e+02 9.449e+01 +3.536e-01 4.538e-02 2.984e-02 7.958e-02 1.193e+02 9.449e+01 +3.536e-01 4.537e-02 3.011e-02 7.824e-02 1.179e+02 9.449e+01 +3.537e-01 4.536e-02 3.039e-02 7.708e-02 1.166e+02 9.449e+01 +3.537e-01 4.535e-02 3.067e-02 7.584e-02 1.154e+02 9.449e+01 +3.538e-01 4.534e-02 3.096e-02 7.474e-02 1.142e+02 9.449e+01 +3.539e-01 4.533e-02 3.124e-02 7.367e-02 1.130e+02 9.449e+01 +3.539e-01 4.533e-02 3.153e-02 7.258e-02 1.119e+02 9.449e+01 +3.540e-01 4.532e-02 3.182e-02 7.158e-02 1.108e+02 9.449e+01 +3.541e-01 4.531e-02 3.212e-02 7.061e-02 1.098e+02 9.449e+01 +3.541e-01 4.530e-02 3.241e-02 6.967e-02 1.088e+02 9.449e+01 +3.542e-01 4.529e-02 3.271e-02 6.875e-02 1.078e+02 9.449e+01 +3.543e-01 4.528e-02 3.302e-02 6.786e-02 1.069e+02 9.449e+01 +3.543e-01 4.527e-02 3.332e-02 6.700e-02 1.060e+02 9.449e+01 +3.544e-01 4.526e-02 3.363e-02 6.616e-02 1.052e+02 9.449e+01 +3.544e-01 4.525e-02 3.394e-02 6.534e-02 1.044e+02 9.449e+01 +3.545e-01 4.524e-02 3.426e-02 6.455e-02 1.036e+02 9.449e+01 +3.545e-01 4.523e-02 3.457e-02 6.377e-02 1.029e+02 9.449e+01 +3.546e-01 4.522e-02 3.489e-02 6.302e-02 1.022e+02 9.449e+01 +3.546e-01 4.521e-02 3.522e-02 6.228e-02 1.015e+02 9.449e+01 +3.546e-01 4.520e-02 3.554e-02 6.157e-02 1.009e+02 9.449e+01 +3.547e-01 4.519e-02 3.587e-02 6.087e-02 1.003e+02 9.449e+01 +3.547e-01 4.519e-02 3.620e-02 6.019e-02 9.976e+01 9.449e+01 +3.548e-01 4.517e-02 3.654e-02 5.963e-02 9.924e+01 9.449e+01 +3.548e-01 4.516e-02 3.688e-02 5.898e-02 9.874e+01 9.449e+01 +3.548e-01 4.515e-02 3.722e-02 5.833e-02 9.828e+01 9.449e+01 +3.549e-01 4.515e-02 3.756e-02 5.771e-02 9.784e+01 9.449e+01 +3.549e-01 4.514e-02 3.791e-02 5.710e-02 9.744e+01 9.449e+01 +3.549e-01 4.513e-02 3.826e-02 5.651e-02 9.706e+01 9.449e+01 +3.549e-01 4.512e-02 3.861e-02 5.593e-02 9.672e+01 9.449e+01 +3.549e-01 4.512e-02 3.897e-02 5.536e-02 9.640e+01 9.449e+01 +3.549e-01 4.512e-02 3.933e-02 5.474e-02 9.610e+01 9.449e+01 +3.549e-01 4.511e-02 3.970e-02 5.420e-02 9.583e+01 9.449e+01 +3.549e-01 4.510e-02 4.006e-02 5.369e-02 9.559e+01 9.449e+01 +3.549e-01 4.509e-02 4.043e-02 5.318e-02 9.538e+01 9.449e+01 +3.549e-01 4.509e-02 4.081e-02 5.268e-02 9.519e+01 9.449e+01 +3.550e-01 4.507e-02 4.119e-02 5.226e-02 9.502e+01 9.449e+01 +3.550e-01 4.506e-02 4.157e-02 5.178e-02 9.488e+01 9.449e+01 +3.550e-01 4.506e-02 4.195e-02 5.131e-02 9.476e+01 9.449e+01 +3.550e-01 4.505e-02 4.234e-02 5.084e-02 9.466e+01 9.449e+01 +3.549e-01 4.505e-02 4.273e-02 5.035e-02 9.459e+01 9.449e+01 +3.550e-01 4.504e-02 4.313e-02 4.998e-02 9.453e+01 9.449e+01 +3.549e-01 4.504e-02 4.353e-02 4.948e-02 9.450e+01 9.449e+01 +3.549e-01 4.503e-02 4.393e-02 4.910e-02 9.449e+01 9.449e+01 +3.549e-01 4.502e-02 4.433e-02 4.872e-02 9.450e+01 9.449e+01 +3.549e-01 4.502e-02 4.474e-02 4.826e-02 9.453e+01 9.449e+01 +3.549e-01 4.500e-02 4.516e-02 4.792e-02 9.458e+01 9.449e+01 +3.549e-01 4.500e-02 4.558e-02 4.754e-02 9.465e+01 9.449e+01 +3.548e-01 4.499e-02 4.600e-02 4.717e-02 9.474e+01 9.449e+01 +3.548e-01 4.499e-02 4.642e-02 4.676e-02 9.485e+01 9.449e+01 +3.548e-01 4.498e-02 4.685e-02 4.640e-02 9.497e+01 9.449e+01 +3.547e-01 4.498e-02 4.729e-02 4.605e-02 9.511e+01 9.449e+01 +3.547e-01 4.497e-02 4.772e-02 4.571e-02 9.527e+01 9.449e+01 +3.547e-01 4.496e-02 4.817e-02 4.537e-02 9.544e+01 9.449e+01 +3.546e-01 4.496e-02 4.861e-02 4.504e-02 9.564e+01 9.449e+01 +3.546e-01 4.495e-02 4.906e-02 4.472e-02 9.584e+01 9.449e+01 +3.546e-01 4.494e-02 4.952e-02 4.440e-02 9.607e+01 9.449e+01 +3.545e-01 4.494e-02 4.997e-02 4.409e-02 9.631e+01 9.449e+01 +3.545e-01 4.493e-02 5.044e-02 4.379e-02 9.656e+01 9.449e+01 +3.545e-01 4.493e-02 5.090e-02 4.349e-02 9.683e+01 9.449e+01 +3.544e-01 4.492e-02 5.137e-02 4.320e-02 9.711e+01 9.449e+01 +3.544e-01 4.491e-02 5.185e-02 4.291e-02 9.741e+01 9.449e+01 +3.544e-01 4.491e-02 5.233e-02 4.263e-02 9.772e+01 9.449e+01 +3.543e-01 4.490e-02 5.281e-02 4.236e-02 9.805e+01 9.449e+01 +3.543e-01 4.489e-02 5.330e-02 4.209e-02 9.838e+01 9.449e+01 +3.542e-01 4.489e-02 5.379e-02 4.182e-02 9.873e+01 9.449e+01 +3.542e-01 4.488e-02 5.429e-02 4.156e-02 9.910e+01 9.449e+01 +3.541e-01 4.488e-02 5.479e-02 4.131e-02 9.947e+01 9.449e+01 +3.541e-01 4.487e-02 5.530e-02 4.106e-02 9.986e+01 9.449e+01 +3.540e-01 4.486e-02 5.581e-02 4.082e-02 1.003e+02 9.449e+01 +3.540e-01 4.486e-02 5.633e-02 4.058e-02 1.007e+02 9.449e+01 +3.539e-01 4.485e-02 5.685e-02 4.034e-02 1.011e+02 9.449e+01 +3.539e-01 4.484e-02 5.738e-02 4.011e-02 1.015e+02 9.449e+01 +3.538e-01 4.484e-02 5.791e-02 3.988e-02 1.020e+02 9.449e+01 +3.537e-01 4.483e-02 5.844e-02 3.966e-02 1.024e+02 9.449e+01 +3.537e-01 4.483e-02 5.898e-02 3.944e-02 1.029e+02 9.449e+01 +3.536e-01 4.482e-02 5.953e-02 3.923e-02 1.034e+02 9.449e+01 +3.536e-01 4.481e-02 6.008e-02 3.902e-02 1.038e+02 9.449e+01 +3.535e-01 4.481e-02 6.064e-02 3.881e-02 1.043e+02 9.449e+01 +3.534e-01 4.480e-02 6.120e-02 3.861e-02 1.048e+02 9.449e+01 +3.534e-01 4.479e-02 6.176e-02 3.841e-02 1.053e+02 9.449e+01 +3.533e-01 4.479e-02 6.234e-02 3.821e-02 1.059e+02 9.449e+01 +3.532e-01 4.478e-02 6.291e-02 3.802e-02 1.064e+02 9.449e+01 +3.532e-01 4.478e-02 6.349e-02 3.783e-02 1.069e+02 9.449e+01 +3.531e-01 4.477e-02 6.408e-02 3.765e-02 1.075e+02 9.449e+01 +3.530e-01 4.476e-02 6.468e-02 3.746e-02 1.080e+02 9.449e+01 +3.530e-01 4.476e-02 6.527e-02 3.728e-02 1.086e+02 9.449e+01 +3.529e-01 4.475e-02 6.588e-02 3.711e-02 1.092e+02 9.449e+01 +3.528e-01 4.475e-02 6.649e-02 3.693e-02 1.097e+02 9.449e+01 +3.527e-01 4.474e-02 6.710e-02 3.676e-02 1.103e+02 9.449e+01 +3.527e-01 4.474e-02 6.772e-02 3.660e-02 1.109e+02 9.449e+01 +3.526e-01 4.473e-02 6.835e-02 3.643e-02 1.115e+02 9.449e+01 +3.525e-01 4.472e-02 6.898e-02 3.627e-02 1.121e+02 9.449e+01 +3.524e-01 4.472e-02 6.962e-02 3.611e-02 1.127e+02 9.449e+01 diff --git a/results/fitnup107nup358/k1_ip_CI_coeffMax.txt b/results/fitnup107nup358/k1_ip_CI_coeffMax.txt new file mode 100644 index 0000000000000000000000000000000000000000..8f9931c9e3c79f5c60345c6051c97bb70663ca7d --- /dev/null +++ b/results/fitnup107nup358/k1_ip_CI_coeffMax.txt @@ -0,0 +1,101 @@ +%k1_pm k2_pm k1_ip k2_ip norm normBest +0.3433 0.04389 0.02379 0.03893 1.012 1.025 1.068 1.13 86.03 78.32 +0.3436 0.04389 0.02401 0.03869 1.012 1.025 1.067 1.129 85.76 78.32 +0.3439 0.04389 0.02423 0.03844 1.011 1.025 1.066 1.128 85.49 78.32 +0.3442 0.04389 0.02446 0.03819 1.011 1.025 1.064 1.127 85.22 78.32 NUP107 Nup358 +0.3445 0.04389 0.02468 0.03795 1.011 1.025 1.063 1.126 84.95 78.32 Postmitotic 0.354 0.0437 +0.3448 0.04389 0.02491 0.0377 1.011 1.024 1.062 1.125 84.69 78.32 Interphase 0.0374 0.0276 +0.3451 0.04389 0.02514 0.03745 1.011 1.024 1.06 1.125 84.43 78.32 +0.3454 0.04389 0.02537 0.0372 1.011 1.024 1.059 1.124 84.18 78.32 +0.3457 0.04388 0.02561 0.03702 1.01 1.024 1.058 1.123 83.93 78.32 Postmitotic 0.023 0.0041 +0.346 0.04388 0.02585 0.03671 1.01 1.024 1.056 1.122 83.68 78.32 Interphase 0.0012 0.0072 +0.3462 0.04388 0.02608 0.0364 1.01 1.024 1.055 1.122 83.44 78.32 +0.3465 0.04388 0.02633 0.03616 1.01 1.024 1.054 1.121 83.2 78.32 +0.3468 0.04388 0.02657 0.03591 1.01 1.024 1.053 1.12 82.97 78.32 +0.347 0.04388 0.02682 0.03567 1.009 1.024 1.052 1.12 82.74 78.32 +0.3473 0.04387 0.02706 0.03543 1.009 1.024 1.05 1.119 82.52 78.32 +0.3476 0.04387 0.02731 0.03538 1.009 1.023 1.049 1.118 82.3 78.32 0.0229 +0.3479 0.04387 0.02757 0.03513 1.009 1.023 1.048 1.117 82.08 78.32 0.0012 +0.3482 0.04386 0.02782 0.03488 1.009 1.023 1.047 1.117 81.87 78.32 +0.3484 0.04386 0.02808 0.03463 1.009 1.023 1.046 1.116 81.67 78.32 +0.3487 0.04386 0.02834 0.03438 1.008 1.023 1.045 1.116 81.47 78.32 +0.3489 0.04385 0.0286 0.03413 1.008 1.023 1.044 1.115 81.27 78.32 +0.3491 0.04385 0.02887 0.03389 1.008 1.023 1.042 1.115 81.08 78.32 +0.3493 0.04385 0.02913 0.03349 1.008 1.023 1.041 1.115 80.9 78.32 +0.3496 0.04384 0.0294 0.03326 1.008 1.023 1.04 1.114 80.72 78.32 +0.3498 0.04384 0.02967 0.03303 1.008 1.023 1.039 1.114 80.55 78.32 +0.3501 0.04383 0.02995 0.03296 1.007 1.023 1.038 1.113 80.38 78.32 +0.3503 0.04383 0.03023 0.03271 1.007 1.023 1.037 1.112 80.22 78.32 +0.3505 0.04382 0.03051 0.03246 1.007 1.023 1.036 1.112 80.06 78.32 +0.3507 0.04382 0.03079 0.03209 1.007 1.023 1.035 1.112 79.91 78.32 +0.3509 0.04381 0.03107 0.03187 1.007 1.023 1.034 1.112 79.77 78.32 +0.3512 0.0438 0.03136 0.03178 1.007 1.023 1.033 1.111 79.63 78.32 +0.3514 0.0438 0.03165 0.03153 1.007 1.023 1.032 1.111 79.5 78.32 +0.3516 0.04379 0.03194 0.03119 1.006 1.023 1.031 1.111 79.38 78.32 +0.3518 0.04379 0.03224 0.03098 1.006 1.023 1.03 1.111 79.26 78.32 +0.352 0.04378 0.03254 0.03078 1.006 1.023 1.029 1.111 79.15 78.32 +0.3522 0.04377 0.03284 0.03057 1.006 1.023 1.029 1.11 79.05 78.32 +0.3524 0.04376 0.03314 0.03037 1.006 1.022 1.028 1.11 78.95 78.32 +0.3526 0.04376 0.03345 0.03016 1.006 1.022 1.027 1.11 78.86 78.32 +0.3528 0.04375 0.03376 0.02996 1.006 1.022 1.026 1.11 78.77 78.32 +0.353 0.04374 0.03407 0.02975 1.006 1.022 1.025 1.109 78.7 78.32 +0.3531 0.04373 0.03439 0.02942 1.005 1.023 1.024 1.11 78.63 78.32 +0.3533 0.04373 0.0347 0.0292 1.005 1.023 1.023 1.11 78.56 78.32 +0.3535 0.04372 0.03503 0.02898 1.005 1.023 1.022 1.11 78.51 78.32 +0.3537 0.04371 0.03535 0.02876 1.005 1.023 1.022 1.11 78.46 78.32 +0.3538 0.0437 0.03568 0.02855 1.005 1.023 1.021 1.11 78.42 78.32 +0.354 0.04369 0.03601 0.02835 1.005 1.023 1.02 1.11 78.38 78.32 +0.3542 0.04368 0.03634 0.02814 1.005 1.023 1.019 1.11 78.36 78.32 +0.3544 0.04367 0.03668 0.02807 1.005 1.023 1.019 1.109 78.34 78.32 +0.3545 0.04366 0.03702 0.02773 1.004 1.023 1.018 1.11 78.33 78.32 +0.3547 0.04365 0.03736 0.02759 1.004 1.023 1.017 1.109 78.32 78.32 +0.3549 0.04364 0.0377 0.02745 1.004 1.023 1.016 1.109 78.33 78.32 +0.355 0.04363 0.03805 0.02724 1.004 1.023 1.016 1.109 78.34 78.32 +0.3552 0.04362 0.03841 0.02701 1.004 1.023 1.015 1.11 78.36 78.32 +0.3553 0.04361 0.03876 0.02683 1.004 1.023 1.014 1.11 78.38 78.32 +0.3555 0.0436 0.03912 0.02668 1.004 1.023 1.013 1.109 78.42 78.32 +0.3557 0.04359 0.03948 0.02648 1.004 1.023 1.013 1.11 78.46 78.32 +0.3557 0.04358 0.03985 0.02617 1.004 1.023 1.012 1.11 78.51 78.32 +0.3559 0.04357 0.04022 0.0261 1.003 1.023 1.011 1.11 78.57 78.32 +0.3561 0.04355 0.04059 0.02592 1.003 1.023 1.011 1.11 78.64 78.32 +0.3562 0.04354 0.04096 0.02574 1.003 1.023 1.01 1.11 78.71 78.32 +0.3563 0.04353 0.04134 0.02556 1.003 1.023 1.009 1.11 78.79 78.32 +0.3565 0.04352 0.04172 0.02534 1.003 1.023 1.009 1.111 78.88 78.32 +0.3566 0.04351 0.04211 0.02516 1.003 1.023 1.008 1.111 78.98 78.32 +0.3567 0.04349 0.0425 0.02498 1.003 1.023 1.008 1.111 79.08 78.32 +0.3568 0.04348 0.04289 0.02481 1.003 1.023 1.007 1.112 79.2 78.32 +0.357 0.04347 0.04329 0.02464 1.003 1.023 1.006 1.112 79.32 78.32 +0.3571 0.04345 0.04369 0.02447 1.003 1.023 1.006 1.112 79.45 78.32 +0.3572 0.04344 0.04409 0.0243 1.003 1.023 1.005 1.112 79.59 78.32 +0.3573 0.04343 0.0445 0.02413 1.002 1.024 1.005 1.113 79.73 78.32 +0.3574 0.04341 0.04491 0.02396 1.002 1.024 1.004 1.113 79.88 78.32 +0.3575 0.0434 0.04533 0.0238 1.002 1.024 1.004 1.113 80.05 78.32 +0.3576 0.04339 0.04575 0.0236 1.002 1.024 1.003 1.114 80.21 78.32 +0.3577 0.04337 0.04617 0.02343 1.002 1.024 1.003 1.114 80.39 78.32 +0.3578 0.04336 0.0466 0.02327 1.002 1.024 1.002 1.115 80.58 78.32 +0.3579 0.04334 0.04703 0.02312 1.002 1.024 1.002 1.115 80.77 78.32 +0.358 0.04333 0.04747 0.02296 1.002 1.024 1.001 1.115 80.97 78.32 +0.3581 0.04331 0.04791 0.02281 1.002 1.024 1 1.116 81.18 78.32 +0.3582 0.0433 0.04835 0.02266 1.002 1.024 1 1.116 81.39 78.32 +0.3583 0.04329 0.0488 0.02256 1.002 1.024 0.9995 1.116 81.62 78.32 +0.3584 0.04327 0.04925 0.02241 1.002 1.024 0.9991 1.117 81.85 78.32 +0.3585 0.04326 0.0497 0.02227 1.002 1.025 0.9986 1.117 82.08 78.32 +0.3586 0.04324 0.05016 0.02213 1.001 1.025 0.9981 1.117 82.33 78.32 +0.3587 0.04323 0.05063 0.02199 1.001 1.025 0.9977 1.118 82.58 78.32 +0.3588 0.04321 0.0511 0.02185 1.001 1.025 0.9972 1.118 82.85 78.32 +0.3588 0.0432 0.05157 0.02172 1.001 1.025 0.9968 1.118 83.11 78.32 +0.3589 0.04318 0.05205 0.02159 1.001 1.025 0.9964 1.119 83.39 78.32 +0.359 0.04316 0.05253 0.02146 1.001 1.025 0.9959 1.119 83.67 78.32 +0.3591 0.04315 0.05301 0.02133 1.001 1.025 0.9955 1.12 83.96 78.32 +0.3592 0.04313 0.0535 0.02121 1.001 1.025 0.9951 1.12 84.26 78.32 +0.3592 0.04312 0.054 0.02108 1.001 1.025 0.9947 1.12 84.56 78.32 +0.3593 0.0431 0.0545 0.02096 1.001 1.025 0.9943 1.121 84.87 78.32 +0.3594 0.04309 0.055 0.02084 1.001 1.025 0.9939 1.121 85.19 78.32 +0.3595 0.04307 0.05551 0.02072 1.001 1.026 0.9935 1.122 85.52 78.32 +0.3595 0.04305 0.05603 0.0205 1.001 1.026 0.9931 1.123 85.85 78.32 +0.3596 0.04304 0.05654 0.02043 1.001 1.026 0.9927 1.123 86.19 78.32 +0.3596 0.04302 0.05707 0.02031 1.001 1.026 0.9924 1.123 86.53 78.32 +0.3597 0.04301 0.0576 0.02019 1.001 1.026 0.992 1.124 86.88 78.32 +0.3597 0.04299 0.05813 0.02007 1.001 1.026 0.9916 1.124 87.24 78.32 +0.3598 0.04297 0.05867 0.01996 1 1.026 0.9913 1.125 87.61 78.32 +0.3598 0.04296 0.05921 0.01985 1 1.026 0.9909 1.125 87.98 78.32 diff --git a/results/fitnup107nup358/k1_pm_CI.txt b/results/fitnup107nup358/k1_pm_CI.txt new file mode 100644 index 0000000000000000000000000000000000000000..39d657d5b830238a4e0e90b1bf9ced94f2fd856f --- /dev/null +++ b/results/fitnup107nup358/k1_pm_CI.txt @@ -0,0 +1,101 @@ +%k1_pm k2_pm k1_ip k2_ip norm normBest +2.260e-01 5.010e-02 4.450e-02 4.540e-02 1.256e+02 9.449e+01 +2.281e-01 4.996e-02 4.447e-02 4.550e-02 1.243e+02 9.449e+01 +2.302e-01 4.983e-02 4.444e-02 4.558e-02 1.230e+02 9.449e+01 +2.323e-01 4.969e-02 4.442e-02 4.567e-02 1.217e+02 9.449e+01 +2.345e-01 4.956e-02 4.439e-02 4.576e-02 1.205e+02 9.449e+01 +2.367e-01 4.943e-02 4.437e-02 4.585e-02 1.193e+02 9.449e+01 +2.388e-01 4.930e-02 4.435e-02 4.594e-02 1.181e+02 9.449e+01 +2.411e-01 4.918e-02 4.433e-02 4.602e-02 1.170e+02 9.449e+01 +2.433e-01 4.905e-02 4.431e-02 4.611e-02 1.159e+02 9.449e+01 +2.455e-01 4.893e-02 4.429e-02 4.620e-02 1.148e+02 9.449e+01 +2.478e-01 4.881e-02 4.427e-02 4.628e-02 1.137e+02 9.449e+01 +2.501e-01 4.869e-02 4.425e-02 4.636e-02 1.127e+02 9.449e+01 +2.524e-01 4.857e-02 4.423e-02 4.645e-02 1.117e+02 9.449e+01 +2.547e-01 4.845e-02 4.421e-02 4.653e-02 1.108e+02 9.449e+01 +2.571e-01 4.834e-02 4.419e-02 4.661e-02 1.098e+02 9.449e+01 +2.595e-01 4.822e-02 4.418e-02 4.669e-02 1.089e+02 9.449e+01 +2.619e-01 4.811e-02 4.416e-02 4.677e-02 1.081e+02 9.449e+01 +2.643e-01 4.800e-02 4.415e-02 4.685e-02 1.072e+02 9.449e+01 +2.668e-01 4.789e-02 4.413e-02 4.693e-02 1.064e+02 9.449e+01 +2.692e-01 4.778e-02 4.412e-02 4.701e-02 1.056e+02 9.449e+01 +2.717e-01 4.767e-02 4.411e-02 4.709e-02 1.049e+02 9.449e+01 +2.742e-01 4.756e-02 4.409e-02 4.717e-02 1.041e+02 9.449e+01 +2.768e-01 4.746e-02 4.408e-02 4.724e-02 1.034e+02 9.449e+01 +2.793e-01 4.735e-02 4.407e-02 4.732e-02 1.028e+02 9.449e+01 +2.819e-01 4.725e-02 4.406e-02 4.740e-02 1.021e+02 9.449e+01 +2.845e-01 4.715e-02 4.405e-02 4.747e-02 1.015e+02 9.449e+01 +2.872e-01 4.705e-02 4.404e-02 4.755e-02 1.009e+02 9.449e+01 +2.898e-01 4.695e-02 4.403e-02 4.762e-02 1.004e+02 9.449e+01 +2.925e-01 4.685e-02 4.402e-02 4.769e-02 9.982e+01 9.449e+01 +2.952e-01 4.675e-02 4.401e-02 4.777e-02 9.931e+01 9.449e+01 +2.979e-01 4.666e-02 4.400e-02 4.784e-02 9.883e+01 9.449e+01 +3.007e-01 4.656e-02 4.400e-02 4.791e-02 9.838e+01 9.449e+01 +3.035e-01 4.647e-02 4.399e-02 4.798e-02 9.795e+01 9.449e+01 +3.063e-01 4.638e-02 4.398e-02 4.805e-02 9.755e+01 9.449e+01 +3.091e-01 4.629e-02 4.397e-02 4.812e-02 9.717e+01 9.449e+01 +3.120e-01 4.620e-02 4.397e-02 4.819e-02 9.682e+01 9.449e+01 +3.149e-01 4.611e-02 4.396e-02 4.826e-02 9.649e+01 9.449e+01 +3.178e-01 4.602e-02 4.396e-02 4.833e-02 9.619e+01 9.449e+01 +3.207e-01 4.593e-02 4.395e-02 4.840e-02 9.592e+01 9.449e+01 +3.237e-01 4.584e-02 4.396e-02 4.844e-02 9.567e+01 9.449e+01 +3.267e-01 4.576e-02 4.395e-02 4.851e-02 9.544e+01 9.449e+01 +3.297e-01 4.567e-02 4.395e-02 4.858e-02 9.524e+01 9.449e+01 +3.327e-01 4.559e-02 4.395e-02 4.864e-02 9.506e+01 9.449e+01 +3.358e-01 4.551e-02 4.394e-02 4.871e-02 9.491e+01 9.449e+01 +3.389e-01 4.542e-02 4.394e-02 4.877e-02 9.478e+01 9.449e+01 +3.421e-01 4.534e-02 4.394e-02 4.884e-02 9.468e+01 9.449e+01 +3.452e-01 4.526e-02 4.389e-02 4.902e-02 9.460e+01 9.449e+01 +3.484e-01 4.518e-02 4.390e-02 4.905e-02 9.454e+01 9.449e+01 +3.517e-01 4.510e-02 4.392e-02 4.907e-02 9.450e+01 9.449e+01 +3.549e-01 4.503e-02 4.393e-02 4.910e-02 9.449e+01 9.449e+01 +3.582e-01 4.495e-02 4.393e-02 4.916e-02 9.450e+01 9.449e+01 +3.615e-01 4.487e-02 4.392e-02 4.923e-02 9.454e+01 9.449e+01 +3.648e-01 4.480e-02 4.392e-02 4.930e-02 9.460e+01 9.449e+01 +3.682e-01 4.472e-02 4.392e-02 4.936e-02 9.467e+01 9.449e+01 +3.716e-01 4.465e-02 4.392e-02 4.943e-02 9.477e+01 9.449e+01 +3.751e-01 4.458e-02 4.392e-02 4.949e-02 9.490e+01 9.449e+01 +3.785e-01 4.450e-02 4.392e-02 4.955e-02 9.504e+01 9.449e+01 +3.820e-01 4.443e-02 4.392e-02 4.962e-02 9.521e+01 9.449e+01 +3.856e-01 4.436e-02 4.392e-02 4.968e-02 9.540e+01 9.449e+01 +3.891e-01 4.429e-02 4.392e-02 4.974e-02 9.561e+01 9.449e+01 +3.927e-01 4.422e-02 4.392e-02 4.980e-02 9.584e+01 9.449e+01 +3.964e-01 4.415e-02 4.392e-02 4.986e-02 9.609e+01 9.449e+01 +4.000e-01 4.408e-02 4.392e-02 4.992e-02 9.636e+01 9.449e+01 +4.037e-01 4.402e-02 4.392e-02 4.999e-02 9.665e+01 9.449e+01 +4.075e-01 4.395e-02 4.392e-02 5.005e-02 9.696e+01 9.449e+01 +4.113e-01 4.388e-02 4.392e-02 5.010e-02 9.729e+01 9.449e+01 +4.151e-01 4.382e-02 4.393e-02 5.016e-02 9.765e+01 9.449e+01 +4.189e-01 4.375e-02 4.393e-02 5.022e-02 9.802e+01 9.449e+01 +4.228e-01 4.369e-02 4.393e-02 5.028e-02 9.841e+01 9.449e+01 +4.267e-01 4.363e-02 4.393e-02 5.034e-02 9.882e+01 9.449e+01 +4.306e-01 4.356e-02 4.393e-02 5.040e-02 9.925e+01 9.449e+01 +4.346e-01 4.350e-02 4.394e-02 5.046e-02 9.970e+01 9.449e+01 +4.386e-01 4.344e-02 4.394e-02 5.051e-02 1.002e+02 9.449e+01 +4.427e-01 4.338e-02 4.394e-02 5.057e-02 1.007e+02 9.449e+01 +4.468e-01 4.332e-02 4.394e-02 5.063e-02 1.012e+02 9.449e+01 +4.509e-01 4.326e-02 4.395e-02 5.068e-02 1.017e+02 9.449e+01 +4.551e-01 4.320e-02 4.395e-02 5.075e-02 1.022e+02 9.449e+01 +4.593e-01 4.314e-02 4.395e-02 5.080e-02 1.028e+02 9.449e+01 +4.636e-01 4.308e-02 4.395e-02 5.086e-02 1.034e+02 9.449e+01 +4.679e-01 4.302e-02 4.396e-02 5.091e-02 1.040e+02 9.449e+01 +4.722e-01 4.297e-02 4.396e-02 5.097e-02 1.046e+02 9.449e+01 +4.766e-01 4.291e-02 4.396e-02 5.102e-02 1.052e+02 9.449e+01 +4.810e-01 4.285e-02 4.397e-02 5.108e-02 1.059e+02 9.449e+01 +4.854e-01 4.280e-02 4.397e-02 5.113e-02 1.065e+02 9.449e+01 +4.899e-01 4.274e-02 4.397e-02 5.119e-02 1.072e+02 9.449e+01 +4.944e-01 4.269e-02 4.398e-02 5.121e-02 1.079e+02 9.449e+01 +4.990e-01 4.264e-02 4.399e-02 5.127e-02 1.086e+02 9.449e+01 +5.036e-01 4.258e-02 4.399e-02 5.132e-02 1.093e+02 9.449e+01 +5.083e-01 4.253e-02 4.399e-02 5.137e-02 1.101e+02 9.449e+01 +5.130e-01 4.248e-02 4.400e-02 5.143e-02 1.109e+02 9.449e+01 +5.177e-01 4.243e-02 4.400e-02 5.148e-02 1.116e+02 9.449e+01 +5.225e-01 4.238e-02 4.401e-02 5.153e-02 1.124e+02 9.449e+01 +5.274e-01 4.233e-02 4.401e-02 5.158e-02 1.132e+02 9.449e+01 +5.322e-01 4.227e-02 4.401e-02 5.163e-02 1.141e+02 9.449e+01 +5.372e-01 4.223e-02 4.402e-02 5.168e-02 1.149e+02 9.449e+01 +5.421e-01 4.218e-02 4.402e-02 5.173e-02 1.158e+02 9.449e+01 +5.472e-01 4.213e-02 4.403e-02 5.178e-02 1.166e+02 9.449e+01 +5.522e-01 4.208e-02 4.403e-02 5.183e-02 1.175e+02 9.449e+01 +5.573e-01 4.203e-02 4.404e-02 5.188e-02 1.184e+02 9.449e+01 +5.625e-01 4.198e-02 4.404e-02 5.193e-02 1.193e+02 9.449e+01 diff --git a/results/fitnup107nup358/k1_pm_CI_coeffMax.txt b/results/fitnup107nup358/k1_pm_CI_coeffMax.txt new file mode 100644 index 0000000000000000000000000000000000000000..64dbda42226ce2799126402cce6df0739e3ace73 --- /dev/null +++ b/results/fitnup107nup358/k1_pm_CI_coeffMax.txt @@ -0,0 +1,101 @@ +%k1_pm k2_pm k1_ip k2_ip norm normBest +2.259e-01 4.952e-02 3.548e-02 2.426e-02 1.007e+00 1.019e+00 1.024e+00 1.133e+00 1.083e+02 7.832e+01 +2.280e-01 4.937e-02 3.550e-02 2.432e-02 1.007e+00 1.019e+00 1.024e+00 1.133e+00 1.070e+02 7.832e+01 +2.301e-01 4.921e-02 3.553e-02 2.438e-02 1.007e+00 1.019e+00 1.024e+00 1.132e+00 1.057e+02 7.832e+01 +2.322e-01 4.906e-02 3.555e-02 2.445e-02 1.007e+00 1.019e+00 1.024e+00 1.132e+00 1.045e+02 7.832e+01 +2.344e-01 4.891e-02 3.558e-02 2.449e-02 1.007e+00 1.019e+00 1.024e+00 1.131e+00 1.033e+02 7.832e+01 +2.365e-01 4.877e-02 3.560e-02 2.455e-02 1.006e+00 1.019e+00 1.024e+00 1.131e+00 1.022e+02 7.832e+01 +2.387e-01 4.862e-02 3.563e-02 2.462e-02 1.006e+00 1.019e+00 1.024e+00 1.131e+00 1.010e+02 7.832e+01 +2.409e-01 4.848e-02 3.565e-02 2.468e-02 1.006e+00 1.019e+00 1.024e+00 1.130e+00 9.993e+01 7.832e+01 +2.431e-01 4.833e-02 3.568e-02 2.475e-02 1.006e+00 1.019e+00 1.023e+00 1.130e+00 9.886e+01 7.832e+01 +2.454e-01 4.819e-02 3.571e-02 2.481e-02 1.006e+00 1.019e+00 1.023e+00 1.129e+00 9.782e+01 7.832e+01 +2.477e-01 4.805e-02 3.574e-02 2.488e-02 1.006e+00 1.020e+00 1.023e+00 1.129e+00 9.681e+01 7.832e+01 +2.500e-01 4.792e-02 3.577e-02 2.495e-02 1.006e+00 1.020e+00 1.023e+00 1.128e+00 9.583e+01 7.832e+01 +2.523e-01 4.778e-02 3.580e-02 2.499e-02 1.006e+00 1.020e+00 1.023e+00 1.128e+00 9.487e+01 7.832e+01 +2.546e-01 4.765e-02 3.583e-02 2.506e-02 1.006e+00 1.020e+00 1.023e+00 1.127e+00 9.395e+01 7.832e+01 +2.570e-01 4.751e-02 3.587e-02 2.513e-02 1.006e+00 1.020e+00 1.023e+00 1.127e+00 9.306e+01 7.832e+01 +2.593e-01 4.738e-02 3.590e-02 2.520e-02 1.006e+00 1.020e+00 1.022e+00 1.126e+00 9.219e+01 7.832e+01 +2.617e-01 4.725e-02 3.594e-02 2.532e-02 1.006e+00 1.020e+00 1.022e+00 1.125e+00 9.135e+01 7.832e+01 +2.642e-01 4.713e-02 3.598e-02 2.539e-02 1.006e+00 1.020e+00 1.022e+00 1.125e+00 9.054e+01 7.832e+01 +2.666e-01 4.700e-02 3.601e-02 2.545e-02 1.006e+00 1.020e+00 1.022e+00 1.124e+00 8.976e+01 7.832e+01 +2.691e-01 4.687e-02 3.605e-02 2.552e-02 1.006e+00 1.020e+00 1.022e+00 1.124e+00 8.900e+01 7.832e+01 +2.716e-01 4.675e-02 3.609e-02 2.558e-02 1.006e+00 1.020e+00 1.022e+00 1.124e+00 8.828e+01 7.832e+01 +2.741e-01 4.663e-02 3.613e-02 2.565e-02 1.006e+00 1.020e+00 1.022e+00 1.123e+00 8.758e+01 7.832e+01 +2.766e-01 4.650e-02 3.617e-02 2.572e-02 1.006e+00 1.020e+00 1.021e+00 1.123e+00 8.691e+01 7.832e+01 +2.792e-01 4.638e-02 3.621e-02 2.579e-02 1.006e+00 1.021e+00 1.021e+00 1.122e+00 8.626e+01 7.832e+01 +2.818e-01 4.626e-02 3.624e-02 2.586e-02 1.006e+00 1.021e+00 1.021e+00 1.122e+00 8.564e+01 7.832e+01 +2.844e-01 4.615e-02 3.628e-02 2.593e-02 1.005e+00 1.021e+00 1.021e+00 1.121e+00 8.505e+01 7.832e+01 +2.870e-01 4.603e-02 3.632e-02 2.600e-02 1.005e+00 1.021e+00 1.021e+00 1.120e+00 8.448e+01 7.832e+01 +2.896e-01 4.591e-02 3.637e-02 2.601e-02 1.005e+00 1.021e+00 1.021e+00 1.120e+00 8.394e+01 7.832e+01 +2.923e-01 4.580e-02 3.641e-02 2.608e-02 1.005e+00 1.021e+00 1.020e+00 1.120e+00 8.343e+01 7.832e+01 +2.950e-01 4.569e-02 3.645e-02 2.615e-02 1.005e+00 1.021e+00 1.020e+00 1.119e+00 8.294e+01 7.832e+01 +2.978e-01 4.557e-02 3.649e-02 2.622e-02 1.005e+00 1.021e+00 1.020e+00 1.119e+00 8.248e+01 7.832e+01 +3.005e-01 4.546e-02 3.654e-02 2.629e-02 1.005e+00 1.021e+00 1.020e+00 1.118e+00 8.205e+01 7.832e+01 +3.033e-01 4.535e-02 3.658e-02 2.636e-02 1.005e+00 1.021e+00 1.020e+00 1.118e+00 8.163e+01 7.832e+01 +3.061e-01 4.525e-02 3.662e-02 2.643e-02 1.005e+00 1.021e+00 1.020e+00 1.117e+00 8.125e+01 7.832e+01 +3.089e-01 4.514e-02 3.667e-02 2.650e-02 1.005e+00 1.021e+00 1.019e+00 1.117e+00 8.089e+01 7.832e+01 +3.118e-01 4.503e-02 3.671e-02 2.657e-02 1.005e+00 1.022e+00 1.019e+00 1.116e+00 8.055e+01 7.832e+01 +3.147e-01 4.493e-02 3.676e-02 2.664e-02 1.005e+00 1.022e+00 1.019e+00 1.116e+00 8.024e+01 7.832e+01 +3.176e-01 4.482e-02 3.680e-02 2.671e-02 1.005e+00 1.022e+00 1.019e+00 1.115e+00 7.995e+01 7.832e+01 +3.205e-01 4.472e-02 3.685e-02 2.678e-02 1.005e+00 1.022e+00 1.019e+00 1.115e+00 7.969e+01 7.832e+01 +3.235e-01 4.462e-02 3.689e-02 2.685e-02 1.005e+00 1.022e+00 1.019e+00 1.114e+00 7.945e+01 7.832e+01 +3.265e-01 4.452e-02 3.694e-02 2.693e-02 1.005e+00 1.022e+00 1.019e+00 1.114e+00 7.923e+01 7.832e+01 +3.295e-01 4.442e-02 3.698e-02 2.700e-02 1.005e+00 1.022e+00 1.018e+00 1.113e+00 7.904e+01 7.832e+01 +3.326e-01 4.432e-02 3.703e-02 2.707e-02 1.005e+00 1.022e+00 1.018e+00 1.113e+00 7.887e+01 7.832e+01 +3.356e-01 4.422e-02 3.708e-02 2.714e-02 1.005e+00 1.022e+00 1.018e+00 1.112e+00 7.872e+01 7.832e+01 +3.387e-01 4.412e-02 3.713e-02 2.720e-02 1.005e+00 1.022e+00 1.018e+00 1.112e+00 7.860e+01 7.832e+01 +3.419e-01 4.403e-02 3.716e-02 2.738e-02 1.004e+00 1.022e+00 1.018e+00 1.111e+00 7.850e+01 7.832e+01 +3.450e-01 4.393e-02 3.721e-02 2.744e-02 1.004e+00 1.022e+00 1.018e+00 1.110e+00 7.842e+01 7.832e+01 +3.482e-01 4.384e-02 3.727e-02 2.738e-02 1.004e+00 1.023e+00 1.017e+00 1.111e+00 7.837e+01 7.832e+01 +3.515e-01 4.374e-02 3.731e-02 2.756e-02 1.004e+00 1.023e+00 1.017e+00 1.110e+00 7.833e+01 7.832e+01 +3.547e-01 4.365e-02 3.736e-02 2.758e-02 1.004e+00 1.023e+00 1.017e+00 1.109e+00 7.832e+01 7.832e+01 +3.580e-01 4.356e-02 3.741e-02 2.763e-02 1.004e+00 1.023e+00 1.017e+00 1.109e+00 7.833e+01 7.832e+01 +3.613e-01 4.347e-02 3.746e-02 2.778e-02 1.004e+00 1.023e+00 1.017e+00 1.108e+00 7.837e+01 7.832e+01 +3.646e-01 4.338e-02 3.751e-02 2.776e-02 1.004e+00 1.023e+00 1.016e+00 1.108e+00 7.842e+01 7.832e+01 +3.680e-01 4.329e-02 3.756e-02 2.783e-02 1.004e+00 1.023e+00 1.016e+00 1.108e+00 7.850e+01 7.832e+01 +3.714e-01 4.320e-02 3.761e-02 2.790e-02 1.004e+00 1.023e+00 1.016e+00 1.107e+00 7.859e+01 7.832e+01 +3.749e-01 4.311e-02 3.766e-02 2.797e-02 1.004e+00 1.023e+00 1.016e+00 1.107e+00 7.871e+01 7.832e+01 +3.783e-01 4.303e-02 3.770e-02 2.804e-02 1.004e+00 1.023e+00 1.016e+00 1.106e+00 7.885e+01 7.832e+01 +3.818e-01 4.294e-02 3.775e-02 2.820e-02 1.004e+00 1.023e+00 1.016e+00 1.105e+00 7.901e+01 7.832e+01 +3.854e-01 4.286e-02 3.780e-02 2.827e-02 1.004e+00 1.023e+00 1.015e+00 1.105e+00 7.919e+01 7.832e+01 +3.889e-01 4.277e-02 3.786e-02 2.834e-02 1.004e+00 1.023e+00 1.015e+00 1.104e+00 7.939e+01 7.832e+01 +3.925e-01 4.269e-02 3.791e-02 2.840e-02 1.004e+00 1.023e+00 1.015e+00 1.104e+00 7.961e+01 7.832e+01 +3.962e-01 4.261e-02 3.796e-02 2.847e-02 1.004e+00 1.023e+00 1.015e+00 1.104e+00 7.985e+01 7.832e+01 +3.998e-01 4.252e-02 3.802e-02 2.853e-02 1.004e+00 1.024e+00 1.015e+00 1.103e+00 8.011e+01 7.832e+01 +4.035e-01 4.244e-02 3.807e-02 2.860e-02 1.004e+00 1.024e+00 1.015e+00 1.103e+00 8.039e+01 7.832e+01 +4.073e-01 4.236e-02 3.813e-02 2.866e-02 1.004e+00 1.024e+00 1.014e+00 1.102e+00 8.068e+01 7.832e+01 +4.110e-01 4.228e-02 3.818e-02 2.873e-02 1.004e+00 1.024e+00 1.014e+00 1.102e+00 8.100e+01 7.832e+01 +4.148e-01 4.220e-02 3.824e-02 2.879e-02 1.004e+00 1.024e+00 1.014e+00 1.101e+00 8.134e+01 7.832e+01 +4.187e-01 4.212e-02 3.829e-02 2.885e-02 1.003e+00 1.024e+00 1.014e+00 1.101e+00 8.170e+01 7.832e+01 +4.225e-01 4.205e-02 3.828e-02 2.900e-02 1.003e+00 1.024e+00 1.014e+00 1.100e+00 8.207e+01 7.832e+01 +4.265e-01 4.197e-02 3.833e-02 2.907e-02 1.003e+00 1.024e+00 1.014e+00 1.100e+00 8.246e+01 7.832e+01 +4.304e-01 4.189e-02 3.838e-02 2.915e-02 1.003e+00 1.024e+00 1.014e+00 1.099e+00 8.287e+01 7.832e+01 +4.344e-01 4.182e-02 3.843e-02 2.923e-02 1.003e+00 1.024e+00 1.013e+00 1.099e+00 8.330e+01 7.832e+01 +4.384e-01 4.174e-02 3.847e-02 2.931e-02 1.003e+00 1.024e+00 1.013e+00 1.098e+00 8.375e+01 7.832e+01 +4.425e-01 4.167e-02 3.859e-02 2.928e-02 1.003e+00 1.024e+00 1.013e+00 1.098e+00 8.422e+01 7.832e+01 +4.466e-01 4.159e-02 3.865e-02 2.934e-02 1.003e+00 1.024e+00 1.013e+00 1.098e+00 8.470e+01 7.832e+01 +4.507e-01 4.152e-02 3.870e-02 2.941e-02 1.003e+00 1.025e+00 1.013e+00 1.098e+00 8.520e+01 7.832e+01 +4.549e-01 4.145e-02 3.875e-02 2.948e-02 1.003e+00 1.025e+00 1.012e+00 1.097e+00 8.572e+01 7.832e+01 +4.591e-01 4.138e-02 3.881e-02 2.954e-02 1.003e+00 1.025e+00 1.012e+00 1.097e+00 8.626e+01 7.832e+01 +4.633e-01 4.131e-02 3.879e-02 2.972e-02 1.003e+00 1.025e+00 1.012e+00 1.096e+00 8.681e+01 7.832e+01 +4.676e-01 4.123e-02 3.884e-02 2.979e-02 1.003e+00 1.025e+00 1.012e+00 1.096e+00 8.738e+01 7.832e+01 +4.719e-01 4.116e-02 3.889e-02 2.986e-02 1.003e+00 1.025e+00 1.012e+00 1.095e+00 8.796e+01 7.832e+01 +4.763e-01 4.110e-02 3.894e-02 2.993e-02 1.003e+00 1.025e+00 1.012e+00 1.095e+00 8.857e+01 7.832e+01 +4.807e-01 4.103e-02 3.899e-02 3.001e-02 1.003e+00 1.025e+00 1.011e+00 1.094e+00 8.919e+01 7.832e+01 +4.851e-01 4.096e-02 3.904e-02 3.008e-02 1.003e+00 1.025e+00 1.011e+00 1.094e+00 8.982e+01 7.832e+01 +4.896e-01 4.089e-02 3.909e-02 3.015e-02 1.003e+00 1.025e+00 1.011e+00 1.093e+00 9.047e+01 7.832e+01 +4.942e-01 4.082e-02 3.921e-02 3.012e-02 1.003e+00 1.025e+00 1.011e+00 1.093e+00 9.114e+01 7.832e+01 +4.987e-01 4.076e-02 3.927e-02 3.018e-02 1.003e+00 1.025e+00 1.011e+00 1.093e+00 9.182e+01 7.832e+01 +5.034e-01 4.069e-02 3.932e-02 3.025e-02 1.003e+00 1.025e+00 1.010e+00 1.093e+00 9.252e+01 7.832e+01 +5.080e-01 4.063e-02 3.938e-02 3.032e-02 1.003e+00 1.025e+00 1.010e+00 1.092e+00 9.324e+01 7.832e+01 +5.127e-01 4.056e-02 3.943e-02 3.038e-02 1.003e+00 1.026e+00 1.010e+00 1.092e+00 9.397e+01 7.832e+01 +5.175e-01 4.050e-02 3.948e-02 3.045e-02 1.003e+00 1.026e+00 1.010e+00 1.091e+00 9.471e+01 7.832e+01 +5.222e-01 4.043e-02 3.954e-02 3.051e-02 1.003e+00 1.026e+00 1.010e+00 1.091e+00 9.547e+01 7.832e+01 +5.271e-01 4.037e-02 3.959e-02 3.058e-02 1.003e+00 1.026e+00 1.010e+00 1.091e+00 9.625e+01 7.832e+01 +5.320e-01 4.031e-02 3.965e-02 3.064e-02 1.002e+00 1.026e+00 1.009e+00 1.090e+00 9.704e+01 7.832e+01 +5.369e-01 4.025e-02 3.962e-02 3.081e-02 1.002e+00 1.026e+00 1.009e+00 1.090e+00 9.784e+01 7.832e+01 +5.418e-01 4.018e-02 3.967e-02 3.088e-02 1.002e+00 1.026e+00 1.009e+00 1.089e+00 9.866e+01 7.832e+01 +5.469e-01 4.012e-02 3.972e-02 3.095e-02 1.002e+00 1.026e+00 1.009e+00 1.089e+00 9.950e+01 7.832e+01 +5.519e-01 4.006e-02 3.977e-02 3.102e-02 1.002e+00 1.026e+00 1.009e+00 1.088e+00 1.003e+02 7.832e+01 +5.570e-01 4.000e-02 3.982e-02 3.109e-02 1.002e+00 1.026e+00 1.009e+00 1.088e+00 1.012e+02 7.832e+01 +5.622e-01 3.994e-02 3.987e-02 3.115e-02 1.002e+00 1.026e+00 1.009e+00 1.088e+00 1.021e+02 7.832e+01 diff --git a/results/fitnup107nup358/k2_ip_CI.txt b/results/fitnup107nup358/k2_ip_CI.txt new file mode 100644 index 0000000000000000000000000000000000000000..a4f74d3d64030bea5ba8f7c78b8e1848bd054ebe --- /dev/null +++ b/results/fitnup107nup358/k2_ip_CI.txt @@ -0,0 +1,101 @@ +%k1_pm k2_pm k1_ip k2_ip norm normBest +3.440e-01 4.664e-02 5.439e-02 3.127e-02 1.089e+02 9.449e+01 +3.442e-01 4.660e-02 5.410e-02 3.156e-02 1.082e+02 9.449e+01 +3.445e-01 4.656e-02 5.382e-02 3.185e-02 1.075e+02 9.449e+01 +3.448e-01 4.653e-02 5.354e-02 3.214e-02 1.069e+02 9.449e+01 +3.450e-01 4.649e-02 5.326e-02 3.244e-02 1.063e+02 9.449e+01 +3.453e-01 4.645e-02 5.299e-02 3.274e-02 1.057e+02 9.449e+01 +3.456e-01 4.642e-02 5.272e-02 3.304e-02 1.051e+02 9.449e+01 +3.458e-01 4.638e-02 5.245e-02 3.335e-02 1.045e+02 9.449e+01 +3.461e-01 4.634e-02 5.219e-02 3.366e-02 1.040e+02 9.449e+01 +3.464e-01 4.631e-02 5.193e-02 3.397e-02 1.035e+02 9.449e+01 +3.466e-01 4.627e-02 5.167e-02 3.428e-02 1.030e+02 9.449e+01 +3.469e-01 4.624e-02 5.142e-02 3.460e-02 1.025e+02 9.449e+01 +3.471e-01 4.620e-02 5.117e-02 3.492e-02 1.020e+02 9.449e+01 +3.474e-01 4.617e-02 5.090e-02 3.524e-02 1.016e+02 9.449e+01 +3.476e-01 4.613e-02 5.065e-02 3.557e-02 1.011e+02 9.449e+01 +3.478e-01 4.610e-02 5.042e-02 3.590e-02 1.007e+02 9.449e+01 +3.481e-01 4.606e-02 5.018e-02 3.623e-02 1.003e+02 9.449e+01 +3.483e-01 4.603e-02 4.995e-02 3.657e-02 9.991e+01 9.449e+01 +3.486e-01 4.600e-02 4.972e-02 3.690e-02 9.954e+01 9.449e+01 +3.488e-01 4.596e-02 4.949e-02 3.725e-02 9.918e+01 9.449e+01 +3.490e-01 4.593e-02 4.927e-02 3.759e-02 9.884e+01 9.449e+01 +3.493e-01 4.589e-02 4.905e-02 3.794e-02 9.852e+01 9.449e+01 +3.495e-01 4.586e-02 4.883e-02 3.829e-02 9.821e+01 9.449e+01 +3.497e-01 4.583e-02 4.861e-02 3.864e-02 9.791e+01 9.449e+01 +3.500e-01 4.579e-02 4.840e-02 3.900e-02 9.763e+01 9.449e+01 +3.502e-01 4.576e-02 4.819e-02 3.936e-02 9.737e+01 9.449e+01 +3.504e-01 4.573e-02 4.799e-02 3.973e-02 9.711e+01 9.449e+01 +3.506e-01 4.570e-02 4.778e-02 4.009e-02 9.687e+01 9.449e+01 +3.508e-01 4.566e-02 4.758e-02 4.047e-02 9.664e+01 9.449e+01 +3.511e-01 4.563e-02 4.738e-02 4.084e-02 9.643e+01 9.449e+01 +3.513e-01 4.560e-02 4.719e-02 4.122e-02 9.623e+01 9.449e+01 +3.515e-01 4.557e-02 4.699e-02 4.160e-02 9.604e+01 9.449e+01 +3.517e-01 4.554e-02 4.680e-02 4.198e-02 9.586e+01 9.449e+01 +3.519e-01 4.550e-02 4.662e-02 4.237e-02 9.570e+01 9.449e+01 +3.521e-01 4.547e-02 4.643e-02 4.276e-02 9.554e+01 9.449e+01 +3.523e-01 4.544e-02 4.625e-02 4.316e-02 9.540e+01 9.449e+01 +3.525e-01 4.541e-02 4.607e-02 4.356e-02 9.527e+01 9.449e+01 +3.527e-01 4.538e-02 4.591e-02 4.396e-02 9.515e+01 9.449e+01 +3.529e-01 4.535e-02 4.573e-02 4.437e-02 9.504e+01 9.449e+01 +3.531e-01 4.532e-02 4.556e-02 4.478e-02 9.494e+01 9.449e+01 +3.533e-01 4.529e-02 4.538e-02 4.519e-02 9.485e+01 9.449e+01 +3.535e-01 4.526e-02 4.521e-02 4.561e-02 9.478e+01 9.449e+01 +3.537e-01 4.523e-02 4.505e-02 4.603e-02 9.471e+01 9.449e+01 +3.538e-01 4.520e-02 4.488e-02 4.646e-02 9.465e+01 9.449e+01 +3.540e-01 4.517e-02 4.471e-02 4.689e-02 9.460e+01 9.449e+01 +3.542e-01 4.514e-02 4.454e-02 4.732e-02 9.456e+01 9.449e+01 +3.544e-01 4.511e-02 4.438e-02 4.776e-02 9.453e+01 9.449e+01 +3.546e-01 4.508e-02 4.423e-02 4.820e-02 9.451e+01 9.449e+01 +3.547e-01 4.505e-02 4.409e-02 4.865e-02 9.450e+01 9.449e+01 +3.549e-01 4.503e-02 4.393e-02 4.910e-02 9.449e+01 9.449e+01 +3.551e-01 4.500e-02 4.377e-02 4.955e-02 9.450e+01 9.449e+01 +3.552e-01 4.497e-02 4.364e-02 5.001e-02 9.451e+01 9.449e+01 +3.554e-01 4.494e-02 4.349e-02 5.048e-02 9.453e+01 9.449e+01 +3.556e-01 4.491e-02 4.335e-02 5.094e-02 9.456e+01 9.449e+01 +3.557e-01 4.489e-02 4.319e-02 5.141e-02 9.459e+01 9.449e+01 +3.559e-01 4.486e-02 4.305e-02 5.189e-02 9.464e+01 9.449e+01 +3.560e-01 4.483e-02 4.291e-02 5.237e-02 9.469e+01 9.449e+01 +3.562e-01 4.481e-02 4.277e-02 5.285e-02 9.475e+01 9.449e+01 +3.564e-01 4.478e-02 4.264e-02 5.334e-02 9.481e+01 9.449e+01 +3.565e-01 4.475e-02 4.250e-02 5.384e-02 9.488e+01 9.449e+01 +3.567e-01 4.473e-02 4.237e-02 5.434e-02 9.496e+01 9.449e+01 +3.568e-01 4.470e-02 4.224e-02 5.484e-02 9.505e+01 9.449e+01 +3.570e-01 4.467e-02 4.211e-02 5.535e-02 9.514e+01 9.449e+01 +3.571e-01 4.465e-02 4.199e-02 5.586e-02 9.524e+01 9.449e+01 +3.573e-01 4.462e-02 4.188e-02 5.637e-02 9.534e+01 9.449e+01 +3.574e-01 4.459e-02 4.175e-02 5.690e-02 9.545e+01 9.449e+01 +3.575e-01 4.457e-02 4.163e-02 5.742e-02 9.556e+01 9.449e+01 +3.577e-01 4.454e-02 4.151e-02 5.795e-02 9.568e+01 9.449e+01 +3.578e-01 4.452e-02 4.139e-02 5.849e-02 9.581e+01 9.449e+01 +3.579e-01 4.449e-02 4.128e-02 5.903e-02 9.594e+01 9.449e+01 +3.581e-01 4.447e-02 4.116e-02 5.958e-02 9.608e+01 9.449e+01 +3.582e-01 4.444e-02 4.105e-02 6.013e-02 9.622e+01 9.449e+01 +3.583e-01 4.442e-02 4.094e-02 6.069e-02 9.637e+01 9.449e+01 +3.585e-01 4.439e-02 4.083e-02 6.125e-02 9.652e+01 9.449e+01 +3.586e-01 4.437e-02 4.072e-02 6.181e-02 9.668e+01 9.449e+01 +3.587e-01 4.435e-02 4.061e-02 6.239e-02 9.684e+01 9.449e+01 +3.588e-01 4.432e-02 4.050e-02 6.296e-02 9.700e+01 9.449e+01 +3.589e-01 4.430e-02 4.040e-02 6.355e-02 9.717e+01 9.449e+01 +3.591e-01 4.427e-02 4.029e-02 6.413e-02 9.735e+01 9.449e+01 +3.592e-01 4.425e-02 4.019e-02 6.473e-02 9.753e+01 9.449e+01 +3.593e-01 4.423e-02 4.009e-02 6.533e-02 9.771e+01 9.449e+01 +3.594e-01 4.420e-02 3.999e-02 6.593e-02 9.790e+01 9.449e+01 +3.595e-01 4.418e-02 3.989e-02 6.654e-02 9.809e+01 9.449e+01 +3.596e-01 4.416e-02 3.979e-02 6.716e-02 9.828e+01 9.449e+01 +3.597e-01 4.414e-02 3.969e-02 6.778e-02 9.848e+01 9.449e+01 +3.598e-01 4.411e-02 3.960e-02 6.840e-02 9.868e+01 9.449e+01 +3.599e-01 4.409e-02 3.951e-02 6.904e-02 9.888e+01 9.449e+01 +3.600e-01 4.407e-02 3.941e-02 6.968e-02 9.909e+01 9.449e+01 +3.601e-01 4.405e-02 3.932e-02 7.032e-02 9.930e+01 9.449e+01 +3.602e-01 4.403e-02 3.923e-02 7.097e-02 9.951e+01 9.449e+01 +3.603e-01 4.400e-02 3.914e-02 7.163e-02 9.973e+01 9.449e+01 +3.604e-01 4.398e-02 3.905e-02 7.229e-02 9.995e+01 9.449e+01 +3.605e-01 4.396e-02 3.896e-02 7.296e-02 1.002e+02 9.449e+01 +3.606e-01 4.394e-02 3.888e-02 7.363e-02 1.004e+02 9.449e+01 +3.607e-01 4.392e-02 3.879e-02 7.432e-02 1.006e+02 9.449e+01 +3.608e-01 4.390e-02 3.871e-02 7.500e-02 1.009e+02 9.449e+01 +3.609e-01 4.388e-02 3.863e-02 7.570e-02 1.011e+02 9.449e+01 +3.610e-01 4.386e-02 3.854e-02 7.640e-02 1.013e+02 9.449e+01 +3.610e-01 4.384e-02 3.846e-02 7.711e-02 1.016e+02 9.449e+01 +3.611e-01 4.381e-02 3.838e-02 7.782e-02 1.018e+02 9.449e+01 diff --git a/results/fitnup107nup358/k2_ip_CI_coeffMax.txt b/results/fitnup107nup358/k2_ip_CI_coeffMax.txt new file mode 100644 index 0000000000000000000000000000000000000000..f2c192d0faf857628bf3a22cbaf740464bca186d --- /dev/null +++ b/results/fitnup107nup358/k2_ip_CI_coeffMax.txt @@ -0,0 +1,101 @@ +%k1_pm k2_pm k1_ip k2_ip norm normBest +3.522e-01 4.317e-02 3.913e-02 1.757e-02 1.004e+00 1.034e+00 1.013e+00 1.193e+00 7.979e+01 7.832e+01 +3.522e-01 4.318e-02 3.910e-02 1.774e-02 1.004e+00 1.034e+00 1.014e+00 1.191e+00 7.973e+01 7.832e+01 +3.523e-01 4.320e-02 3.908e-02 1.790e-02 1.004e+00 1.034e+00 1.014e+00 1.189e+00 7.968e+01 7.832e+01 +3.523e-01 4.322e-02 3.905e-02 1.807e-02 1.004e+00 1.034e+00 1.014e+00 1.187e+00 7.963e+01 7.832e+01 +3.524e-01 4.324e-02 3.903e-02 1.823e-02 1.004e+00 1.033e+00 1.014e+00 1.185e+00 7.958e+01 7.832e+01 +3.524e-01 4.326e-02 3.901e-02 1.840e-02 1.004e+00 1.033e+00 1.014e+00 1.183e+00 7.953e+01 7.832e+01 +3.524e-01 4.327e-02 3.899e-02 1.857e-02 1.004e+00 1.033e+00 1.014e+00 1.181e+00 7.948e+01 7.832e+01 +3.525e-01 4.329e-02 3.900e-02 1.875e-02 1.004e+00 1.032e+00 1.014e+00 1.179e+00 7.943e+01 7.832e+01 +3.525e-01 4.330e-02 3.894e-02 1.892e-02 1.004e+00 1.032e+00 1.014e+00 1.177e+00 7.938e+01 7.832e+01 +3.525e-01 4.332e-02 3.891e-02 1.909e-02 1.004e+00 1.032e+00 1.014e+00 1.175e+00 7.934e+01 7.832e+01 +3.526e-01 4.334e-02 3.889e-02 1.927e-02 1.004e+00 1.032e+00 1.014e+00 1.173e+00 7.929e+01 7.832e+01 +3.526e-01 4.335e-02 3.886e-02 1.945e-02 1.004e+00 1.031e+00 1.014e+00 1.171e+00 7.925e+01 7.832e+01 +3.527e-01 4.337e-02 3.883e-02 1.963e-02 1.004e+00 1.031e+00 1.014e+00 1.169e+00 7.920e+01 7.832e+01 +3.527e-01 4.338e-02 3.880e-02 1.981e-02 1.004e+00 1.031e+00 1.014e+00 1.167e+00 7.916e+01 7.832e+01 +3.528e-01 4.339e-02 3.878e-02 1.999e-02 1.004e+00 1.030e+00 1.014e+00 1.165e+00 7.912e+01 7.832e+01 +3.528e-01 4.341e-02 3.875e-02 2.018e-02 1.004e+00 1.030e+00 1.014e+00 1.164e+00 7.907e+01 7.832e+01 +3.528e-01 4.342e-02 3.872e-02 2.037e-02 1.004e+00 1.030e+00 1.014e+00 1.162e+00 7.903e+01 7.832e+01 +3.529e-01 4.343e-02 3.869e-02 2.055e-02 1.004e+00 1.030e+00 1.014e+00 1.160e+00 7.899e+01 7.832e+01 +3.529e-01 4.345e-02 3.866e-02 2.074e-02 1.004e+00 1.029e+00 1.014e+00 1.158e+00 7.896e+01 7.832e+01 +3.530e-01 4.346e-02 3.866e-02 2.094e-02 1.004e+00 1.029e+00 1.014e+00 1.156e+00 7.892e+01 7.832e+01 +3.531e-01 4.347e-02 3.863e-02 2.113e-02 1.004e+00 1.029e+00 1.014e+00 1.154e+00 7.888e+01 7.832e+01 +3.531e-01 4.348e-02 3.859e-02 2.133e-02 1.004e+00 1.029e+00 1.014e+00 1.153e+00 7.885e+01 7.832e+01 +3.531e-01 4.349e-02 3.856e-02 2.152e-02 1.004e+00 1.028e+00 1.015e+00 1.151e+00 7.881e+01 7.832e+01 +3.532e-01 4.350e-02 3.852e-02 2.172e-02 1.004e+00 1.028e+00 1.015e+00 1.149e+00 7.878e+01 7.832e+01 +3.533e-01 4.351e-02 3.849e-02 2.192e-02 1.004e+00 1.028e+00 1.015e+00 1.147e+00 7.875e+01 7.832e+01 +3.533e-01 4.352e-02 3.845e-02 2.213e-02 1.004e+00 1.028e+00 1.015e+00 1.146e+00 7.871e+01 7.832e+01 +3.534e-01 4.353e-02 3.842e-02 2.233e-02 1.004e+00 1.027e+00 1.015e+00 1.144e+00 7.868e+01 7.832e+01 +3.534e-01 4.354e-02 3.838e-02 2.254e-02 1.004e+00 1.027e+00 1.015e+00 1.142e+00 7.865e+01 7.832e+01 +3.535e-01 4.355e-02 3.834e-02 2.275e-02 1.004e+00 1.027e+00 1.015e+00 1.141e+00 7.863e+01 7.832e+01 +3.535e-01 4.356e-02 3.830e-02 2.296e-02 1.004e+00 1.027e+00 1.015e+00 1.139e+00 7.860e+01 7.832e+01 +3.536e-01 4.357e-02 3.826e-02 2.317e-02 1.004e+00 1.027e+00 1.015e+00 1.137e+00 7.857e+01 7.832e+01 +3.536e-01 4.357e-02 3.822e-02 2.338e-02 1.004e+00 1.026e+00 1.015e+00 1.136e+00 7.855e+01 7.832e+01 +3.537e-01 4.358e-02 3.818e-02 2.360e-02 1.004e+00 1.026e+00 1.015e+00 1.134e+00 7.853e+01 7.832e+01 +3.537e-01 4.359e-02 3.814e-02 2.382e-02 1.004e+00 1.026e+00 1.015e+00 1.133e+00 7.850e+01 7.832e+01 +3.538e-01 4.360e-02 3.810e-02 2.404e-02 1.004e+00 1.026e+00 1.015e+00 1.131e+00 7.848e+01 7.832e+01 +3.539e-01 4.360e-02 3.805e-02 2.426e-02 1.004e+00 1.025e+00 1.016e+00 1.129e+00 7.846e+01 7.832e+01 +3.539e-01 4.361e-02 3.801e-02 2.448e-02 1.004e+00 1.025e+00 1.016e+00 1.128e+00 7.844e+01 7.832e+01 +3.539e-01 4.361e-02 3.793e-02 2.471e-02 1.004e+00 1.025e+00 1.016e+00 1.126e+00 7.843e+01 7.832e+01 +3.540e-01 4.362e-02 3.789e-02 2.494e-02 1.004e+00 1.025e+00 1.016e+00 1.125e+00 7.841e+01 7.832e+01 +3.541e-01 4.362e-02 3.784e-02 2.517e-02 1.004e+00 1.025e+00 1.016e+00 1.123e+00 7.840e+01 7.832e+01 +3.541e-01 4.363e-02 3.780e-02 2.540e-02 1.004e+00 1.024e+00 1.016e+00 1.122e+00 7.838e+01 7.832e+01 +3.542e-01 4.363e-02 3.775e-02 2.564e-02 1.004e+00 1.024e+00 1.016e+00 1.120e+00 7.837e+01 7.832e+01 +3.543e-01 4.364e-02 3.771e-02 2.588e-02 1.004e+00 1.024e+00 1.016e+00 1.119e+00 7.836e+01 7.832e+01 +3.543e-01 4.364e-02 3.769e-02 2.612e-02 1.004e+00 1.024e+00 1.016e+00 1.117e+00 7.835e+01 7.832e+01 +3.544e-01 4.364e-02 3.764e-02 2.636e-02 1.004e+00 1.024e+00 1.016e+00 1.116e+00 7.834e+01 7.832e+01 +3.545e-01 4.364e-02 3.758e-02 2.660e-02 1.004e+00 1.023e+00 1.017e+00 1.115e+00 7.833e+01 7.832e+01 +3.545e-01 4.365e-02 3.750e-02 2.685e-02 1.004e+00 1.023e+00 1.017e+00 1.113e+00 7.833e+01 7.832e+01 +3.546e-01 4.365e-02 3.749e-02 2.710e-02 1.004e+00 1.023e+00 1.017e+00 1.112e+00 7.833e+01 7.832e+01 +3.546e-01 4.365e-02 3.740e-02 2.735e-02 1.004e+00 1.023e+00 1.017e+00 1.111e+00 7.832e+01 7.832e+01 +3.547e-01 4.365e-02 3.736e-02 2.760e-02 1.004e+00 1.023e+00 1.017e+00 1.109e+00 7.832e+01 7.832e+01 +3.548e-01 4.365e-02 3.731e-02 2.785e-02 1.004e+00 1.022e+00 1.017e+00 1.108e+00 7.832e+01 7.832e+01 +3.548e-01 4.365e-02 3.727e-02 2.811e-02 1.004e+00 1.022e+00 1.017e+00 1.107e+00 7.833e+01 7.832e+01 +3.549e-01 4.365e-02 3.718e-02 2.837e-02 1.004e+00 1.022e+00 1.017e+00 1.105e+00 7.833e+01 7.832e+01 +3.550e-01 4.365e-02 3.712e-02 2.863e-02 1.004e+00 1.022e+00 1.018e+00 1.104e+00 7.834e+01 7.832e+01 +3.550e-01 4.365e-02 3.706e-02 2.890e-02 1.004e+00 1.022e+00 1.018e+00 1.103e+00 7.834e+01 7.832e+01 +3.551e-01 4.365e-02 3.700e-02 2.917e-02 1.004e+00 1.022e+00 1.018e+00 1.102e+00 7.835e+01 7.832e+01 +3.551e-01 4.365e-02 3.694e-02 2.944e-02 1.004e+00 1.021e+00 1.018e+00 1.101e+00 7.836e+01 7.832e+01 +3.552e-01 4.365e-02 3.691e-02 2.971e-02 1.004e+00 1.021e+00 1.018e+00 1.099e+00 7.837e+01 7.832e+01 +3.553e-01 4.365e-02 3.685e-02 2.998e-02 1.004e+00 1.021e+00 1.018e+00 1.098e+00 7.839e+01 7.832e+01 +3.554e-01 4.365e-02 3.679e-02 3.026e-02 1.004e+00 1.021e+00 1.018e+00 1.097e+00 7.840e+01 7.832e+01 +3.554e-01 4.364e-02 3.673e-02 3.054e-02 1.005e+00 1.021e+00 1.018e+00 1.096e+00 7.842e+01 7.832e+01 +3.555e-01 4.364e-02 3.667e-02 3.082e-02 1.005e+00 1.021e+00 1.019e+00 1.095e+00 7.844e+01 7.832e+01 +3.556e-01 4.364e-02 3.661e-02 3.111e-02 1.005e+00 1.020e+00 1.019e+00 1.093e+00 7.846e+01 7.832e+01 +3.556e-01 4.363e-02 3.654e-02 3.140e-02 1.005e+00 1.020e+00 1.019e+00 1.092e+00 7.848e+01 7.832e+01 +3.557e-01 4.363e-02 3.648e-02 3.169e-02 1.005e+00 1.020e+00 1.019e+00 1.091e+00 7.850e+01 7.832e+01 +3.558e-01 4.363e-02 3.641e-02 3.198e-02 1.005e+00 1.020e+00 1.019e+00 1.090e+00 7.853e+01 7.832e+01 +3.558e-01 4.362e-02 3.635e-02 3.228e-02 1.005e+00 1.020e+00 1.019e+00 1.089e+00 7.856e+01 7.832e+01 +3.559e-01 4.362e-02 3.628e-02 3.258e-02 1.005e+00 1.020e+00 1.019e+00 1.088e+00 7.859e+01 7.832e+01 +3.560e-01 4.361e-02 3.621e-02 3.288e-02 1.005e+00 1.020e+00 1.020e+00 1.087e+00 7.862e+01 7.832e+01 +3.560e-01 4.361e-02 3.611e-02 3.318e-02 1.005e+00 1.020e+00 1.020e+00 1.086e+00 7.865e+01 7.832e+01 +3.561e-01 4.360e-02 3.604e-02 3.349e-02 1.005e+00 1.019e+00 1.020e+00 1.085e+00 7.868e+01 7.832e+01 +3.561e-01 4.360e-02 3.598e-02 3.380e-02 1.005e+00 1.019e+00 1.020e+00 1.084e+00 7.872e+01 7.832e+01 +3.562e-01 4.359e-02 3.591e-02 3.411e-02 1.005e+00 1.019e+00 1.020e+00 1.083e+00 7.876e+01 7.832e+01 +3.563e-01 4.358e-02 3.584e-02 3.443e-02 1.005e+00 1.019e+00 1.020e+00 1.082e+00 7.880e+01 7.832e+01 +3.563e-01 4.358e-02 3.577e-02 3.474e-02 1.005e+00 1.019e+00 1.021e+00 1.081e+00 7.884e+01 7.832e+01 +3.564e-01 4.357e-02 3.570e-02 3.507e-02 1.005e+00 1.019e+00 1.021e+00 1.080e+00 7.888e+01 7.832e+01 +3.565e-01 4.356e-02 3.564e-02 3.539e-02 1.005e+00 1.019e+00 1.021e+00 1.079e+00 7.893e+01 7.832e+01 +3.565e-01 4.355e-02 3.557e-02 3.572e-02 1.005e+00 1.019e+00 1.021e+00 1.078e+00 7.898e+01 7.832e+01 +3.566e-01 4.355e-02 3.550e-02 3.605e-02 1.005e+00 1.018e+00 1.021e+00 1.078e+00 7.903e+01 7.832e+01 +3.567e-01 4.354e-02 3.543e-02 3.638e-02 1.005e+00 1.018e+00 1.021e+00 1.077e+00 7.908e+01 7.832e+01 +3.567e-01 4.353e-02 3.536e-02 3.672e-02 1.005e+00 1.018e+00 1.022e+00 1.076e+00 7.913e+01 7.832e+01 +3.568e-01 4.352e-02 3.528e-02 3.706e-02 1.005e+00 1.018e+00 1.022e+00 1.075e+00 7.918e+01 7.832e+01 +3.569e-01 4.351e-02 3.521e-02 3.740e-02 1.005e+00 1.018e+00 1.022e+00 1.074e+00 7.924e+01 7.832e+01 +3.569e-01 4.351e-02 3.510e-02 3.775e-02 1.005e+00 1.018e+00 1.022e+00 1.073e+00 7.930e+01 7.832e+01 +3.570e-01 4.350e-02 3.503e-02 3.810e-02 1.005e+00 1.018e+00 1.022e+00 1.073e+00 7.936e+01 7.832e+01 +3.570e-01 4.349e-02 3.495e-02 3.845e-02 1.005e+00 1.018e+00 1.023e+00 1.072e+00 7.942e+01 7.832e+01 +3.571e-01 4.348e-02 3.487e-02 3.881e-02 1.005e+00 1.018e+00 1.023e+00 1.071e+00 7.949e+01 7.832e+01 +3.571e-01 4.347e-02 3.479e-02 3.916e-02 1.005e+00 1.018e+00 1.023e+00 1.070e+00 7.955e+01 7.832e+01 +3.572e-01 4.346e-02 3.474e-02 3.953e-02 1.005e+00 1.018e+00 1.023e+00 1.069e+00 7.962e+01 7.832e+01 +3.573e-01 4.345e-02 3.466e-02 3.989e-02 1.005e+00 1.017e+00 1.023e+00 1.069e+00 7.969e+01 7.832e+01 +3.574e-01 4.344e-02 3.458e-02 4.026e-02 1.005e+00 1.017e+00 1.024e+00 1.068e+00 7.976e+01 7.832e+01 +3.574e-01 4.342e-02 3.450e-02 4.063e-02 1.005e+00 1.017e+00 1.024e+00 1.067e+00 7.984e+01 7.832e+01 +3.575e-01 4.341e-02 3.443e-02 4.101e-02 1.005e+00 1.017e+00 1.024e+00 1.067e+00 7.991e+01 7.832e+01 +3.575e-01 4.340e-02 3.435e-02 4.139e-02 1.005e+00 1.017e+00 1.024e+00 1.066e+00 7.999e+01 7.832e+01 +3.576e-01 4.339e-02 3.427e-02 4.177e-02 1.005e+00 1.017e+00 1.024e+00 1.065e+00 8.007e+01 7.832e+01 +3.576e-01 4.338e-02 3.419e-02 4.216e-02 1.005e+00 1.017e+00 1.025e+00 1.065e+00 8.015e+01 7.832e+01 +3.577e-01 4.337e-02 3.411e-02 4.255e-02 1.005e+00 1.017e+00 1.025e+00 1.064e+00 8.023e+01 7.832e+01 +3.578e-01 4.336e-02 3.403e-02 4.294e-02 1.005e+00 1.017e+00 1.025e+00 1.063e+00 8.032e+01 7.832e+01 +3.578e-01 4.335e-02 3.391e-02 4.334e-02 1.006e+00 1.017e+00 1.025e+00 1.063e+00 8.041e+01 7.832e+01 +3.578e-01 4.333e-02 3.383e-02 4.374e-02 1.006e+00 1.017e+00 1.026e+00 1.062e+00 8.049e+01 7.832e+01 diff --git a/results/fitnup107nup358/k2_pm_CI.txt b/results/fitnup107nup358/k2_pm_CI.txt new file mode 100644 index 0000000000000000000000000000000000000000..9bda277e14e64c46786e08240bc64a500e94b176 --- /dev/null +++ b/results/fitnup107nup358/k2_pm_CI.txt @@ -0,0 +1,101 @@ +%k1_pm k2_pm k1_ip k2_ip norm normBest +5.630e-01 2.867e-02 4.830e-02 1.414e-01 3.538e+02 9.449e+01 +5.587e-01 2.894e-02 4.821e-02 1.360e-01 3.413e+02 9.449e+01 +5.544e-01 2.921e-02 4.811e-02 1.310e-01 3.293e+02 9.449e+01 +5.501e-01 2.948e-02 4.801e-02 1.264e-01 3.176e+02 9.449e+01 +5.459e-01 2.975e-02 4.791e-02 1.222e-01 3.063e+02 9.449e+01 +5.416e-01 3.002e-02 4.780e-02 1.182e-01 2.955e+02 9.449e+01 +5.373e-01 3.030e-02 4.769e-02 1.145e-01 2.849e+02 9.449e+01 +5.330e-01 3.058e-02 4.758e-02 1.110e-01 2.748e+02 9.449e+01 +5.287e-01 3.087e-02 4.746e-02 1.077e-01 2.649e+02 9.449e+01 +5.244e-01 3.115e-02 4.735e-02 1.046e-01 2.555e+02 9.449e+01 +5.200e-01 3.144e-02 4.723e-02 1.017e-01 2.463e+02 9.449e+01 +5.157e-01 3.173e-02 4.712e-02 9.888e-02 2.376e+02 9.449e+01 +5.114e-01 3.202e-02 4.700e-02 9.624e-02 2.291e+02 9.449e+01 +5.071e-01 3.232e-02 4.689e-02 9.373e-02 2.209e+02 9.449e+01 +5.028e-01 3.262e-02 4.677e-02 9.133e-02 2.131e+02 9.449e+01 +4.985e-01 3.292e-02 4.666e-02 8.905e-02 2.056e+02 9.449e+01 +4.942e-01 3.323e-02 4.655e-02 8.687e-02 1.984e+02 9.449e+01 +4.899e-01 3.353e-02 4.636e-02 8.509e-02 1.914e+02 9.449e+01 +4.855e-01 3.384e-02 4.624e-02 8.312e-02 1.848e+02 9.449e+01 +4.812e-01 3.416e-02 4.613e-02 8.124e-02 1.784e+02 9.449e+01 +4.769e-01 3.447e-02 4.602e-02 7.944e-02 1.723e+02 9.449e+01 +4.727e-01 3.479e-02 4.601e-02 7.739e-02 1.665e+02 9.449e+01 +4.684e-01 3.511e-02 4.591e-02 7.571e-02 1.610e+02 9.449e+01 +4.640e-01 3.544e-02 4.573e-02 7.438e-02 1.557e+02 9.449e+01 +4.597e-01 3.577e-02 4.563e-02 7.287e-02 1.507e+02 9.449e+01 +4.555e-01 3.610e-02 4.552e-02 7.141e-02 1.459e+02 9.449e+01 +4.512e-01 3.643e-02 4.550e-02 6.982e-02 1.414e+02 9.449e+01 +4.469e-01 3.677e-02 4.542e-02 6.843e-02 1.371e+02 9.449e+01 +4.426e-01 3.711e-02 4.530e-02 6.720e-02 1.330e+02 9.449e+01 +4.384e-01 3.745e-02 4.522e-02 6.595e-02 1.292e+02 9.449e+01 +4.341e-01 3.780e-02 4.514e-02 6.476e-02 1.256e+02 9.449e+01 +4.298e-01 3.815e-02 4.505e-02 6.361e-02 1.222e+02 9.449e+01 +4.256e-01 3.850e-02 4.498e-02 6.251e-02 1.190e+02 9.449e+01 +4.213e-01 3.886e-02 4.490e-02 6.144e-02 1.161e+02 9.449e+01 +4.171e-01 3.922e-02 4.482e-02 6.044e-02 1.133e+02 9.449e+01 +4.128e-01 3.958e-02 4.474e-02 5.951e-02 1.108e+02 9.449e+01 +4.086e-01 3.995e-02 4.467e-02 5.855e-02 1.084e+02 9.449e+01 +4.044e-01 4.032e-02 4.460e-02 5.765e-02 1.063e+02 9.449e+01 +4.002e-01 4.069e-02 4.455e-02 5.671e-02 1.043e+02 9.449e+01 +3.960e-01 4.106e-02 4.449e-02 5.588e-02 1.026e+02 9.449e+01 +3.918e-01 4.144e-02 4.443e-02 5.506e-02 1.010e+02 9.449e+01 +3.877e-01 4.183e-02 4.436e-02 5.430e-02 9.959e+01 9.449e+01 +3.835e-01 4.222e-02 4.430e-02 5.356e-02 9.837e+01 9.449e+01 +3.794e-01 4.261e-02 4.425e-02 5.284e-02 9.732e+01 9.449e+01 +3.752e-01 4.300e-02 4.419e-02 5.216e-02 9.644e+01 9.449e+01 +3.711e-01 4.340e-02 4.414e-02 5.148e-02 9.573e+01 9.449e+01 +3.671e-01 4.380e-02 4.408e-02 5.085e-02 9.518e+01 9.449e+01 +3.630e-01 4.420e-02 4.403e-02 5.024e-02 9.480e+01 9.449e+01 +3.589e-01 4.461e-02 4.398e-02 4.966e-02 9.457e+01 9.449e+01 +3.549e-01 4.503e-02 4.393e-02 4.910e-02 9.449e+01 9.449e+01 +3.509e-01 4.544e-02 4.389e-02 4.852e-02 9.457e+01 9.449e+01 +3.469e-01 4.586e-02 4.382e-02 4.806e-02 9.479e+01 9.449e+01 +3.429e-01 4.629e-02 4.377e-02 4.757e-02 9.515e+01 9.449e+01 +3.390e-01 4.672e-02 4.372e-02 4.710e-02 9.566e+01 9.449e+01 +3.351e-01 4.715e-02 4.370e-02 4.656e-02 9.631e+01 9.449e+01 +3.312e-01 4.759e-02 4.366e-02 4.611e-02 9.708e+01 9.449e+01 +3.273e-01 4.803e-02 4.362e-02 4.567e-02 9.799e+01 9.449e+01 +3.235e-01 4.847e-02 4.358e-02 4.525e-02 9.903e+01 9.449e+01 +3.197e-01 4.892e-02 4.354e-02 4.485e-02 1.002e+02 9.449e+01 +3.159e-01 4.937e-02 4.351e-02 4.445e-02 1.015e+02 9.449e+01 +3.122e-01 4.983e-02 4.348e-02 4.408e-02 1.029e+02 9.449e+01 +3.085e-01 5.029e-02 4.344e-02 4.373e-02 1.044e+02 9.449e+01 +3.048e-01 5.075e-02 4.341e-02 4.337e-02 1.060e+02 9.449e+01 +3.011e-01 5.122e-02 4.338e-02 4.303e-02 1.078e+02 9.449e+01 +2.975e-01 5.170e-02 4.335e-02 4.270e-02 1.096e+02 9.449e+01 +2.940e-01 5.218e-02 4.332e-02 4.238e-02 1.116e+02 9.449e+01 +2.904e-01 5.266e-02 4.329e-02 4.207e-02 1.137e+02 9.449e+01 +2.869e-01 5.315e-02 4.327e-02 4.177e-02 1.158e+02 9.449e+01 +2.834e-01 5.364e-02 4.324e-02 4.148e-02 1.181e+02 9.449e+01 +2.800e-01 5.413e-02 4.322e-02 4.120e-02 1.204e+02 9.449e+01 +2.766e-01 5.463e-02 4.319e-02 4.095e-02 1.229e+02 9.449e+01 +2.732e-01 5.514e-02 4.317e-02 4.068e-02 1.254e+02 9.449e+01 +2.699e-01 5.565e-02 4.319e-02 4.036e-02 1.280e+02 9.449e+01 +2.665e-01 5.617e-02 4.315e-02 4.016e-02 1.307e+02 9.449e+01 +2.633e-01 5.669e-02 4.314e-02 3.992e-02 1.335e+02 9.449e+01 +2.601e-01 5.721e-02 4.313e-02 3.968e-02 1.363e+02 9.449e+01 +2.569e-01 5.774e-02 4.312e-02 3.945e-02 1.392e+02 9.449e+01 +2.537e-01 5.827e-02 4.307e-02 3.931e-02 1.422e+02 9.449e+01 +2.506e-01 5.881e-02 4.309e-02 3.902e-02 1.453e+02 9.449e+01 +2.475e-01 5.936e-02 4.309e-02 3.882e-02 1.484e+02 9.449e+01 +2.445e-01 5.991e-02 4.308e-02 3.862e-02 1.516e+02 9.449e+01 +2.415e-01 6.046e-02 4.305e-02 3.847e-02 1.549e+02 9.449e+01 +2.385e-01 6.102e-02 4.306e-02 3.828e-02 1.582e+02 9.449e+01 +2.356e-01 6.158e-02 4.306e-02 3.809e-02 1.616e+02 9.449e+01 +2.327e-01 6.215e-02 4.307e-02 3.790e-02 1.650e+02 9.449e+01 +2.299e-01 6.273e-02 4.307e-02 3.772e-02 1.685e+02 9.449e+01 +2.270e-01 6.331e-02 4.308e-02 3.755e-02 1.720e+02 9.449e+01 +2.243e-01 6.390e-02 4.309e-02 3.738e-02 1.756e+02 9.449e+01 +2.215e-01 6.449e-02 4.311e-02 3.721e-02 1.792e+02 9.449e+01 +2.188e-01 6.508e-02 4.312e-02 3.705e-02 1.828e+02 9.449e+01 +2.162e-01 6.569e-02 4.313e-02 3.690e-02 1.865e+02 9.449e+01 +2.136e-01 6.629e-02 4.315e-02 3.675e-02 1.903e+02 9.449e+01 +2.110e-01 6.691e-02 4.317e-02 3.660e-02 1.941e+02 9.449e+01 +2.084e-01 6.753e-02 4.318e-02 3.645e-02 1.979e+02 9.449e+01 +2.059e-01 6.815e-02 4.320e-02 3.631e-02 2.017e+02 9.449e+01 +2.035e-01 6.878e-02 4.323e-02 3.618e-02 2.056e+02 9.449e+01 +2.010e-01 6.942e-02 4.325e-02 3.605e-02 2.095e+02 9.449e+01 +1.986e-01 7.006e-02 4.327e-02 3.592e-02 2.134e+02 9.449e+01 +1.963e-01 7.071e-02 4.330e-02 3.579e-02 2.173e+02 9.449e+01 +1.940e-01 7.136e-02 4.333e-02 3.567e-02 2.213e+02 9.449e+01 diff --git a/results/fitnup107nup358/k2_pm_CI_coeffMax.txt b/results/fitnup107nup358/k2_pm_CI_coeffMax.txt new file mode 100644 index 0000000000000000000000000000000000000000..dcf20f83572ce5e41221dbe9d84f920552358da6 --- /dev/null +++ b/results/fitnup107nup358/k2_pm_CI_coeffMax.txt @@ -0,0 +1,101 @@ +%k1_pm k2_pm k1_ip k2_ip norm normBest +4.749e-01 2.780e-02 3.060e-02 2.000e+00 1.006e+00 1.081e+00 1.034e+00 1.035e+00 1.869e+02 7.832e+01 +4.924e-01 2.805e-02 3.504e-02 1.878e-01 1.004e+00 1.078e+00 1.020e+00 1.034e+00 1.843e+02 7.832e+01 +4.929e-01 2.831e-02 3.587e-02 1.547e-01 1.004e+00 1.076e+00 1.018e+00 1.034e+00 1.806e+02 7.832e+01 +4.921e-01 2.858e-02 3.629e-02 1.385e-01 1.004e+00 1.074e+00 1.017e+00 1.034e+00 1.768e+02 7.832e+01 +4.908e-01 2.884e-02 3.659e-02 1.269e-01 1.004e+00 1.072e+00 1.016e+00 1.034e+00 1.731e+02 7.832e+01 +4.893e-01 2.911e-02 3.682e-02 1.178e-01 1.004e+00 1.070e+00 1.016e+00 1.035e+00 1.695e+02 7.832e+01 +4.876e-01 2.938e-02 3.701e-02 1.103e-01 1.003e+00 1.068e+00 1.015e+00 1.035e+00 1.659e+02 7.832e+01 +4.858e-01 2.965e-02 3.717e-02 1.039e-01 1.003e+00 1.066e+00 1.015e+00 1.035e+00 1.623e+02 7.832e+01 +4.838e-01 2.992e-02 3.730e-02 9.834e-02 1.003e+00 1.065e+00 1.015e+00 1.036e+00 1.588e+02 7.832e+01 +4.818e-01 3.020e-02 3.742e-02 9.324e-02 1.003e+00 1.063e+00 1.015e+00 1.036e+00 1.554e+02 7.832e+01 +4.797e-01 3.048e-02 3.753e-02 8.873e-02 1.003e+00 1.061e+00 1.014e+00 1.037e+00 1.520e+02 7.832e+01 +4.775e-01 3.076e-02 3.762e-02 8.459e-02 1.003e+00 1.060e+00 1.014e+00 1.037e+00 1.486e+02 7.832e+01 +4.753e-01 3.105e-02 3.771e-02 8.077e-02 1.003e+00 1.058e+00 1.014e+00 1.038e+00 1.453e+02 7.832e+01 +4.730e-01 3.133e-02 3.779e-02 7.722e-02 1.003e+00 1.056e+00 1.014e+00 1.039e+00 1.421e+02 7.832e+01 +4.706e-01 3.162e-02 3.787e-02 7.391e-02 1.003e+00 1.055e+00 1.014e+00 1.040e+00 1.390e+02 7.832e+01 +4.681e-01 3.192e-02 3.795e-02 7.080e-02 1.003e+00 1.054e+00 1.014e+00 1.041e+00 1.359e+02 7.832e+01 +4.657e-01 3.221e-02 3.802e-02 6.783e-02 1.003e+00 1.052e+00 1.014e+00 1.042e+00 1.328e+02 7.832e+01 +4.631e-01 3.251e-02 3.809e-02 6.507e-02 1.003e+00 1.051e+00 1.014e+00 1.043e+00 1.298e+02 7.832e+01 +4.605e-01 3.281e-02 3.815e-02 6.245e-02 1.003e+00 1.049e+00 1.013e+00 1.045e+00 1.269e+02 7.832e+01 +4.578e-01 3.311e-02 3.821e-02 5.996e-02 1.003e+00 1.048e+00 1.013e+00 1.046e+00 1.241e+02 7.832e+01 +4.551e-01 3.342e-02 3.827e-02 5.758e-02 1.003e+00 1.047e+00 1.013e+00 1.048e+00 1.213e+02 7.832e+01 +4.523e-01 3.373e-02 3.833e-02 5.531e-02 1.003e+00 1.046e+00 1.013e+00 1.049e+00 1.186e+02 7.832e+01 +4.494e-01 3.404e-02 3.839e-02 5.313e-02 1.003e+00 1.045e+00 1.013e+00 1.051e+00 1.159e+02 7.832e+01 +4.465e-01 3.436e-02 3.844e-02 5.104e-02 1.003e+00 1.044e+00 1.013e+00 1.053e+00 1.134e+02 7.832e+01 +4.436e-01 3.467e-02 3.850e-02 4.903e-02 1.003e+00 1.043e+00 1.013e+00 1.056e+00 1.109e+02 7.832e+01 +4.406e-01 3.499e-02 3.855e-02 4.709e-02 1.003e+00 1.042e+00 1.013e+00 1.058e+00 1.085e+02 7.832e+01 +4.375e-01 3.532e-02 3.860e-02 4.522e-02 1.003e+00 1.041e+00 1.013e+00 1.061e+00 1.061e+02 7.832e+01 +4.343e-01 3.565e-02 3.865e-02 4.342e-02 1.003e+00 1.040e+00 1.013e+00 1.064e+00 1.039e+02 7.832e+01 +4.311e-01 3.598e-02 3.869e-02 4.171e-02 1.003e+00 1.039e+00 1.013e+00 1.067e+00 1.017e+02 7.832e+01 +4.278e-01 3.631e-02 3.873e-02 4.003e-02 1.003e+00 1.038e+00 1.013e+00 1.070e+00 9.960e+01 7.832e+01 +4.245e-01 3.664e-02 3.877e-02 3.837e-02 1.003e+00 1.038e+00 1.013e+00 1.074e+00 9.759e+01 7.832e+01 +4.211e-01 3.698e-02 3.880e-02 3.681e-02 1.003e+00 1.037e+00 1.013e+00 1.078e+00 9.567e+01 7.832e+01 +4.176e-01 3.733e-02 3.883e-02 3.533e-02 1.003e+00 1.036e+00 1.013e+00 1.082e+00 9.383e+01 7.832e+01 +4.140e-01 3.767e-02 3.884e-02 3.393e-02 1.003e+00 1.036e+00 1.013e+00 1.087e+00 9.208e+01 7.832e+01 +4.103e-01 3.802e-02 3.884e-02 3.263e-02 1.003e+00 1.035e+00 1.013e+00 1.091e+00 9.043e+01 7.832e+01 +4.066e-01 3.837e-02 3.882e-02 3.144e-02 1.003e+00 1.035e+00 1.013e+00 1.096e+00 8.888e+01 7.832e+01 +4.029e-01 3.873e-02 3.881e-02 3.041e-02 1.003e+00 1.034e+00 1.013e+00 1.100e+00 8.742e+01 7.832e+01 +3.991e-01 3.908e-02 3.875e-02 2.960e-02 1.004e+00 1.034e+00 1.013e+00 1.103e+00 8.607e+01 7.832e+01 +3.953e-01 3.945e-02 3.867e-02 2.896e-02 1.004e+00 1.033e+00 1.014e+00 1.106e+00 8.483e+01 7.832e+01 +3.915e-01 3.981e-02 3.858e-02 2.847e-02 1.004e+00 1.032e+00 1.014e+00 1.108e+00 8.369e+01 7.832e+01 +3.877e-01 4.018e-02 3.847e-02 2.811e-02 1.004e+00 1.031e+00 1.014e+00 1.110e+00 8.267e+01 7.832e+01 +3.840e-01 4.055e-02 3.836e-02 2.786e-02 1.004e+00 1.030e+00 1.014e+00 1.111e+00 8.175e+01 7.832e+01 +3.803e-01 4.093e-02 3.821e-02 2.771e-02 1.004e+00 1.029e+00 1.015e+00 1.111e+00 8.094e+01 7.832e+01 +3.766e-01 4.131e-02 3.808e-02 2.761e-02 1.004e+00 1.029e+00 1.015e+00 1.112e+00 8.024e+01 7.832e+01 +3.729e-01 4.169e-02 3.796e-02 2.754e-02 1.004e+00 1.028e+00 1.015e+00 1.111e+00 7.965e+01 7.832e+01 +3.692e-01 4.207e-02 3.784e-02 2.751e-02 1.004e+00 1.027e+00 1.016e+00 1.111e+00 7.917e+01 7.832e+01 +3.655e-01 4.246e-02 3.772e-02 2.750e-02 1.004e+00 1.026e+00 1.016e+00 1.111e+00 7.880e+01 7.832e+01 +3.619e-01 4.286e-02 3.762e-02 2.754e-02 1.004e+00 1.025e+00 1.016e+00 1.110e+00 7.853e+01 7.832e+01 +3.583e-01 4.325e-02 3.749e-02 2.756e-02 1.004e+00 1.024e+00 1.017e+00 1.110e+00 7.838e+01 7.832e+01 +3.547e-01 4.365e-02 3.736e-02 2.759e-02 1.004e+00 1.023e+00 1.017e+00 1.109e+00 7.832e+01 7.832e+01 +3.511e-01 4.406e-02 3.724e-02 2.763e-02 1.004e+00 1.022e+00 1.017e+00 1.109e+00 7.838e+01 7.832e+01 +3.475e-01 4.446e-02 3.712e-02 2.766e-02 1.004e+00 1.021e+00 1.018e+00 1.108e+00 7.853e+01 7.832e+01 +3.440e-01 4.488e-02 3.701e-02 2.780e-02 1.005e+00 1.020e+00 1.018e+00 1.107e+00 7.879e+01 7.832e+01 +3.404e-01 4.529e-02 3.691e-02 2.780e-02 1.005e+00 1.019e+00 1.018e+00 1.107e+00 7.915e+01 7.832e+01 +3.369e-01 4.571e-02 3.680e-02 2.789e-02 1.005e+00 1.018e+00 1.019e+00 1.106e+00 7.962e+01 7.832e+01 +3.335e-01 4.613e-02 3.670e-02 2.804e-02 1.005e+00 1.017e+00 1.019e+00 1.105e+00 8.018e+01 7.832e+01 +3.300e-01 4.656e-02 3.660e-02 2.807e-02 1.005e+00 1.016e+00 1.019e+00 1.104e+00 8.084e+01 7.832e+01 +3.265e-01 4.699e-02 3.650e-02 2.815e-02 1.005e+00 1.015e+00 1.020e+00 1.103e+00 8.160e+01 7.832e+01 +3.231e-01 4.742e-02 3.640e-02 2.823e-02 1.005e+00 1.014e+00 1.020e+00 1.103e+00 8.245e+01 7.832e+01 +3.197e-01 4.786e-02 3.630e-02 2.831e-02 1.005e+00 1.014e+00 1.020e+00 1.102e+00 8.340e+01 7.832e+01 +3.163e-01 4.831e-02 3.621e-02 2.838e-02 1.005e+00 1.013e+00 1.020e+00 1.101e+00 8.445e+01 7.832e+01 +3.129e-01 4.875e-02 3.612e-02 2.846e-02 1.005e+00 1.012e+00 1.021e+00 1.101e+00 8.559e+01 7.832e+01 +3.095e-01 4.920e-02 3.603e-02 2.855e-02 1.005e+00 1.011e+00 1.021e+00 1.100e+00 8.681e+01 7.832e+01 +3.062e-01 4.966e-02 3.598e-02 2.859e-02 1.005e+00 1.010e+00 1.021e+00 1.099e+00 8.813e+01 7.832e+01 +3.029e-01 5.012e-02 3.590e-02 2.867e-02 1.005e+00 1.010e+00 1.021e+00 1.098e+00 8.954e+01 7.832e+01 +2.996e-01 5.058e-02 3.580e-02 2.879e-02 1.005e+00 1.009e+00 1.022e+00 1.097e+00 9.103e+01 7.832e+01 +2.964e-01 5.105e-02 3.573e-02 2.887e-02 1.006e+00 1.008e+00 1.022e+00 1.097e+00 9.261e+01 7.832e+01 +2.931e-01 5.152e-02 3.565e-02 2.895e-02 1.006e+00 1.007e+00 1.022e+00 1.096e+00 9.427e+01 7.832e+01 +2.899e-01 5.200e-02 3.560e-02 2.902e-02 1.006e+00 1.007e+00 1.022e+00 1.095e+00 9.602e+01 7.832e+01 +2.867e-01 5.248e-02 3.554e-02 2.909e-02 1.006e+00 1.006e+00 1.023e+00 1.095e+00 9.785e+01 7.832e+01 +2.836e-01 5.297e-02 3.547e-02 2.919e-02 1.006e+00 1.005e+00 1.023e+00 1.094e+00 9.975e+01 7.832e+01 +2.805e-01 5.346e-02 3.541e-02 2.927e-02 1.006e+00 1.005e+00 1.023e+00 1.093e+00 1.017e+02 7.832e+01 +2.774e-01 5.395e-02 3.534e-02 2.933e-02 1.006e+00 1.004e+00 1.023e+00 1.092e+00 1.038e+02 7.832e+01 +2.743e-01 5.445e-02 3.526e-02 2.946e-02 1.006e+00 1.003e+00 1.024e+00 1.092e+00 1.059e+02 7.832e+01 +2.712e-01 5.495e-02 3.521e-02 2.951e-02 1.006e+00 1.003e+00 1.024e+00 1.091e+00 1.081e+02 7.832e+01 +2.682e-01 5.546e-02 3.517e-02 2.957e-02 1.006e+00 1.002e+00 1.024e+00 1.090e+00 1.104e+02 7.832e+01 +2.652e-01 5.598e-02 3.513e-02 2.964e-02 1.006e+00 1.001e+00 1.024e+00 1.090e+00 1.128e+02 7.832e+01 +2.623e-01 5.649e-02 3.506e-02 2.967e-02 1.006e+00 1.001e+00 1.024e+00 1.089e+00 1.152e+02 7.832e+01 +2.594e-01 5.702e-02 3.504e-02 2.978e-02 1.006e+00 1.000e+00 1.025e+00 1.088e+00 1.177e+02 7.832e+01 +2.565e-01 5.754e-02 3.498e-02 2.989e-02 1.006e+00 9.996e-01 1.025e+00 1.088e+00 1.202e+02 7.832e+01 +2.536e-01 5.808e-02 3.495e-02 2.993e-02 1.006e+00 9.990e-01 1.025e+00 1.087e+00 1.228e+02 7.832e+01 +2.508e-01 5.861e-02 3.492e-02 2.996e-02 1.006e+00 9.985e-01 1.025e+00 1.087e+00 1.255e+02 7.832e+01 +2.479e-01 5.916e-02 3.486e-02 3.007e-02 1.007e+00 9.979e-01 1.025e+00 1.086e+00 1.282e+02 7.832e+01 +2.451e-01 5.970e-02 3.483e-02 3.012e-02 1.007e+00 9.974e-01 1.026e+00 1.085e+00 1.310e+02 7.832e+01 +2.424e-01 6.026e-02 3.481e-02 3.017e-02 1.007e+00 9.968e-01 1.026e+00 1.085e+00 1.338e+02 7.832e+01 +2.397e-01 6.081e-02 3.476e-02 3.025e-02 1.007e+00 9.963e-01 1.026e+00 1.084e+00 1.367e+02 7.832e+01 +2.370e-01 6.138e-02 3.478e-02 3.027e-02 1.007e+00 9.958e-01 1.026e+00 1.084e+00 1.397e+02 7.832e+01 +2.343e-01 6.194e-02 3.473e-02 3.035e-02 1.007e+00 9.953e-01 1.026e+00 1.083e+00 1.427e+02 7.832e+01 +2.317e-01 6.252e-02 3.470e-02 3.041e-02 1.007e+00 9.948e-01 1.026e+00 1.082e+00 1.457e+02 7.832e+01 +2.292e-01 6.310e-02 3.469e-02 3.046e-02 1.007e+00 9.944e-01 1.026e+00 1.082e+00 1.488e+02 7.832e+01 +2.265e-01 6.368e-02 3.467e-02 3.049e-02 1.007e+00 9.939e-01 1.026e+00 1.081e+00 1.520e+02 7.832e+01 +2.240e-01 6.427e-02 3.465e-02 3.057e-02 1.007e+00 9.935e-01 1.027e+00 1.081e+00 1.551e+02 7.832e+01 +2.215e-01 6.486e-02 3.465e-02 3.060e-02 1.007e+00 9.930e-01 1.027e+00 1.080e+00 1.584e+02 7.832e+01 +2.190e-01 6.546e-02 3.463e-02 3.065e-02 1.007e+00 9.926e-01 1.027e+00 1.080e+00 1.616e+02 7.832e+01 +2.167e-01 6.607e-02 3.463e-02 3.067e-02 1.007e+00 9.922e-01 1.027e+00 1.079e+00 1.649e+02 7.832e+01 +2.142e-01 6.668e-02 3.461e-02 3.076e-02 1.007e+00 9.917e-01 1.027e+00 1.079e+00 1.683e+02 7.832e+01 +2.118e-01 6.730e-02 3.460e-02 3.081e-02 1.007e+00 9.913e-01 1.027e+00 1.078e+00 1.717e+02 7.832e+01 +2.095e-01 6.792e-02 3.460e-02 3.086e-02 1.007e+00 9.909e-01 1.027e+00 1.078e+00 1.751e+02 7.832e+01 +2.072e-01 6.855e-02 3.459e-02 3.092e-02 1.008e+00 9.905e-01 1.027e+00 1.077e+00 1.785e+02 7.832e+01 +2.049e-01 6.918e-02 3.460e-02 3.092e-02 1.008e+00 9.902e-01 1.027e+00 1.077e+00 1.820e+02 7.832e+01 diff --git a/results/fitnup107nup358/modelNup107_con.png b/results/fitnup107nup358/modelNup107_con.png new file mode 100644 index 0000000000000000000000000000000000000000..09123d8bf5503d59154a0aa300ce4f819665369c Binary files /dev/null and b/results/fitnup107nup358/modelNup107_con.png differ diff --git a/results/fitnup107nup358/modelNup107_con_10.png b/results/fitnup107nup358/modelNup107_con_10.png new file mode 100644 index 0000000000000000000000000000000000000000..9f786d74a5bd9554c666e59c156ffdfb2621f377 Binary files /dev/null and b/results/fitnup107nup358/modelNup107_con_10.png differ diff --git a/results/fitnup107nup358/modelNup107_uncon.png b/results/fitnup107nup358/modelNup107_uncon.png new file mode 100644 index 0000000000000000000000000000000000000000..95c46d82681f07e7f7a9b3baf6664e6cf7cfbf1b Binary files /dev/null and b/results/fitnup107nup358/modelNup107_uncon.png differ diff --git a/results/fitnup107nup358/modelNup358_con.png b/results/fitnup107nup358/modelNup358_con.png new file mode 100644 index 0000000000000000000000000000000000000000..4f9ec808f2e667649919ee257e861603a7be267a Binary files /dev/null and b/results/fitnup107nup358/modelNup358_con.png differ diff --git a/results/fitnup107nup358/modelNup358_con_10.png b/results/fitnup107nup358/modelNup358_con_10.png new file mode 100644 index 0000000000000000000000000000000000000000..774b40d4ac23d61ecf3785251c2c8581f2e9bcb6 Binary files /dev/null and b/results/fitnup107nup358/modelNup358_con_10.png differ diff --git a/results/fitnup107nup358/modelNup358_uncon.png b/results/fitnup107nup358/modelNup358_uncon.png new file mode 100644 index 0000000000000000000000000000000000000000..62caca0d87bb22b7c8ff7494d5f826a26edbdad9 Binary files /dev/null and b/results/fitnup107nup358/modelNup358_uncon.png differ diff --git a/results/fitnup107nup358/plotResults.R b/results/fitnup107nup358/plotResults.R new file mode 100644 index 0000000000000000000000000000000000000000..4ffc1651da72fa6c10e24d0d1fd3debd796d4757 --- /dev/null +++ b/results/fitnup107nup358/plotResults.R @@ -0,0 +1,112 @@ +library("Hmisc") + +# ----- +setwd('/Users/toni/Dropbox/Shotaro_PoreMaturation/Simulations/LiveCellImaging') +simulNUP107<-read.csv('simulated_NUP107.csv', skip=3, sep='\t') +simulNUP358<-read.csv('simulated_NUP358.csv', skip=3, sep='\t') +datain <-read.csv('/Users/toni/Dropbox/Shotaro_PoreMaturation/ExpData/LiveCellImaging/HeLa4D.csv', sep='\t') + +# ----- +plotDataF <-function(d, simul, names, cols, alpha, fname) { + # function to plot data no separation between interphase and postmitotic + pdf(fname, width = 2.8, height=2, useDingbats = F,pointsize = 8, paper="special", fonts="Helvetica") + lw<- 2 + par(mar=c(4, 4, 1,1), mgp=c(2.5,1,0)) + plot(d$x, d$y, xlab = "", ylab = "Normalized total intensity", + type = "n", yaxt = 'n', xaxt = 'n', xlim = c(0, 120), ylim = c(0, 1.4), bty='l') + axis(1, tck = -0.05, at=c(0,20,40,60,80,100,120)) + axis(2, las=2, tck =-0.05, at=c(0,0.2,0.4,0.6,0.8,1,1.2)) + mtext(side=1, text="Time after AO (min)", line=2) + cols1_t <-adjustcolor(cols[1], alpha.f = alpha) + polygon(c(d$x, rev(d$x)), c(d$upper, rev(d$lower)), col = cols1_t, border = NA) + points(d$x, d$y, pch = 17, col = cols[1], cex=0.6) + colfun <- colorRampPalette(c(cols[1], 'gray60')) + lines(simul[,names[1]], simul[,names[2]], lwd = lw, col = colfun(50)[40]) + lines(simul[,names[1]], simul[,names[3]], lwd = lw, lty = 2, col = cols[2]) + lines(simul[,names[1]], simul[,names[4]], lwd = lw, col=cols[2]) + par(xpd=TRUE) + legend(0,1.65, c("data", "post-mitotic+interphase", "post-mitotic", "interphase"), + lwd = c(0,1.5,1.5,1.5),pch = c(17,-1,-1,-1), lty = c(0,1,2,1), + col=c(cols[1], colfun(50)[40], cols[2], cols[2]), horiz = F, bty='n', + y.intersp = 0.88, cex = 0.88) + dev.off() +} + + +# ----- +plotDataAll <-function(d1,d2, simul, names, cols, alpha, fname) { + # function to plot data no separation between interphase and postmitotic + pdf(fname, width = 2.8, height=2, useDingbats = F,pointsize = 8, paper="special", fonts="Helvetica") + lw<- 2 + par(mar=c(4, 4, 1,1), mgp=c(2.5,1,0)) + plot(d1$x, d1$y, xlab = "", ylab = "Normalized total intensity", + type = "n", yaxt = 'n', xaxt = 'n', xlim = c(0, 120), ylim = c(0, 1.4), bty='l') + #title('NUP107 inner-core region') + axis(1, tck = -0.05, at=c(0,20,40,60,80,100,120)) + axis(2, las=2, tck =-0.05, at=c(0,0.2,0.4,0.6,0.8,1,1.2)) + mtext(side=1, text="Time after AO (min)", line=2) + #transparent colors + cols1_t <-adjustcolor(cols[1], alpha.f = alpha) + cols2_t <-adjustcolor(cols[2], alpha.f = alpha) + #plot data for Non-core and core region + polygon(c(d1$x, rev(d1$x)), c(d1$upper, rev(d1$lower)), col = cols1_t, border = NA) + polygon(c(d1$x, rev(d1$x)), c(d2$upper, rev(d2$lower)), col = cols2_t, border = NA) + points(d1$x, d1$y, pch = 17, col = cols[1], cex=0.6) + points(d1$x, d2$y, pch = 17, col = cols[2], cex=0.6) + + #plot simulations + lines(simul[,names[1]], simul[,names[2]], lwd = lw, col = cols[1]) + lines(simul[,names[1]], simul[,names[3]], lwd = lw, col = cols[2]) + par(xpd=TRUE) + legend(0,1.65, c("Non-core", "Inner-core","Non-core model", "Inner-core model"), + lwd = c(0,0, 1.5,1.5), pch = c(17,17, -1,-1), lty = c(0,0, 1,1), + col= c(cols[1], cols[2], cols[1], cols[2]), horiz = F, bty='n', + y.intersp = 0.88, cex = 0.88) + + dev.off() +} + +# ----- +plotData_IP <-function(simul1, simul2, names, cols, fname) { + pdf(fname, width = 2.8, height=2, useDingbats = F,pointsize = 8, paper="special", fonts="Helvetica") + lw<- 2 + par(mar=c(4, 4, 1,1), mgp=c(2.5,1,0)) + plot(simul1[,names[1]], simul1[,names[2]], xlab = "", ylab = "Normalized interphase density", + type = "n", yaxt = 'n', xaxt = 'n', xlim = c(0, 120), ylim = c(0, 1.1), bty='l') + axis(1, tck = -0.05, at=c(0,20,40,60,80,100,120)) + axis(2, las=2, tck =-0.05, at=c(0,0.2,0.4,0.6,0.8,1)) + mtext(side=1, text="Time after AO (min)", line=2) + lines(simul1[,names[1]], simul1[,names[2]], lwd = lw, lty = 1, col = cols[1]) + lines(simul2[,names[1]], simul2[,names[2]], lwd = lw, lty = 1, col = cols[2]) + par(xpd=TRUE) + legend(0,1.1, c("Nup107", "Nup358"), + lwd = c(1.5,1.5), + col=c(cols[1], cols[2]), horiz = F, bty='n', + y.intersp = 0.88, cex = 0.88) + + dev.off() +} + + +# ---- +nc107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_NC, 'upper'=datain$NUP107_NC + datain$NUP107_NC_std, + 'lower' = datain$NUP107_NC - datain$NUP107_NC_std) +c107 <- data.frame('x'= datain$time, 'y' = datain$NUP107_C, 'upper'=datain$NUP107_C + datain$NUP107_C_std, + 'lower' = datain$NUP107_C - datain$NUP107_C_std) + +nc358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_NC, 'upper'=datain$NUP358_NC + datain$NUP358_NC_std, + 'lower' = datain$NUP358_NC - datain$NUP358_NC_std) +c358 <- data.frame('x'= datain$time, 'y' = datain$NUP358_C, 'upper'=datain$NUP358_C + datain$NUP358_C_std, + 'lower' = datain$NUP358_C - datain$NUP358_C_std) + +plotDataAll(nc107, c107, simulNUP107, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP107.pdf') +plotDataAll(nc358, c358, simulNUP358, c('time','NC', 'C'), c('#20BCEF', '#A2C964' ), 0.4, 'NUP358.pdf') + +plotDataF(nc107, simulNUP107, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray20'), 0.4, 'NUP107_NC.pdf') +plotDataF(c107, simulNUP107, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray20'), 0.4, 'NUP107_C.pdf') + +plotDataF(nc358, simulNUP358, c('time','NC', 'NC_pm', 'NC_ip'), c('#20BCEF', 'gray20' ), 0.4, 'NUP358_NC.pdf') +plotDataF(c358, simulNUP358, c('time','C', 'C_pm', 'C_ip'), c('#A2C964', 'gray20'), 0.4, 'NUP358_C.pdf') + +plotData_IP(simulNUP107, simulNUP358, c('time','C_ip_density'), c('#A2C964' , '#d01c8b'), 'NUP107_358_C_den.pdf' ) +#c('#A2C964', '#C1FF00') \ No newline at end of file diff --git a/results/fitnup107nup358/simulated_NUP107.csv b/results/fitnup107nup358/simulated_NUP107.csv new file mode 100644 index 0000000000000000000000000000000000000000..a250d641381649797286f58e07eff2a4c3bf2a08 --- /dev/null +++ b/results/fitnup107nup358/simulated_NUP107.csv @@ -0,0 +1,126 @@ +% Parameters k1_pm k2_pm k1_ip k2_ip NC_frac_pm C_frac_pm factProd_NC_p1 factPro_NC_p2 factProd_C_p1 factPro_C_p2 (only the first 4 parameters are fixed) +0.3547 0.0437 0.0374 0.0276 0.9200 0.5000 1.0043 1.0226 1.0170 1.1092 +tpm = [4 4], tip = [10 10], N1 = 1, N2 = 1 +time NC NC_pm NC_ip C C_pm C_ip C_ip_pure C_ip_density +4 0.000e+00 0.000e+00 0.000e+00 0.000e+00 0.000e+00 0.000e+00 0.000e+00 0.000e+00 +5 2.759e-01 2.759e-01 0.000e+00 1.519e-01 1.519e-01 0.000e+00 0.000e+00 0.000e+00 +6 4.695e-01 4.695e-01 0.000e+00 2.584e-01 2.584e-01 0.000e+00 0.000e+00 0.000e+00 +7 6.052e-01 6.052e-01 0.000e+00 3.331e-01 3.331e-01 0.000e+00 0.000e+00 0.000e+00 +8 7.004e-01 7.004e-01 0.000e+00 3.855e-01 3.855e-01 0.000e+00 0.000e+00 0.000e+00 +9 7.671e-01 7.671e-01 0.000e+00 4.222e-01 4.222e-01 0.000e+00 0.000e+00 0.000e+00 +10 8.140e-01 8.140e-01 0.000e+00 4.480e-01 4.480e-01 0.000e+00 0.000e+00 0.000e+00 +11 8.497e-01 8.468e-01 2.946e-03 4.847e-01 4.660e-01 1.864e-02 3.667e-02 6.382e-02 +12 8.756e-01 8.699e-01 5.784e-03 5.153e-01 4.787e-01 3.661e-02 7.199e-02 1.159e-01 +13 8.945e-01 8.860e-01 8.517e-03 5.415e-01 4.876e-01 5.391e-02 1.060e-01 1.606e-01 +14 9.085e-01 8.973e-01 1.115e-02 5.644e-01 4.939e-01 7.058e-02 1.388e-01 2.006e-01 +15 9.190e-01 9.053e-01 1.369e-02 5.849e-01 4.982e-01 8.663e-02 1.704e-01 2.379e-01 +16 9.270e-01 9.109e-01 1.613e-02 6.034e-01 5.013e-01 1.021e-01 2.008e-01 2.723e-01 +17 9.333e-01 9.148e-01 1.849e-02 6.205e-01 5.035e-01 1.170e-01 2.301e-01 3.040e-01 +18 9.383e-01 9.175e-01 2.075e-02 6.363e-01 5.050e-01 1.314e-01 2.583e-01 3.334e-01 +19 9.424e-01 9.194e-01 2.294e-02 6.512e-01 5.060e-01 1.452e-01 2.855e-01 3.608e-01 +20 9.458e-01 9.208e-01 2.504e-02 6.653e-01 5.068e-01 1.585e-01 3.117e-01 3.866e-01 +21 9.488e-01 9.217e-01 2.707e-02 6.786e-01 5.073e-01 1.713e-01 3.369e-01 4.111e-01 +22 9.514e-01 9.224e-01 2.902e-02 6.914e-01 5.077e-01 1.837e-01 3.613e-01 4.345e-01 +23 9.538e-01 9.229e-01 3.091e-02 7.035e-01 5.079e-01 1.956e-01 3.847e-01 4.569e-01 +24 9.559e-01 9.232e-01 3.272e-02 7.152e-01 5.081e-01 2.071e-01 4.072e-01 4.780e-01 +25 9.579e-01 9.234e-01 3.447e-02 7.264e-01 5.082e-01 2.181e-01 4.290e-01 4.980e-01 +26 9.597e-01 9.236e-01 3.615e-02 7.371e-01 5.083e-01 2.288e-01 4.499e-01 5.170e-01 +27 9.615e-01 9.237e-01 3.777e-02 7.474e-01 5.084e-01 2.390e-01 4.701e-01 5.352e-01 +28 9.631e-01 9.238e-01 3.933e-02 7.573e-01 5.084e-01 2.489e-01 4.895e-01 5.525e-01 +29 9.647e-01 9.238e-01 4.083e-02 7.669e-01 5.084e-01 2.584e-01 5.082e-01 5.691e-01 +30 9.662e-01 9.239e-01 4.228e-02 7.761e-01 5.085e-01 2.676e-01 5.263e-01 5.850e-01 +31 9.676e-01 9.239e-01 4.368e-02 7.849e-01 5.085e-01 2.764e-01 5.436e-01 6.004e-01 +32 9.689e-01 9.239e-01 4.502e-02 7.934e-01 5.085e-01 2.849e-01 5.604e-01 6.152e-01 +33 9.702e-01 9.239e-01 4.632e-02 8.016e-01 5.085e-01 2.931e-01 5.765e-01 6.296e-01 +34 9.715e-01 9.239e-01 4.756e-02 8.095e-01 5.085e-01 3.010e-01 5.920e-01 6.437e-01 +35 9.727e-01 9.239e-01 4.877e-02 8.172e-01 5.085e-01 3.087e-01 6.070e-01 6.574e-01 +36 9.739e-01 9.239e-01 4.992e-02 8.245e-01 5.085e-01 3.160e-01 6.214e-01 6.708e-01 +37 9.750e-01 9.240e-01 5.104e-02 8.316e-01 5.085e-01 3.230e-01 6.353e-01 6.840e-01 +38 9.761e-01 9.240e-01 5.211e-02 8.384e-01 5.085e-01 3.298e-01 6.486e-01 6.970e-01 +39 9.771e-01 9.240e-01 5.315e-02 8.449e-01 5.085e-01 3.364e-01 6.615e-01 7.098e-01 +40 9.781e-01 9.240e-01 5.415e-02 8.512e-01 5.085e-01 3.427e-01 6.739e-01 7.220e-01 +41 9.791e-01 9.240e-01 5.511e-02 8.573e-01 5.085e-01 3.488e-01 6.859e-01 7.338e-01 +42 9.800e-01 9.240e-01 5.603e-02 8.632e-01 5.085e-01 3.546e-01 6.974e-01 7.450e-01 +43 9.809e-01 9.240e-01 5.692e-02 8.688e-01 5.085e-01 3.603e-01 7.085e-01 7.558e-01 +44 9.817e-01 9.240e-01 5.778e-02 8.742e-01 5.085e-01 3.657e-01 7.192e-01 7.662e-01 +45 9.826e-01 9.240e-01 5.861e-02 8.795e-01 5.085e-01 3.709e-01 7.295e-01 7.761e-01 +46 9.834e-01 9.240e-01 5.941e-02 8.845e-01 5.085e-01 3.760e-01 7.394e-01 7.857e-01 +47 9.841e-01 9.240e-01 6.017e-02 8.894e-01 5.085e-01 3.809e-01 7.490e-01 7.948e-01 +48 9.849e-01 9.240e-01 6.091e-02 8.940e-01 5.085e-01 3.855e-01 7.582e-01 8.036e-01 +49 9.856e-01 9.240e-01 6.163e-02 8.986e-01 5.085e-01 3.900e-01 7.670e-01 8.121e-01 +50 9.863e-01 9.240e-01 6.231e-02 9.029e-01 5.085e-01 3.944e-01 7.756e-01 8.202e-01 +51 9.869e-01 9.240e-01 6.297e-02 9.071e-01 5.085e-01 3.986e-01 7.838e-01 8.279e-01 +52 9.876e-01 9.240e-01 6.361e-02 9.111e-01 5.085e-01 4.026e-01 7.917e-01 8.354e-01 +53 9.882e-01 9.240e-01 6.422e-02 9.150e-01 5.085e-01 4.065e-01 7.994e-01 8.426e-01 +54 9.888e-01 9.240e-01 6.482e-02 9.187e-01 5.085e-01 4.102e-01 8.067e-01 8.495e-01 +55 9.893e-01 9.240e-01 6.538e-02 9.223e-01 5.085e-01 4.138e-01 8.138e-01 8.561e-01 +56 9.899e-01 9.240e-01 6.593e-02 9.258e-01 5.085e-01 4.173e-01 8.206e-01 8.624e-01 +57 9.904e-01 9.240e-01 6.646e-02 9.292e-01 5.085e-01 4.206e-01 8.272e-01 8.685e-01 +58 9.909e-01 9.240e-01 6.697e-02 9.324e-01 5.085e-01 4.239e-01 8.335e-01 8.744e-01 +59 9.914e-01 9.240e-01 6.746e-02 9.355e-01 5.085e-01 4.270e-01 8.396e-01 8.801e-01 +60 9.919e-01 9.240e-01 6.793e-02 9.385e-01 5.085e-01 4.300e-01 8.455e-01 8.855e-01 +61 9.923e-01 9.240e-01 6.839e-02 9.414e-01 5.085e-01 4.328e-01 8.512e-01 8.907e-01 +62 9.928e-01 9.240e-01 6.883e-02 9.441e-01 5.085e-01 4.356e-01 8.566e-01 8.958e-01 +63 9.932e-01 9.240e-01 6.925e-02 9.468e-01 5.085e-01 4.383e-01 8.619e-01 9.006e-01 +64 9.936e-01 9.240e-01 6.966e-02 9.494e-01 5.085e-01 4.409e-01 8.670e-01 9.053e-01 +65 9.940e-01 9.240e-01 7.005e-02 9.519e-01 5.085e-01 4.433e-01 8.718e-01 9.097e-01 +66 9.944e-01 9.240e-01 7.043e-02 9.542e-01 5.085e-01 4.457e-01 8.765e-01 9.140e-01 +67 9.947e-01 9.240e-01 7.079e-02 9.565e-01 5.085e-01 4.480e-01 8.811e-01 9.181e-01 +68 9.951e-01 9.240e-01 7.114e-02 9.588e-01 5.085e-01 4.503e-01 8.854e-01 9.220e-01 +69 9.954e-01 9.240e-01 7.148e-02 9.609e-01 5.085e-01 4.524e-01 8.896e-01 9.258e-01 +70 9.958e-01 9.240e-01 7.180e-02 9.630e-01 5.085e-01 4.544e-01 8.937e-01 9.293e-01 +71 9.961e-01 9.240e-01 7.211e-02 9.649e-01 5.085e-01 4.564e-01 8.976e-01 9.328e-01 +72 9.964e-01 9.240e-01 7.242e-02 9.669e-01 5.085e-01 4.583e-01 9.013e-01 9.360e-01 +73 9.967e-01 9.240e-01 7.271e-02 9.687e-01 5.085e-01 4.602e-01 9.049e-01 9.391e-01 +74 9.969e-01 9.240e-01 7.299e-02 9.705e-01 5.085e-01 4.620e-01 9.084e-01 9.421e-01 +75 9.972e-01 9.240e-01 7.326e-02 9.722e-01 5.085e-01 4.637e-01 9.118e-01 9.449e-01 +76 9.975e-01 9.240e-01 7.352e-02 9.738e-01 5.085e-01 4.653e-01 9.150e-01 9.476e-01 +77 9.977e-01 9.240e-01 7.377e-02 9.754e-01 5.085e-01 4.669e-01 9.181e-01 9.502e-01 +78 9.980e-01 9.240e-01 7.401e-02 9.769e-01 5.085e-01 4.684e-01 9.211e-01 9.526e-01 +79 9.982e-01 9.240e-01 7.424e-02 9.784e-01 5.085e-01 4.699e-01 9.240e-01 9.550e-01 +80 9.984e-01 9.240e-01 7.446e-02 9.798e-01 5.085e-01 4.713e-01 9.268e-01 9.572e-01 +81 9.986e-01 9.240e-01 7.468e-02 9.812e-01 5.085e-01 4.727e-01 9.295e-01 9.593e-01 +82 9.988e-01 9.240e-01 7.489e-02 9.825e-01 5.085e-01 4.740e-01 9.321e-01 9.613e-01 +83 9.990e-01 9.240e-01 7.509e-02 9.838e-01 5.085e-01 4.752e-01 9.346e-01 9.632e-01 +84 9.992e-01 9.240e-01 7.528e-02 9.850e-01 5.085e-01 4.765e-01 9.370e-01 9.650e-01 +85 9.994e-01 9.240e-01 7.547e-02 9.862e-01 5.085e-01 4.776e-01 9.393e-01 9.667e-01 +86 9.996e-01 9.240e-01 7.565e-02 9.873e-01 5.085e-01 4.788e-01 9.415e-01 9.683e-01 +87 9.998e-01 9.240e-01 7.582e-02 9.884e-01 5.085e-01 4.799e-01 9.437e-01 9.698e-01 +88 9.999e-01 9.240e-01 7.598e-02 9.894e-01 5.085e-01 4.809e-01 9.457e-01 9.712e-01 +89 1.000e+00 9.240e-01 7.614e-02 9.904e-01 5.085e-01 4.819e-01 9.477e-01 9.726e-01 +90 1.000e+00 9.240e-01 7.630e-02 9.914e-01 5.085e-01 4.829e-01 9.496e-01 9.739e-01 +91 1.000e+00 9.240e-01 7.645e-02 9.924e-01 5.085e-01 4.838e-01 9.515e-01 9.751e-01 +92 1.001e+00 9.240e-01 7.659e-02 9.933e-01 5.085e-01 4.847e-01 9.533e-01 9.762e-01 +93 1.001e+00 9.240e-01 7.673e-02 9.941e-01 5.085e-01 4.856e-01 9.550e-01 9.773e-01 +94 1.001e+00 9.240e-01 7.686e-02 9.950e-01 5.085e-01 4.865e-01 9.566e-01 9.783e-01 +95 1.001e+00 9.240e-01 7.699e-02 9.958e-01 5.085e-01 4.873e-01 9.582e-01 9.792e-01 +96 1.001e+00 9.240e-01 7.711e-02 9.966e-01 5.085e-01 4.880e-01 9.597e-01 9.801e-01 +97 1.001e+00 9.240e-01 7.723e-02 9.973e-01 5.085e-01 4.888e-01 9.612e-01 9.809e-01 +98 1.001e+00 9.240e-01 7.734e-02 9.980e-01 5.085e-01 4.895e-01 9.626e-01 9.816e-01 +99 1.001e+00 9.240e-01 7.745e-02 9.987e-01 5.085e-01 4.902e-01 9.640e-01 9.823e-01 +100 1.002e+00 9.240e-01 7.756e-02 9.994e-01 5.085e-01 4.909e-01 9.653e-01 9.830e-01 +101 1.002e+00 9.240e-01 7.766e-02 1.000e+00 5.085e-01 4.915e-01 9.666e-01 9.836e-01 +102 1.002e+00 9.240e-01 7.776e-02 1.001e+00 5.085e-01 4.922e-01 9.678e-01 9.841e-01 +103 1.002e+00 9.240e-01 7.785e-02 1.001e+00 5.085e-01 4.928e-01 9.690e-01 9.846e-01 +104 1.002e+00 9.240e-01 7.795e-02 1.002e+00 5.085e-01 4.933e-01 9.701e-01 9.851e-01 +105 1.002e+00 9.240e-01 7.803e-02 1.002e+00 5.085e-01 4.939e-01 9.712e-01 9.855e-01 +106 1.002e+00 9.240e-01 7.812e-02 1.003e+00 5.085e-01 4.944e-01 9.723e-01 9.858e-01 +107 1.002e+00 9.240e-01 7.820e-02 1.003e+00 5.085e-01 4.949e-01 9.733e-01 9.861e-01 +108 1.002e+00 9.240e-01 7.828e-02 1.004e+00 5.085e-01 4.954e-01 9.743e-01 9.864e-01 +109 1.002e+00 9.240e-01 7.835e-02 1.004e+00 5.085e-01 4.959e-01 9.752e-01 9.867e-01 +110 1.002e+00 9.240e-01 7.843e-02 1.005e+00 5.085e-01 4.964e-01 9.761e-01 9.869e-01 +111 1.002e+00 9.240e-01 7.850e-02 1.005e+00 5.085e-01 4.968e-01 9.770e-01 9.870e-01 +112 1.003e+00 9.240e-01 7.857e-02 1.006e+00 5.085e-01 4.973e-01 9.779e-01 9.872e-01 +113 1.003e+00 9.240e-01 7.863e-02 1.006e+00 5.085e-01 4.977e-01 9.787e-01 9.873e-01 +114 1.003e+00 9.240e-01 7.869e-02 1.007e+00 5.085e-01 4.981e-01 9.794e-01 9.874e-01 +115 1.003e+00 9.240e-01 7.875e-02 1.007e+00 5.085e-01 4.984e-01 9.802e-01 9.874e-01 +116 1.003e+00 9.240e-01 7.881e-02 1.007e+00 5.085e-01 4.988e-01 9.809e-01 9.874e-01 +117 1.003e+00 9.240e-01 7.887e-02 1.008e+00 5.085e-01 4.992e-01 9.816e-01 9.874e-01 +118 1.003e+00 9.240e-01 7.892e-02 1.008e+00 5.085e-01 4.995e-01 9.823e-01 9.874e-01 +119 1.003e+00 9.240e-01 7.897e-02 1.008e+00 5.085e-01 4.998e-01 9.829e-01 9.873e-01 +120 1.003e+00 9.240e-01 7.902e-02 1.009e+00 5.085e-01 5.002e-01 9.836e-01 9.872e-01 +121 1.003e+00 9.240e-01 7.907e-02 1.009e+00 5.085e-01 5.005e-01 9.842e-01 9.871e-01 +122 1.003e+00 9.240e-01 7.912e-02 1.009e+00 5.085e-01 5.008e-01 9.848e-01 9.869e-01 +123 1.003e+00 9.240e-01 7.916e-02 1.010e+00 5.085e-01 5.010e-01 9.853e-01 9.868e-01 +124 1.003e+00 9.240e-01 7.921e-02 1.010e+00 5.085e-01 5.013e-01 9.859e-01 9.866e-01 +125 1.003e+00 9.240e-01 7.925e-02 1.010e+00 5.085e-01 5.016e-01 9.864e-01 9.864e-01 diff --git a/results/fitnup107nup358/simulated_NUP358.csv b/results/fitnup107nup358/simulated_NUP358.csv new file mode 100644 index 0000000000000000000000000000000000000000..95186ef693ca66aeb914a141cfeca241dffbcd58 --- /dev/null +++ b/results/fitnup107nup358/simulated_NUP358.csv @@ -0,0 +1,126 @@ +% Parameters k1_pm k2_pm k1_ip k2_ip NC_frac_pm C_frac_pm factProd_NC_p1 factPro_NC_p2 factProd_C_p1 factPro_C_p2 (only the first 4 parameters are fixed) +0.3547 0.0437 0.0374 0.0276 0.9200 0.5000 1.0043 1.0226 1.0170 1.1092 +tpm = [4 4], tip = [10 10], N1 = 1, N2 = 1 +time NC NC_pm NC_ip C C_pm C_ip C_ip_pure C_ip_density +4 0.000e+00 0.000e+00 0.000e+00 0.000e+00 0.000e+00 0.000e+00 0.000e+00 0.000e+00 +5 6.397e-03 6.397e-03 0.000e+00 3.771e-03 3.771e-03 0.000e+00 0.000e+00 0.000e+00 +6 2.261e-02 2.261e-02 0.000e+00 1.333e-02 1.333e-02 0.000e+00 0.000e+00 0.000e+00 +7 4.521e-02 4.521e-02 0.000e+00 2.665e-02 2.665e-02 0.000e+00 0.000e+00 0.000e+00 +8 7.180e-02 7.180e-02 0.000e+00 4.233e-02 4.233e-02 0.000e+00 0.000e+00 0.000e+00 +9 1.007e-01 1.007e-01 0.000e+00 5.939e-02 5.939e-02 0.000e+00 0.000e+00 0.000e+00 +10 1.309e-01 1.309e-01 0.000e+00 7.716e-02 7.716e-02 0.000e+00 0.000e+00 0.000e+00 +11 1.615e-01 1.615e-01 4.130e-05 9.547e-02 9.519e-02 2.800e-04 5.049e-04 8.789e-04 +12 1.921e-01 1.919e-01 1.617e-04 1.142e-01 1.131e-01 1.096e-03 1.976e-03 3.183e-03 +13 2.223e-01 2.219e-01 3.561e-04 1.332e-01 1.308e-01 2.414e-03 4.352e-03 6.593e-03 +14 2.519e-01 2.513e-01 6.196e-04 1.523e-01 1.481e-01 4.200e-03 7.574e-03 1.095e-02 +15 2.807e-01 2.797e-01 9.476e-04 1.713e-01 1.649e-01 6.424e-03 1.158e-02 1.617e-02 +16 3.086e-01 3.073e-01 1.336e-03 1.902e-01 1.811e-01 9.055e-03 1.633e-02 2.214e-02 +17 3.356e-01 3.339e-01 1.780e-03 2.089e-01 1.968e-01 1.207e-02 2.176e-02 2.874e-02 +18 3.617e-01 3.595e-01 2.276e-03 2.273e-01 2.119e-01 1.543e-02 2.782e-02 3.590e-02 +19 3.869e-01 3.841e-01 2.820e-03 2.455e-01 2.264e-01 1.912e-02 3.448e-02 4.356e-02 +20 4.111e-01 4.077e-01 3.409e-03 2.634e-01 2.403e-01 2.311e-02 4.168e-02 5.169e-02 +21 4.344e-01 4.303e-01 4.040e-03 2.811e-01 2.537e-01 2.739e-02 4.938e-02 6.025e-02 +22 4.568e-01 4.520e-01 4.708e-03 2.984e-01 2.665e-01 3.192e-02 5.755e-02 6.922e-02 +23 4.783e-01 4.729e-01 5.411e-03 3.154e-01 2.788e-01 3.668e-02 6.615e-02 7.856e-02 +24 4.990e-01 4.928e-01 6.146e-03 3.322e-01 2.905e-01 4.167e-02 7.513e-02 8.819e-02 +25 5.188e-01 5.119e-01 6.911e-03 3.486e-01 3.018e-01 4.685e-02 8.448e-02 9.807e-02 +26 5.379e-01 5.302e-01 7.702e-03 3.648e-01 3.126e-01 5.221e-02 9.415e-02 1.082e-01 +27 5.563e-01 5.477e-01 8.518e-03 3.806e-01 3.229e-01 5.774e-02 1.041e-01 1.185e-01 +28 5.739e-01 5.645e-01 9.355e-03 3.962e-01 3.328e-01 6.342e-02 1.144e-01 1.291e-01 +29 5.908e-01 5.806e-01 1.021e-02 4.115e-01 3.423e-01 6.923e-02 1.248e-01 1.398e-01 +30 6.071e-01 5.960e-01 1.109e-02 4.265e-01 3.513e-01 7.517e-02 1.355e-01 1.507e-01 +31 6.227e-01 6.107e-01 1.198e-02 4.412e-01 3.600e-01 8.121e-02 1.464e-01 1.617e-01 +32 6.377e-01 6.248e-01 1.288e-02 4.556e-01 3.683e-01 8.734e-02 1.575e-01 1.729e-01 +33 6.521e-01 6.383e-01 1.380e-02 4.698e-01 3.763e-01 9.355e-02 1.687e-01 1.842e-01 +34 6.659e-01 6.512e-01 1.473e-02 4.837e-01 3.839e-01 9.983e-02 1.800e-01 1.957e-01 +35 6.792e-01 6.636e-01 1.566e-02 4.973e-01 3.912e-01 1.062e-01 1.914e-01 2.073e-01 +36 6.920e-01 6.754e-01 1.660e-02 5.107e-01 3.982e-01 1.126e-01 2.029e-01 2.191e-01 +37 7.043e-01 6.867e-01 1.755e-02 5.238e-01 4.048e-01 1.190e-01 2.145e-01 2.310e-01 +38 7.161e-01 6.976e-01 1.850e-02 5.367e-01 4.112e-01 1.254e-01 2.262e-01 2.431e-01 +39 7.274e-01 7.080e-01 1.946e-02 5.493e-01 4.174e-01 1.319e-01 2.379e-01 2.552e-01 +40 7.383e-01 7.179e-01 2.042e-02 5.616e-01 4.232e-01 1.384e-01 2.496e-01 2.674e-01 +41 7.488e-01 7.274e-01 2.138e-02 5.737e-01 4.288e-01 1.449e-01 2.613e-01 2.795e-01 +42 7.589e-01 7.366e-01 2.234e-02 5.856e-01 4.342e-01 1.514e-01 2.730e-01 2.917e-01 +43 7.686e-01 7.453e-01 2.330e-02 5.973e-01 4.393e-01 1.579e-01 2.848e-01 3.038e-01 +44 7.779e-01 7.536e-01 2.425e-02 6.087e-01 4.443e-01 1.644e-01 2.964e-01 3.158e-01 +45 7.868e-01 7.616e-01 2.521e-02 6.199e-01 4.490e-01 1.709e-01 3.081e-01 3.278e-01 +46 7.954e-01 7.693e-01 2.616e-02 6.308e-01 4.535e-01 1.773e-01 3.197e-01 3.397e-01 +47 8.037e-01 7.766e-01 2.710e-02 6.415e-01 4.578e-01 1.837e-01 3.313e-01 3.516e-01 +48 8.117e-01 7.836e-01 2.804e-02 6.521e-01 4.619e-01 1.901e-01 3.428e-01 3.633e-01 +49 8.193e-01 7.903e-01 2.898e-02 6.624e-01 4.659e-01 1.965e-01 3.542e-01 3.750e-01 +50 8.267e-01 7.968e-01 2.991e-02 6.724e-01 4.697e-01 2.028e-01 3.656e-01 3.866e-01 +51 8.337e-01 8.029e-01 3.083e-02 6.823e-01 4.733e-01 2.090e-01 3.769e-01 3.981e-01 +52 8.405e-01 8.088e-01 3.175e-02 6.920e-01 4.768e-01 2.152e-01 3.881e-01 4.095e-01 +53 8.471e-01 8.144e-01 3.266e-02 7.015e-01 4.801e-01 2.214e-01 3.992e-01 4.208e-01 +54 8.534e-01 8.198e-01 3.356e-02 7.108e-01 4.833e-01 2.275e-01 4.102e-01 4.319e-01 +55 8.594e-01 8.250e-01 3.445e-02 7.199e-01 4.863e-01 2.335e-01 4.211e-01 4.430e-01 +56 8.653e-01 8.299e-01 3.533e-02 7.288e-01 4.893e-01 2.395e-01 4.319e-01 4.539e-01 +57 8.709e-01 8.347e-01 3.621e-02 7.375e-01 4.920e-01 2.454e-01 4.426e-01 4.647e-01 +58 8.763e-01 8.392e-01 3.707e-02 7.460e-01 4.947e-01 2.513e-01 4.531e-01 4.754e-01 +59 8.815e-01 8.436e-01 3.792e-02 7.544e-01 4.973e-01 2.571e-01 4.636e-01 4.859e-01 +60 8.865e-01 8.477e-01 3.877e-02 7.625e-01 4.997e-01 2.628e-01 4.739e-01 4.963e-01 +61 8.913e-01 8.517e-01 3.960e-02 7.705e-01 5.021e-01 2.685e-01 4.841e-01 5.066e-01 +62 8.959e-01 8.555e-01 4.043e-02 7.784e-01 5.043e-01 2.741e-01 4.942e-01 5.167e-01 +63 9.004e-01 8.591e-01 4.124e-02 7.860e-01 5.065e-01 2.796e-01 5.041e-01 5.268e-01 +64 9.047e-01 8.626e-01 4.204e-02 7.935e-01 5.085e-01 2.850e-01 5.139e-01 5.367e-01 +65 9.088e-01 8.660e-01 4.284e-02 8.009e-01 5.105e-01 2.904e-01 5.236e-01 5.464e-01 +66 9.128e-01 8.692e-01 4.362e-02 8.081e-01 5.124e-01 2.957e-01 5.332e-01 5.560e-01 +67 9.166e-01 8.722e-01 4.439e-02 8.151e-01 5.142e-01 3.009e-01 5.426e-01 5.654e-01 +68 9.203e-01 8.751e-01 4.515e-02 8.220e-01 5.159e-01 3.061e-01 5.519e-01 5.747e-01 +69 9.238e-01 8.779e-01 4.590e-02 8.287e-01 5.175e-01 3.111e-01 5.610e-01 5.838e-01 +70 9.273e-01 8.806e-01 4.663e-02 8.353e-01 5.191e-01 3.161e-01 5.700e-01 5.928e-01 +71 9.306e-01 8.832e-01 4.736e-02 8.417e-01 5.206e-01 3.211e-01 5.789e-01 6.016e-01 +72 9.337e-01 8.857e-01 4.807e-02 8.480e-01 5.221e-01 3.259e-01 5.876e-01 6.103e-01 +73 9.368e-01 8.880e-01 4.878e-02 8.542e-01 5.235e-01 3.307e-01 5.962e-01 6.188e-01 +74 9.397e-01 8.903e-01 4.947e-02 8.602e-01 5.248e-01 3.354e-01 6.047e-01 6.271e-01 +75 9.426e-01 8.924e-01 5.015e-02 8.661e-01 5.261e-01 3.400e-01 6.130e-01 6.353e-01 +76 9.453e-01 8.945e-01 5.082e-02 8.718e-01 5.273e-01 3.445e-01 6.212e-01 6.433e-01 +77 9.479e-01 8.965e-01 5.148e-02 8.774e-01 5.285e-01 3.490e-01 6.293e-01 6.512e-01 +78 9.505e-01 8.984e-01 5.213e-02 8.830e-01 5.296e-01 3.534e-01 6.372e-01 6.589e-01 +79 9.529e-01 9.002e-01 5.276e-02 8.883e-01 5.307e-01 3.577e-01 6.449e-01 6.665e-01 +80 9.553e-01 9.019e-01 5.339e-02 8.936e-01 5.317e-01 3.619e-01 6.526e-01 6.740e-01 +81 9.576e-01 9.036e-01 5.400e-02 8.987e-01 5.327e-01 3.661e-01 6.601e-01 6.812e-01 +82 9.598e-01 9.052e-01 5.460e-02 9.038e-01 5.336e-01 3.702e-01 6.675e-01 6.884e-01 +83 9.619e-01 9.067e-01 5.520e-02 9.087e-01 5.345e-01 3.742e-01 6.747e-01 6.954e-01 +84 9.639e-01 9.081e-01 5.578e-02 9.135e-01 5.353e-01 3.781e-01 6.818e-01 7.022e-01 +85 9.659e-01 9.095e-01 5.635e-02 9.182e-01 5.362e-01 3.820e-01 6.888e-01 7.089e-01 +86 9.678e-01 9.109e-01 5.691e-02 9.228e-01 5.370e-01 3.858e-01 6.957e-01 7.154e-01 +87 9.696e-01 9.121e-01 5.746e-02 9.273e-01 5.377e-01 3.895e-01 7.024e-01 7.219e-01 +88 9.714e-01 9.134e-01 5.800e-02 9.316e-01 5.384e-01 3.932e-01 7.090e-01 7.281e-01 +89 9.731e-01 9.145e-01 5.853e-02 9.359e-01 5.391e-01 3.968e-01 7.155e-01 7.343e-01 +90 9.747e-01 9.157e-01 5.905e-02 9.401e-01 5.398e-01 4.003e-01 7.218e-01 7.403e-01 +91 9.763e-01 9.167e-01 5.956e-02 9.442e-01 5.404e-01 4.038e-01 7.281e-01 7.461e-01 +92 9.778e-01 9.178e-01 6.006e-02 9.482e-01 5.410e-01 4.072e-01 7.342e-01 7.519e-01 +93 9.793e-01 9.187e-01 6.055e-02 9.521e-01 5.416e-01 4.105e-01 7.402e-01 7.575e-01 +94 9.807e-01 9.197e-01 6.103e-02 9.559e-01 5.422e-01 4.137e-01 7.460e-01 7.629e-01 +95 9.821e-01 9.206e-01 6.150e-02 9.596e-01 5.427e-01 4.169e-01 7.518e-01 7.683e-01 +96 9.834e-01 9.215e-01 6.196e-02 9.633e-01 5.432e-01 4.201e-01 7.574e-01 7.735e-01 +97 9.847e-01 9.223e-01 6.242e-02 9.668e-01 5.437e-01 4.231e-01 7.630e-01 7.786e-01 +98 9.859e-01 9.231e-01 6.286e-02 9.703e-01 5.441e-01 4.261e-01 7.684e-01 7.836e-01 +99 9.871e-01 9.238e-01 6.330e-02 9.737e-01 5.446e-01 4.291e-01 7.737e-01 7.884e-01 +100 9.883e-01 9.246e-01 6.372e-02 9.770e-01 5.450e-01 4.320e-01 7.789e-01 7.932e-01 +101 9.894e-01 9.252e-01 6.414e-02 9.802e-01 5.454e-01 4.348e-01 7.840e-01 7.978e-01 +102 9.905e-01 9.259e-01 6.455e-02 9.834e-01 5.458e-01 4.376e-01 7.890e-01 8.023e-01 +103 9.915e-01 9.265e-01 6.495e-02 9.865e-01 5.462e-01 4.403e-01 7.939e-01 8.067e-01 +104 9.925e-01 9.272e-01 6.534e-02 9.895e-01 5.466e-01 4.429e-01 7.987e-01 8.110e-01 +105 9.935e-01 9.277e-01 6.572e-02 9.924e-01 5.469e-01 4.455e-01 8.034e-01 8.151e-01 +106 9.944e-01 9.283e-01 6.610e-02 9.953e-01 5.472e-01 4.481e-01 8.079e-01 8.192e-01 +107 9.953e-01 9.288e-01 6.646e-02 9.981e-01 5.475e-01 4.506e-01 8.124e-01 8.232e-01 +108 9.962e-01 9.293e-01 6.682e-02 1.001e+00 5.478e-01 4.530e-01 8.168e-01 8.270e-01 +109 9.970e-01 9.298e-01 6.718e-02 1.004e+00 5.481e-01 4.554e-01 8.211e-01 8.308e-01 +110 9.978e-01 9.303e-01 6.752e-02 1.006e+00 5.484e-01 4.577e-01 8.253e-01 8.344e-01 +111 9.986e-01 9.307e-01 6.786e-02 1.009e+00 5.487e-01 4.600e-01 8.295e-01 8.380e-01 +112 9.994e-01 9.312e-01 6.819e-02 1.011e+00 5.489e-01 4.623e-01 8.335e-01 8.415e-01 +113 1.000e+00 9.316e-01 6.851e-02 1.014e+00 5.492e-01 4.644e-01 8.374e-01 8.448e-01 +114 1.001e+00 9.320e-01 6.883e-02 1.016e+00 5.494e-01 4.666e-01 8.413e-01 8.481e-01 +115 1.001e+00 9.324e-01 6.913e-02 1.018e+00 5.496e-01 4.687e-01 8.451e-01 8.513e-01 +116 1.002e+00 9.327e-01 6.944e-02 1.021e+00 5.498e-01 4.707e-01 8.488e-01 8.544e-01 +117 1.003e+00 9.331e-01 6.973e-02 1.023e+00 5.500e-01 4.727e-01 8.524e-01 8.574e-01 +118 1.003e+00 9.334e-01 7.002e-02 1.025e+00 5.502e-01 4.747e-01 8.559e-01 8.603e-01 +119 1.004e+00 9.337e-01 7.030e-02 1.027e+00 5.504e-01 4.766e-01 8.594e-01 8.632e-01 +120 1.005e+00 9.340e-01 7.058e-02 1.029e+00 5.506e-01 4.785e-01 8.627e-01 8.659e-01 +121 1.005e+00 9.343e-01 7.085e-02 1.031e+00 5.508e-01 4.803e-01 8.660e-01 8.686e-01 +122 1.006e+00 9.346e-01 7.111e-02 1.033e+00 5.509e-01 4.821e-01 8.693e-01 8.712e-01 +123 1.006e+00 9.348e-01 7.137e-02 1.035e+00 5.511e-01 4.838e-01 8.724e-01 8.737e-01 +124 1.007e+00 9.351e-01 7.162e-02 1.037e+00 5.512e-01 4.855e-01 8.755e-01 8.761e-01 +125 1.007e+00 9.353e-01 7.187e-02 1.039e+00 5.514e-01 4.872e-01 8.785e-01 8.785e-01