2b:null 2c:["$","title",null,{"children":"Interactive SEIR Model | S. Scott Collis | Observable"}] 31:T497,viewof logPlot = { const min = vl.markLine({clip:true}).encode( vl.x().fieldQ("Time").scale({domain:[2,xscale]}), vl.y().fieldQ("Min").axis({tickCount:5}),//.scale({type:"log",base:10,domain:[1,P]}), vl.color().fieldN("Compartment").sort(vl.max("Max").order("descending")) ); const max = vl.markLine({clip:true}).encode( vl.x().fieldQ("Time").scale({domain:[2,xscale]}), vl.y().fieldQ("Max").axis({tickCount:5}),//.scale({type:"log",base:10,domain:[1,P]}), vl.color().fieldN("Compartment").sort(vl.max("Max").order("descending")) ); const area = vl.markArea({clip:true}).encode( vl.x().fieldQ("Time").scale({domain:[2,xscale]}), vl.y().fieldQ("Max").axis({tickCount:5}),//.scale({type:"log",base:10,domain:[1,P]}), vl.y2().fieldQ("Min"), vl.color().fieldN("Compartment").sort(vl.max("Max").order("descending")) ); return vl.data(Urange) .transform({filter:{and:[{field:"Compartment",oneOf:["Rh","Dh"]}, //{field:"Min",gt:0.01},{field:"Max",gt:0.01} ]}}) .layer([min,max,area]) //.width(800) .title("All samples") .config({axisX:{ticks:true}}) .render(); }32:T100e,class SEIR { constructor(P_=P,Io_=Io,rho_=rho,alpha_=5.2,gamma_=3.0,c_=c, to_=0,tf_=tf,nt_=750,Ro_=Ro,Fa_=Fa, stochastic_=stochastic,samples_=samples, k_=k,n_=n) { this.P = P_; this.Io = Io_; this.Ro = Ro_; // FIXME: should be an array this.tk = [0, 30, 80]; this.alpha = alpha_; this.gamma = gamma_; this.rho = rho_; // FIXME: should be an array this.c = c_; this.to = to_; this.tf = tf_; this.nt = nt_; this.dt = (this.tf-this.to)/this.nt; this.Fa = Fa_; this.phase = 0; this.erlang = new Erlang(k,n); this.stochastic = stochastic; this.alphai; // current sample this.gammai; // current sample this.samples = samples_; } get beta() { return this.Ro/this.gammai; // FIXME: this.Ro[this.phase]; } get setIC() { var IC = [this.P-this.Io, 0, 0, this.Io, 0, 0, 0, 0]; return IC; } rk4(yo,yf,to,h) { var f = yo.slice(); var k1 = yo.slice(); var k2 = yo.slice(); var k3 = yo.slice(); var k4 = yo.slice(); var q = yo.slice(); let neq = yo.length; this.rhs(yo,to,f); if (1) { for (let j=0; j