全球范围内的脉冲星测时阵列（PTA）实验在PTA数据集中揭示了支持现有引力波背景（GWB）信号的证据。除了获取更多观测数据外，通过改进噪声建模的准确性也可以提高PTA实验的灵敏度。在PTA数据分析中，噪声建模主要采用贝叶斯统计方法，依赖边缘似然和贝叶斯因子来评估证据。我们引入广义步进采样（GSS）作为PTA-贝叶斯框架中一种高效且准确的边缘似然估计方法。该方法能够以低成本获得高精度估计，特别是在比较昂贵的模型（如Hellings-Downs（HD）模型或重叠归一化函数模型（ORF））时。我们通过重新评估北美纳米赫兹引力波天文台（NANOGrav）15年数据集和欧洲PTA（EPTA）第二数据发布的结果，展示了GSS在模型选择和证据计算中的效率和准确性。我们发现与NANOGrav 15年数据集报告的结果相比，GWB的证据相似；与EPTA第二数据发布的证据相比，我们在所有数据集中发现了显著增加的GWB支持证据。

Globally, Pulsar Timing Array (PTA) experiments have revealed evidence supporting an existing gravitational wave background (GWB) signal in the PTA data set. Apart from acquiring more observations, the sensitivity of PTA experiments can be increased by improving the accuracy of the noise modeling. In PTA data analysis, noise modeling is conducted primarily using Bayesian statistics, relying on the marginal likelihood and the Bayes factor to assess the evidence. We introduce generalized steppingstone (GSS) as an efficient and accurate marginal likelihood estimation method for the PTA-Bayesian framework. This method enables low-cost estimates with high accuracy, especially when comparing expensive models such as the Hellings-Downs (HD) model or the overlap reduction function model (ORF). We demonstrate the efficiency and the accuracy of GSS for model selection and evidence calculation by reevaluating the evidence of previous analyses from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 15 yr data set and the European PTA (EPTA) second data release. We find similar evidence for the GWB compared to the one reported by the NANOGrav 15-year data set. Compared to the evidence reported for the EPTA second data release, we find a substantial increase in evidence supporting GWB across all data sets.