Logistic curve fitting in epidemiology: Difference between revisions

→‎{{header|Wren}}: Now uses new core library method.
(→‎{{header|Wren}}: Now uses new core library method.)
Line 1,158:
=={{header|Wren}}==
{{trans|Phix}}
<lang ecmascript>var K = 7800000000 // approx world population
{{libheader|Wren-math}}
<lang ecmascript>import "/math" for Math
 
var K = 7800000000 // approx world population
var n0 = 27 // number of cases at day 0
 
Line 1,182 ⟶ 1,179:
var sq = 0
for (i in 0...y.count) {
var eri = Math.exp(r*i).exp
var dst = (n0*eri)/(1+n0*(eri-1)/K) - y[i]
sq = sq + dst * dst
Line 1,213 ⟶ 1,210:
 
var r = (solve.call(f, 0.5, 0) * 1e10).round / 1e10
var R0 = (Math.exp(12 * r).exp * 1e8).round / 1e8
System.print("r = %(r), R0 = %(R0)")</lang>
 
9,486

edits