Clarity on calculation of parameters apha and beta

[edger-asiimwe]

The empting rate, kempt(qsto) of the liquid compartment of the stomach in this code is described below which I understand is kqut in this case.

if Dbar > 0:
    aa = 5 / 2 / (1 - params.b) / Dbar
    cc = 5 / 2 / params.d / Dbar
    kgut = params.kmin + (params.kmax - params.kmin) / 2 * (
        np.tanh(aa * (qsto - params.b * Dbar)) -
        np.tanh(cc * (qsto - params.d * Dbar)) + 2)
else:
    kgut = params.kmax

In this paper, A System Model of the Oral Glucose Absorption: Validation on Gold Standard Data shows that the parameters alpha (aa) and beta (cc) are as follows: Translating it to python would yield the code below

aa = 5 / ( 2 * Dbar * (1 - params.b))
cc = 5 / (2 * Dbar * params.d)

I understand both yield the same result. For purposes of consistency and readability of the code with the equations in the @paper, do you think we can adopt the later option for this case?

@jxx123

[jxx123]

Sure, do you mind sending me a Pull Request? I can stamp it.

[edger-asiimwe]

I have updated the parameter alpha and beta via this pull request

Thanks!