In this paper, probabilistic distribution of chassis component fatigue life is determined statistically by applying the design of experiments and the Pearson system. To construct p - ${varepsilon}$ - N curve, the case that fatigue data are random variables is attempted. Probabilistic density function (p.d.f) for fatigue life is obtained by the design of experiment and using this p.d.f fatigue reliability, any aimed fatigue life can be calculated. Lower control arm and rear torsion bar of chassis components are selected as examples for analysis. Component load histories which are obtained by multi-body dynamic simulation for Belsian load history are used. Finite element analysis is performed by using commercial software MSC Nastran and fatigue analysis is performed by using FE Fatigue. When strain-life curve itself is random variable, the probability density function of fatigue life has very little difference from log-normal distribution. And the cases of fatigue data are random variables, probability density functions are approximated to Beta distribution. Each p.d.f is verified by Monte-Carlo simulation.