Pile load tests have been utilized to reduce the uncertainty of pile resistance, thus leading to a higher resistance factor used in the Load and Resistance Factor Design (LRFD). Previous studies have primarily focused on calibrating resistance factors for single piles based on load tests. This calibration hinges upon the resistance bias factor of single piles, defined as the ratio of measured resistance to predicted resistance. Due to the redundancy in the pile group system, it is conventionally assumed that if the individual piles within the group achieve a lower reliability index (e.g., 2.0–2.5), the pile group as a whole attains the target reliability index of 3. However, the approach is empirical as it does not consider system redundancy directly. Moreover, this empirical approach disregards the correlation between resistance bias factors of individual piles, which is inherently influenced by the spatial variability of soils. In this study, the random finite difference method (RFDM) is employed to evaluate the correlation between resistance bias factors of individual piles in spatially variable soils. The resultant correlation matrix is subsequentially employed in Bayes’ theorem to update resistance bias factors using individual pile load test results and their corresponding test locations. The updated resistance bias factors are then used for the direct calibration of resistance factors for pile groups within the framework of LRFD. A pile group subject to vertical loading in undrained clays is adopted for illustration. Comparative analyses between the proposed approach and the empirical approach demonstrate that the latter tends to overestimate the resistance factor. Furthermore, the proposed approach enables the determination of optimal locations for conducting subsequent load tests based on previous test results.