A global study of geometry and material properties of orthopedic screws was performed, considering not only the effect of each single factor (screw pitch, number of threads, fillet angle, etc.) but also their interactions with respect to bone density.The stress patterns resulting from different screw geometries and bone densities were analyzed using finite element techniques, taking into account different levels of osseointegration between the screw and the bone. These numerical models where validated through experimental pull-out tests, where a pull out force of 120 N produced localized failure of the last thread (stresses above 0.42 MPa). The results of the numerical simulations were then summarised using a multi-factorial parametric analysis. This demonstrated the great relevance of the interaction between bone density and screw pitch, showing that the optimal screw pitch can vary by more than 25% for different densities (0.35 g/cm(3) and 0.47 g/cm(3), respectively).The parameters calculated by means of the multi-factorial analysis allow the pull out force to be estimated for different osseointegration levels, different screw geometries and material properties, and for different bone densities. The final objective is to determine the best choice of implant for each individual patient.