In this paper, the procedure to estimate fatigue crack initiation life has been established by considering fretting wear and multiaxial stress states on the contact surface of press-fitted shafts. And a method to calculate the local friction coefficient during the running-in period of fretting wear process has been proposed. The predicted result of worn surface profile in the press-fitted shaft with non-linear local friction coefficient can avoid excessive wear depth estimation compared with that for the case of constant local friction coefficient. Furthermore, the predicted fatigue crack initiation lives based on Smith-Watson-Topper model considering the fretting wear are in good agreement with the experimental data. Consequently, the present method is valid not only for predicting worn surface profile, but also for assessing fatigue crack initiation lives considering the fretting wear during the running-in period in press fits.