Models in lonewolfbio

@hide
// units
k = 10^3
m = 10^6
b = 10^9

// Level Of Experience
n_skill_foreign_discount = beta(49.7510591728389, 2.3576183042092) // 90 to 99%

n_skill_nonexpert = (350 * k to 500 * k) * n_skill_foreign_discount

k = 10^3
m = 10^6
b = 10^9
t= 10^12

n_phd = 20*k to 200*k
n_stem = 5*m to 50*m
n_ee = 50*m to 500*m

b_phd = truncate(0.03 to 0.30,0.01,1)

// h/t to Ozzie!

@name("Human Reliability Assessment Calculator")
@doc("Assumes no correlations between different tasks and attempts")
@startClosed

max_v=1
adjusted_error(error_rate) = mx(truncateRight(error_rate,1),max_v,[cdf(error_rate, 1),1-cdf(error_rate, 1)])
success_at_protocol(error_rate, numberSteps) = exp(-numberSteps * error_rate)
success_at_protocol_binomial(error_rate,numberSteps) = (1-(mx(truncateRight(error_rate,1),max_v,[cdf(error_rate, 1),1-cdf(error_rate, 1)])))^numberSteps

t = 10^12
b = 10^9
m = 10 ^ 6
k = 1000

population = 210 * m

fraction_believers = 0.05% to 0.5%
believers = population * fraction_believers

// Units

t = 10^12
b = 10^9
m = 10 ^ 6
k = 1000

// Skill by time

months = List.upTo(1, 24)