High cycle fatigue tests at a constant positive mean stress have been performed on a Al-Si-Cu cast aluminum alloy. The Random Fatigue Limit (RFL) model was employed to fit the probabilistic S-N curves based on Maximum Likelihood Estimate (MLE). Fractographic studies indicated that fatigue cracks in most specimens initiate from oxide films located at or very close to specimen surface. The RFL model was proved to be able to accurately capture the scatter in fatigue life. The cumulative density function (CDF) of fatigue life determined by RFL fit is found to be approximately equal to the complementary value of the CDF of the near-surface fatigue initiator size.