Parametric model based computer implemented methods for determining the stiffness of a bone and systems for estimating the stiffness of a bone in vivo. The computer implemented methods include determining a complex compliance frequency response function Y(f) and an associated complex stiffness frequency response function H(f) and independently fitting a parametric mathematical model to Y(f) and to H(f). The systems include a device for measuring the stiffness of the bone in vivo and a data analyzer to determine a complex compliance frequency response function Y(f) and an associated complex stiffness frequency response function H(f).