Add new rib fracture risk function to standard functions
Larsson et al. deveolped a new rib fracture risk function (https://doi.org/10.3389/fbioe.2021.677768) which is shown below. It would be great if someone could integrate this into the standard functions.
@staticmethod
def erf_rib_risk_age(x: pint.Quantity, age: float):
# parameters from Larsson et al. 2021:
# https://doi.org/10.3389/fbioe.2021.677768
a = 0.3026
b = -2.9866
c = -0.013
return (0.5+0.5*special.erf((np.log((x))-(b+c*age))/(math.sqrt(2)*a)))
Cheers, Felix