def demand_mass_balance_h(odemand, close_prob, recapture_rate):
"""Solve Demand Mass Balance equation for host-level
Parameters
----------
odemand: int
Observerd demand
close_prob: float
Probability of selecting a closed element from host set H
recapture_rate: float
Estimated recapture rate
Returns
-------
tuple
Estimated demand, spill and recapture
"""
A = np.array([[1, -1, 1], [-close_prob, 1, 0], [0, -recapture_rate, 1]])
B = np.array([odemand, 0, 0])
demand, spill, recapture = solve(A, B)
return demand, spill, recapture
评论列表
文章目录