Data which gathered and used in the field of fatigue and fracture mechanics have a lot of uncertainties. In this case, those uncertainties will make scatter band in evaluation of fatigue life and fracture toughness. Thus, the probabilistic analysis of these data will be needed. For determining the fatigue life in mixed mode, using crack direction law and fatigue crack growth law, the problem is studied as a constrained life minimization. Stress intensity factor(SIF) is computed by approximate solution table(Ewalds/Wanhill 1984) and 0th order PFEM. The variance of fatigue life and SIF are computed by differentiation of tabulated approximate solution and 1st order PFEM. And these are used for criterion of design values, principal parameter determination and modelling. The problem of center cracked plate is solved for checking the PFEM model which is influenced by various parameters like as initial crack length, final crack length, two fatigue parameters in Paris Equation and applied stress.