Published
Edited
Apr 24, 2020
N1_step = (N1) => (1-d)*N1+(d*K)*(f1*N1)/(f1*N1+f2()*(K-N1))
f2 = () => {
var f2 = f1*f2divf1
var delta = 1/2*f2*deltaf2
return f2 + (Math.random() < 0.5 ? -delta : delta)
}
timecourse = {
run_button;
var timecourse = Array(50000);
timecourse[0] = {t:0, N1:N1_0};
for (var i = 1; i < timecourse.length; i++) {
timecourse[i] = {t:i, N1:N1_step(timecourse[i-1].N1)};
};
return timecourse
}
K = 1
f1 = 1.0
delta_N1_per_N1 = {
var data = Array(16);
for (var i = 0; i < 16; i++) {
var N1 = i/15
var f2_high = f1*f2divf1 + 1/2*f1*f2divf1*deltaf2
var f2_low = f1*f2divf1 - 1/2*f1*f2divf1*deltaf2
var N1plus_low = (1-d)*N1+(d*K)*(f1*N1)/(f1*N1+f2_low*(K-N1))
var N1plus_high = (1-d)*N1+(d*K)*(f1*N1)/(f1*N1+f2_high*(K-N1))
var av_N1plus = (N1plus_low+N1plus_high)/2
data[i] = {N1:N1, deltaN1:av_N1plus - N1};
};
return data
}
e2_N1_step = (N1) => {
var f1 = e2_f1()
return (1-e2_d)*N1+(e2_d*K)*(f1*N1)/(f1*N1+e2_f2()*(K-N1))}
e2_f2 = () => {
var f2 = f1*e2_f2divf1
var delta = 1/2*e2_delta
return f2 + (Math.random() < 0.5 ? -delta : delta)
}
e2_f1 = () => {
var delta = 1/2*e2_delta
return f1 + (Math.random() < 0.5 ? -delta : delta)
}
e2_timecourse = {
e2_run_button;
//var epsilon = 10e-9
var timecourse = Array(50000);
timecourse[0] = {t:0, N1:e2_N1_0};
for (var i = 1; i < timecourse.length; i++) {
var next_N2 = e2_N1_step(timecourse[i-1].N1)
/*if (next_N2 < epsilon || timecourse[i-1].N1 < 0){
next_N2 = 0}
else if (next_N2 > (1-epsilon) || timecourse[i-1].N1 > (1-epsilon)) {
next_N2 = 1}*/
timecourse[i] = {t:i, N1:next_N2};
};
return timecourse
}
e2_delta_N1_per_N1 = {
var data = Array(16);
for (var i = 0; i < 16; i++) {
var N1 = i/15
var f2_high = f1*e2_f2divf1 + 1/2*e2_f2divf1*e2_delta
var f2_low = f1*e2_f2divf1 - 1/2*e2_f2divf1*e2_delta
var f1_high = f1 + 1/2*e2_delta
var f1_low = f1 - 1/2*e2_delta
var N1plus = Array(4)
N1plus[0] = (1-d)*N1+(d*K)*(f1_low*N1)/(f1_low*N1+f2_low*(K-N1))
N1plus[1] = (1-d)*N1+(d*K)*(f1_low*N1)/(f1_low*N1+f2_high*(K-N1))
N1plus[3] = (1-d)*N1+(d*K)*(f1_high*N1)/(f1_high*N1+f2_low*(K-N1))
N1plus[4] = (1-d)*N1+(d*K)*(f1_high*N1)/(f1_high*N1+f2_high*(K-N1))
var av_N1plus = _.sum(N1plus)/4
data[i] = {N1:N1, deltaN1:av_N1plus - N1};
};
return data
}
accent_color = "#2980b9"//"dodgerblue"
PD = require('https://bundle.run/probability-distributions@0.9.1')
One platform to build and deploy the best data apps
Experiment and prototype by building visualizations in live JavaScript notebooks. Collaborate with your team and decide which concepts to build out.
Use Observable Framework to build data apps locally. Use data loaders to build in any language or library, including Python, SQL, and R.
Seamlessly deploy to Observable. Test before you ship, use automatic deploy-on-commit, and ensure your projects are always up-to-date.