A parametric model based computer implemented method for determining the stiffness of a bone, comprising: obtaining measurement data of oscillatory forces as functions of time F(t) and obtaining resulting oscillatory acceleration data (a) as functions of time a(t) with a data receiver; wherein said measurement data are derived from applying a superposition of static and oscillatory forces over a range of frequencies (f) to a region of a skin-bone complex thereby exciting oscillatory accelerations (a) over the range of frequencies (f) of the skin-bone complex; wherein said data receiver is communicatively coupled to a controller comprising a processor and a storage medium containing computer readable and executable instructions which, when executed by the processor, cause the controller to automatically: (i) transform a(t) and F(t) to functions of frequency, a(f) and F(f), (ii) reduce a(f) and F(f) to accelerance frequency response function data A(f), (iii) determine, a complex compliance frequency response function, Y(f) and associated complex stiffness frequency response function H(f), (iv) fit a parametric mathematical model to Y(f) to obtain a first set of parameters of the parametric mathematical model, including the stiffness of the bone (Kb), (v) fit the parametric mathematical model to H(f) to obtain a second set of parameters of the parametric mathematical model, including the stiffness of the bone (Kb), (vi) after fitting, determine discrepancies between each of the parameters of the first set of parameters and the second set of parameters as a measure of conformity thereof to the parametric mathematical model, and (vii) save the measure of conformity, the first set of parameters, and the second set of parameters as a parameter set for said skin region; repeating step, wherein the measurements in step are derived from applying the static and oscillatory forces to a lateral shifted region of the skin-bone complex, thereby obtaining a parameter set for the lateral shifted region; repeating step until the conformity worsens, and determining an optimized parameter set based on the measure of conformity; and determining the stiffness of the bone from (kb) values of the optimized parameter set. obtaining measurement data of oscillatory forces as functions of time F(t) and obtaining resulting oscillatory acceleration data (a) as functions of time a(t) with a data receiver; wherein said measurement data are derived from applying a superposition of static and oscillatory forces over a range of frequencies (f) to a region of a skin-bone complex thereby exciting oscillatory accelerations (a) over the range of frequencies (f) of the skin-bone complex; wherein said data receiver is communicatively coupled to a controller comprising a processor and a storage medium containing computer readable and executable instructions which, when executed by the processor, cause the controller to automatically: (i) transform a(t) and F(t) to functions of frequency, a(f) and F(f), (ii) reduce a(f) and F(f) to accelerance frequency response function data A(f), (iii) determine, a complex compliance frequency response function, Y(f) and associated complex stiffness frequency response function H(f), (iv) fit a parametric mathematical model to Y(f) to obtain a first set of parameters of the parametric mathematical model, including the stiffness of the bone (Kb), (v) fit the parametric mathematical model to H(f) to obtain a second set of parameters of the parametric mathematical model, including the stiffness of the bone (Kb), (vi) after fitting, determine discrepancies between each of the parameters of the first set of parameters and the second set of parameters as a measure of conformity thereof to the parametric mathematical model, and (vii) save the measure of conformity, the first set of parameters, and the second set of parameters as a parameter set for said skin region; (i) transform a(t) and F(t) to functions of frequency, a(f) and F(f), (ii) reduce a(f) and F(f) to accelerance frequency response function data A(f), (iii) determine, a complex compliance frequency response function, Y(f) and associated complex stiffness frequency response function H(f), (iv) fit a parametric mathematical model to Y(f) to obtain a first set of parameters of the parametric mathematical model, including the stiffness of the bone (Kb), (v) fit the parametric mathematical model to H(f) to obtain a second set of parameters of the parametric mathematical model, including the stiffness of the bone (Kb), (vi) after fitting, determine discrepancies between each of the parameters of the first set of parameters and the second set of parameters as a measure of conformity thereof to the parametric mathematical model, and (vii) save the measure of conformity, the first set of parameters, and the second set of parameters as a parameter set for said skin region; repeating step, wherein the measurements in step are derived from applying the static and oscillatory forces to a lateral shifted region of the skin-bone complex, thereby obtaining a parameter set for the lateral shifted region; repeating step until the conformity worsens, and determining an optimized parameter set based on the measure of conformity; and determining the stiffness of the bone from (kb) values of the optimized parameter set.