Machine Learning and Frequency–Severity Decomposition for Insurance Pricing

1. Introduction

Accurate estimation of expected losses is a central challenge in automobile insurance pricing. Insurers must determine premiums that reflect the underlying risk of each policyholder while ensuring fairness, competitiveness, and regulatory compliance, as premium estimation directly affects portfolio stability and profitability.

A cornerstone of actuarial practice is the frequency–severity framework, which decomposes the expected loss—or pure premium—into two latent components: the expected number of claims and the expected claim size. By modeling these components separately, actuaries can employ statistical distributions tailored to the distinct stochastic properties of each process before recombining them for final loss estimation. This framework provides both interpretability and flexibility, making it fundamental to modern insurance ratemaking.

Generalized Linear Models (GLMs) have long served as the industry standard for actuarial modeling. Existing literature highlights the capacity of GLMs to provide an interpretable yet flexible structure for non-normal, skewed, and heteroscedastic insurance data [3,7,8]. While Poisson and Negative Binomial models are commonly used for claim frequency, Gamma models are widely applied to claim severity due to the strictly positive and right-skewed nature of loss amounts. Extensions such as Generalized Additive Models (GAMs) allow for nonlinear relationships while preserving model transparency [4,10].

Machine learning (ML) has emerged as a powerful tool in insurance analytics. Tree-based algorithms and gradient boosting methods are particularly effective in capturing complex interactions and nonlinear relationships that traditional parametric models may overlook. Empirical evidence suggests that models such as XGBoost can significantly improve performance in tasks ranging from claim frequency prediction to fraud detection [2,6,9].

Recent studies have explored both direct and component-based modeling strategies using machine learning. Direct modeling of total claim amounts, using methods such as Support Vector Regression (SVR), XGBoost, and neural networks, has demonstrated strong predictive accuracy for aggregate loss estimation [18]. Other studies incorporate machine learning within the classical frequency–severity framework, where gradient boosting models are compared with traditional GLMs [19]. These findings indicate that machine learning consistently improves claim frequency prediction, while results for claim severity remain less conclusive due to the high variability of loss amounts.

Despite these advances, direct modeling approaches do not explicitly account for the structural decomposition of insurance losses into frequency and severity components. Moreover, existing evidence on the relative performance of decomposition-based and direct modeling approaches remains inconclusive, particularly when evaluated under both predictive accuracy and risk segmentation criteria.

Classical actuarial models continue to play an important role due to their interpretability and regulatory transparency. Consequently, recent research has focused on hybrid approaches that integrate machine learning techniques within the frequency–severity framework [11,12]. These approaches aim to enhance predictive performance while preserving the structural foundations of actuarial pricing.

Alternative unified modeling approaches have also been proposed, most notably through the Tweedie distribution, which provides a compound Poisson–Gamma representation of aggregate losses [15,16]. While such models offer strong theoretical appeal, it remains unclear whether they can outperform decomposition-based methods in practical insurance applications.

Motivated by these considerations and the lack of consensus in the literature, this study investigates whether machine learning improves insurance pricing by replacing the classical decomposition framework or by enhancing its individual components. In particular, we evaluate model performance not only in terms of prediction accuracy but also in terms of the ability to identify high-risk policies.

To address these questions, we develop a comprehensive modeling framework that integrates classical actuarial models with modern machine learning techniques. We evaluate several models for claim frequency, including Poisson GLMs, spline-based models, and XGBoost. For claim severity, we compare a Gamma GLM with XGBoost. These models are combined within the frequency–severity decomposition framework to estimate pure premium. In addition, we examine direct modeling approaches, including XGBoost and Tweedie models, for predicting aggregate losses.

The contribution of this study is threefold. First, we demonstrate that machine learning methods significantly improve predictive performance for both claim frequency and claim severity, particularly in capturing nonlinear relationships and interactions. Second, we provide empirical evidence that the frequency–severity decomposition framework remains superior to direct modeling approaches, even when advanced machine learning methods are applied. Third, we identify a trade-off between prediction accuracy and risk segmentation: while models combining machine learning for both components achieve the best overall accuracy, models emphasizing frequency modeling provide superior identification of high-risk policies.

Overall, this study contributes to the literature on actuarial science and financial mathematics by providing a comprehensive comparison of classical and modern approaches to insurance pricing. The results highlight the continued relevance of the frequency–severity decomposition framework while demonstrating how machine learning can be effectively integrated to enhance predictive performance and risk differentiation.

The remainder of the paper is organized as follows. Section 2 describes the dataset and data preparation procedures; Section 3 presents the modeling framework and statistical methods; Section 4 reports the empirical findings; and Section 5 concludes with a discussion of the results and their implications for insurance pricing and risk management.

2. Data Description and Preparation

This study utilizes the French motor third-party liability (freMTPL2) datasets from the CASdatasets package [1]. The raw data comprise two files: freMTPL2freq, containing 677,991 policy-year records, and freMTPL2sev, containing 26,444 individual claim amounts. These datasets provide a comprehensive set of risk features, including driver demographics, vehicle attributes, and geographic factors.

Table 1 summarizes the variables included in the frequency dataset. These covariates represent standard rating factors used in automobile insurance pricing, including driver age, vehicle characteristics, geographic region, and exposure.

The severity dataset contains 26,444 claim-level observations and includes the policy identifier and the corresponding claim amount for each reported claim. Because the frequency data are recorded at the policy level while the severity data are recorded at the claim level, the two datasets must be reconciled before modeling.

2.1. Data Integration and Preprocessing

To reconcile the policy-level frequency data with the claim-level severity data, all individual claim amounts in freMTPL2sev were aggregated by policy identifier (IDpol). This produced, for each policy, the total claim amount incurred during the exposure period and the number of associated claims. Policies with no reported claims were assigned a total loss of zero. The aggregated severity information was then merged with the frequency dataset to create a unified policy-level file containing exposure, claim counts, total incurred losses, and all rating variables.

For policies with at least one claim, an average claim severity variable was also computed for descriptive purposes. Pure premium values were obtained by dividing the total incurred loss by the policy’s exposure. These derived quantities are used only for exploratory analysis in this section; the formal notation and modeling framework are introduced later in Section 3.

2.2. Exploratory Data Analysis

The raw insurance portfolio exhibits the extreme class imbalance and heavy-tailed distributions typical of motor liability risks. Figure 1 illustrates these characteristics using log-scaled axes to visualize the full range of the data.

The frequency distribution (Figure 1, left) is dominated by zero-claim policies, with a rapid decay in frequency as claim counts increase. Notably, the raw data contain rare observations of up to 16 claims per policy, appearing as isolated points in the extreme tail. Similarly, the claim severity and pure premium distributions (Figure 1, center and right) span several orders of magnitude. The severity plot reveals a significant concentration of claims around 1,000 EUR, but also shows an exceptionally long right tail reaching towards ${10}^6$ EUR, representing catastrophic losses.

These visualizations highlight the necessity of the preprocessing steps described in Section 2.3. Specifically, the extreme sparsity of high claim counts and the high-leverage outliers in the severity tail motivate the use of capping (winsorization) to ensure the numerical stability and generalizability of the frequency and severity models.

2.3. Data Cleaning and Preparation

This section details the preprocessing steps implemented to ensure model stability and mitigate the influence of extreme observations. These procedures are critical for both the Poisson frequency models and the Gamma severity models, which are sensitive to high-leverage outliers.

Exposure Filtering Policies with extremely low exposure are prone to producing artificial variance in annualized claim rates. We excluded all records with an exposure below 0.1 years. This threshold ensures that the observations used for training represent a meaningful period of risk. As shown in Table 2, this step improved the mean exposure from 0.529 to 0.632 years.

Capping and Winsorization Given the extreme right-skewness observed in Figure 1, we applied capping to both frequency and severity components:

• Claim Frequency: Observations with more than 4 claims were capped at 4. Although the raw data contained counts as high as 16, these represented a negligible fraction of the portfolio ($<0.01\%$) but could disproportionately influence the maximum likelihood estimation.
• Claim Severity: Individual claim amounts were winsorized at 100,000 EUR. This prevents catastrophic “black swan” events—such as the observed maximum loss of 4.07 million EUR—from distorting the Gamma GLM parameters and the subsequent pure premium calculations.

2.4. Final Modeling Dataset Profile

The resulting dataset provides a stabilized foundation for estimating expected losses. Table 3 summarizes the response variables that will be utilized in Section 3. The disparity between the mean and median values across all metrics reinforces the inherent skewness that persists even after cleaning, necessitating the specialized modeling framework proposed in this study.

3. Methodology

This study investigates alternative approaches to modeling insurance pure premium by comparing the classical frequency–severity decomposition with direct modeling of aggregate claim amounts. Let $Y_i$ denote the total claim amount for policy i, $e_i>0$ the exposure, and ${\mathbf{x}}_i$ the associated covariates. The target of interest is the pure premium,

$${\pi }_i={\displaystyle \frac{\mathbb{E}[Y_i\mid {\mathbf{x}}_i]}{e_i}}.$$

3.1. Frequency–Severity Decomposition Framework

The classical actuarial framework expresses expected aggregate loss as:

$$\mathbb{E}[Y_i\mid {\mathbf{x}}_i]=e_i{\lambda }_i{m}_i,$$

where ${\lambda }_i$ denotes the expected claim frequency per unit exposure and $m_i$ denotes the expected claim severity conditional on at least one claim. Dividing by exposure, the pure premium is given by

$${\pi }_i={\lambda }_i{m}_i.$$

3.2. Claim Frequency Models

We evaluate multiple models for claim frequency and select the best-performing models for use in the decomposition framework.

3.2.1. Poisson Generalized Linear Model

The number of claims $N_i$ is assumed to follow a Poisson distribution:

$$N_i\sim \mathrm{Poisson}\left({\mu }_i\right),$$

with mean ${\mu }_i=e_i{\lambda }_i$. Using a log-link function:

$$log\left({\mu }_i\right)=log\left(e_i\right)+{\mathbf{x}}_i^{\top }\beta ,$$

where $log\left(e_i\right)$ is included as an offset.

3.2.2. Negative Binomial Model

To account for overdispersion, the Negative Binomial model is considered:

$$\mathrm{Var}\left(N_i\right)={\mu }_i+\alpha {\mu }_i^2,$$

with the same log-link structure as the Poisson model.

3.2.3. Poisson GLM with Natural Cubic Splines

To capture nonlinear relationships, selected predictors are modeled using natural cubic splines:

$$log\left({\mu }_i\right)=log\left(e_i\right)+{\beta }_0+\sum _{j=1}^4{f}_j\left(z_{ij}\right)+\sum _{k=1}^p{\gamma }_k{x}_{ik}.$$

3.2.4. Generalized Additive Model (GAM)

The GAM extends this approach by estimating smooth functions:

$$log\left({\mu }_i\right)=log\left(e_i\right)+{\beta }_0+\sum _{j=1}^4{s}_j\left(z_{ij}\right)+\sum _{k=1}^p{\gamma }_k{x}_{ik}.$$

3.2.5. XGBoost for Frequency

To capture complex nonlinearities and interactions, we apply XGBoost:

$${\widehat{\mu }}_i=\sum _{t=1}^T{f}_t\left({\mathbf{x}}_i\right),f_t\in \mathcal{F}.$$

3.2.6. Model Selection for Frequency

Based on predictive performance and risk segmentation metrics, the two best-performing models—Poisson GLM with splines and XGBoost—are selected for use in pure premium estimation.

3.3. Claim Severity Models

We consider two alternative approaches for modeling claim severity.

3.3.1. Gamma GLM

Claim severity $Z_i$ is modeled using a Gamma distribution with log-link:

$$log\left({\theta }_i\right)={\mathbf{x}}_i^{\top }\varphi .$$

3.3.2. XGBoost for Severity

To capture nonlinear relationships, we apply XGBoost to model the logarithm of severity:

$$log\left(Z_i\right)=\sum _{k=1}^K{f}_k\left({\mathbf{x}}_i\right)+{\epsilon }_i.$$

3.4. Pure Premium Estimation

The estimated pure premium under the decomposition approach is:

$${\widehat{\pi }}_i^{\mathrm{decomp}}={\widehat{\lambda }}_i{\widehat{m}}_i.$$

We consider four model combinations:

• Spline–Gamma: Poisson spline frequency + Gamma severity,
• Spline–XGBoost: Poisson spline frequency + XGBoost severity,
• XGBoost–Gamma: XGBoost frequency + Gamma severity,
• XGBoost–XGBoost: XGBoost frequency + XGBoost severity.

3.5. Direct Modeling of Aggregate Loss

To provide a benchmark, we consider models that directly predict aggregate claim amounts.

3.5.1. Tweedie Model

The Tweedie distribution satisfies:

$$\mathrm{Var}(Y_i\mid {\mathbf{x}}_i)=\varphi {\mu }_i^p,1<p<2,$$

corresponding to a compound Poisson–Gamma model. The mean is modeled as:

$$log\left({\mu }_i\right)=log\left(e_i\right)+{\mathbf{x}}_i^{\top }\beta .$$

3.5.2. XGBoost for Total Loss

We also apply XGBoost to directly model total claim amounts:

$$log(1+Y_i)=f\left({\mathbf{x}}_i\right).$$

3.6. Model Evaluation

Model performance is evaluated using both accuracy and risk segmentation metrics.

Prediction accuracy is assessed using Mean Squared Error (MSE), Mean Absolute Error (MAE), and Pearson correlation.

Risk segmentation is evaluated using decile-based analysis, where policies are ranked by predicted values and grouped into ten deciles. A well-performing model should exhibit strong concentration of losses in the highest-risk deciles.

4. Results

This section presents the empirical results for the claim frequency, claim severity, and pure premium models. We begin by evaluating the performance of alternative models for claim frequency, followed by claim severity, and then assess the combined frequency–severity decomposition and direct modeling approaches.

4.1. Claim Frequency Model Performance

Table 4 summarizes the performance of the claim frequency models.

Among all models, XGBoost achieves the best performance across all evaluation metrics, including the lowest MSE and MAE, as well as the highest correlation and lift. In particular, the lift in the top decile indicates a strong ability to identify high-risk policies, highlighting the effectiveness of machine learning in risk segmentation.

Among the parametric models, the Poisson model with spline terms and the GAM provide modest improvements over the standard Poisson GLM, suggesting that incorporating nonlinear effects enhances predictive performance. In contrast, the Negative Binomial model performs similarly to the Poisson GLM, indicating that overdispersion does not substantially affect predictive accuracy in this dataset.

Overall, these results demonstrate that machine learning models provide superior predictive performance and improved risk differentiation compared to classical parametric approaches.

Claim Capture Curves

Figure 2 presents the cumulative claim capture curves for the competing models. The XGBoost model consistently captures a larger proportion of claims in the highest-risk segments, followed by the Poisson model with spline terms and the GAM. The standard Poisson and Negative Binomial models exhibit weaker concentration of claims.

These results are consistent with the lift statistics and confirm that models incorporating nonlinear effects improve risk segmentation, with machine learning providing additional gains.

Decile Lift Analysis

Figure 3 presents the decile lift curves for the competing models. The XGBoost model achieves the highest lift in the top decile, followed by the Poisson model with spline terms and the GAM. Differences across the middle deciles are relatively small, reflecting the low overall claim frequency and limited variability in moderate-risk groups.

Interpretation of Model Components

Variable importance results from the XGBoost model (Figure 4) indicate that the bonus–malus score is the dominant predictor of claim frequency, followed by driver age, population density, and selected vehicle characteristics. These findings are consistent with actuarial intuition and highlight the importance of past claims experience and demographic factors in risk assessment.

The smooth functions estimated by the GAM (Figure 5) reveal nonlinear effects for several continuous predictors. In particular, the bonus–malus variable exhibits a strong nonlinear relationship with claim frequency, while vehicle age and population density display more moderate nonlinear patterns.

Model Selection

The results indicate that XGBoost provides the strongest overall predictive performance, achieving the lowest error measures, the highest correlation, and the greatest lift in the highest-risk decile. In addition, the Poisson model with spline terms delivers competitive performance while retaining the interpretability and transparency of the generalized linear modeling framework widely used in actuarial practice.

To enable a direct comparison between a traditional actuarial approach and a modern machine learning method in subsequent analyses, both models are retained. The Poisson spline model serves as the representative GLM-based specification, while XGBoost represents the machine learning approach. This dual-model strategy allows us to evaluate not only predictive accuracy but also the trade-offs between interpretability and performance in insurance pricing applications. These findings justify the selection of XGBoost and spline-based GLMs for subsequent pure premium modeling.

4.2. Claim Severity Model Performance

Table 5 summarizes the performance of the claim severity models.

The XGBoost model substantially outperforms the Gamma GLM across all evaluation metrics. In particular, it achieves a markedly lower MSE and MAE, along with a higher correlation with observed claim severity. These results indicate that machine learning methods are more effective in capturing the complex nonlinear relationships and interactions inherent in claim size data.

The Gamma GLM, while appropriate for modeling positive and right-skewed outcomes, assumes a log-linear relationship between covariates and severity. This assumption limits its ability to capture nonlinear effects and interactions, leading to inferior predictive performance.

Despite the improvement provided by XGBoost, the correlation remains relatively low, reflecting the inherent variability and stochastic nature of claim severity. This suggests that, although machine learning improves predictive accuracy, a substantial portion of severity variation remains unexplained.

These findings help explain the improved performance of decomposition models that incorporate XGBoost for severity, particularly in terms of prediction accuracy.

Figure 6 illustrates the relationship between actual and predicted claim severity. The Gamma GLM produces highly concentrated predictions, indicating strong shrinkage toward the mean and limited ability to capture variability in claim sizes. In contrast, the XGBoost model exhibits greater dispersion and improved alignment with observed values, highlighting its ability to capture nonlinear relationships and complex interactions.

However, both models display considerable scatter, underscoring the inherently stochastic nature of claim severity and the difficulty of accurately predicting individual claim amounts. Observations with zero or near-zero values were excluded from the log-scale visualization to avoid numerical issues associated with logarithmic transformation.

4.3. Results for Frequency–Severity Decomposition for Pure Premium

Table 6 presents the performance of the frequency–severity decomposition models using both error-based and ranking-based metrics.

Among the models considered, the XGBoost–XGBoost configuration achieves the best overall prediction accuracy, with the lowest MSE and MAE, as well as the highest correlation. This indicates that incorporating machine learning methods for both frequency and severity improves pure premium estimation.

In contrast, the ranking-based results reveal a different pattern. The XGBoost–Gamma model achieves the highest lift, indicating superior performance in identifying high-risk policies. This suggests that improvements in modeling claim frequency have a greater impact on risk segmentation than enhancements in severity modeling.

The difference between error-based and lift-based metrics highlights an important trade-off. While the XGBoost–XGBoost model minimizes prediction error, the XGBoost–Gamma model is more effective in concentrating losses within the highest-risk segment. This indicates that different model configurations may be preferred depending on the objective, such as pricing accuracy versus risk classification.

The classical Spline–Gamma model yields the weakest performance across all metrics, although it remains a useful benchmark due to its interpretability.

Overall, these results demonstrate that machine learning improves predictive performance within the decomposition framework, but its impact differs between accuracy and risk segmentation.

Risk Segmentation Performance

Figure 7 illustrates the decile lift curves for the decomposition models. All models exhibit a strong concentration of losses in the highest-risk decile, indicating effective identification of high-risk policies.

Among the models, the XGBoost–Gamma configuration achieves the highest lift, confirming its superior performance in risk segmentation. In contrast, although the XGBoost–XGBoost model provides the best overall prediction accuracy, its lift is lower, reinforcing the trade-off between prediction accuracy and risk ranking performance.

These findings suggest that claim frequency plays a dominant role in risk differentiation, while claim severity contributes more to overall prediction accuracy. Consequently, improvements in frequency modeling have a greater impact on identifying high-risk policies than enhancements in severity modeling.

The non-monotonic behavior observed in lower-risk deciles reflects the inherent variability of insurance losses, particularly for policies with low exposure. These results provide empirical evidence supporting the continued use of frequency–severity decomposition, particularly when combined with machine learning methods.

4.4. Comparison Between Decomposition and Direct Modeling for Total Loss

To evaluate the effectiveness of direct modeling approaches, we compare the best-performing decomposition model (XGBoost–XGBoost) with direct XGBoost and Tweedie models for total claim amount prediction.

Table 7 summarizes the performance of decomposition and direct modeling approaches.

The decomposition-based XGBoost–XGBoost model consistently outperforms both direct approaches across all evaluation metrics. It achieves the lowest MSE and MAE, as well as the highest correlation, demonstrating the effectiveness of modeling claim frequency and severity separately.

The direct XGBoost model provides competitive performance, indicating that flexible machine learning methods can capture nonlinear relationships in aggregate losses. However, it does not match the accuracy of the decomposition framework, suggesting that modeling aggregate losses directly may overlook important structural information captured by the frequency–severity approach.

The Tweedie model, despite its theoretical foundation as a compound Poisson–Gamma model, exhibits weaker performance across all metrics. This indicates that, although the data exhibit a compound structure, modeling aggregate losses within a single unified framework may be less effective than decomposing the problem into separate components.

The estimated Tweedie parameter ($\widehat{p}=1.52$) confirms that the data follow a compound Poisson–Gamma structure. However, the inferior predictive performance of the Tweedie model suggests that imposing a single functional relationship between covariates and aggregate loss may be overly restrictive. In contrast, the decomposition framework allows for distinct covariate effects on frequency and severity, resulting in greater modeling flexibility and improved predictive accuracy.

Overall, these results highlight the importance of structural decomposition in actuarial modeling. While machine learning enhances predictive performance within individual components, explicitly separating frequency and severity yields superior results. These findings provide empirical support for the continued use of the frequency–severity decomposition framework in modern insurance pricing.

5. Conclusions and Discussion

This study examines the integration of classical actuarial models and modern machine learning techniques for insurance pricing, with particular emphasis on comparing the traditional frequency–severity decomposition framework with direct modeling approaches. The results provide important insights into predictive performance and practical implications for risk management.

The empirical findings show that machine learning methods, particularly XGBoost, substantially improve predictive accuracy for both claim frequency and claim severity. The improvement is especially evident for severity modeling, where the Gamma GLM tends to shrink predictions toward the mean and fails to capture the variability of claim sizes. In contrast, XGBoost effectively models nonlinear relationships and interactions, leading to lower prediction errors. Despite these improvements, correlation values remain relatively low, reflecting the inherently stochastic nature of claim severity.

Within the decomposition framework, the XGBoost–XGBoost model achieves the highest overall predictive accuracy, confirming that machine learning enhances pure premium estimation when applied to both components. However, improvements in error-based metrics are moderate, which highlights the high variability of insurance losses, particularly for policies with low exposure.

An important contribution of this study is the identification of a trade-off between prediction accuracy and risk segmentation. Although the XGBoost–XGBoost model minimizes prediction error, the XGBoost–Gamma model achieves the highest lift and provides superior identification of high-risk policies. This result indicates that claim frequency plays a dominant role in risk differentiation, while claim severity contributes more to overall prediction accuracy. Consequently, improvements in frequency modeling have a greater impact on ranking policies by risk level than enhancements in severity modeling.

The comparison between decomposition and direct modeling approaches provides further insight. While direct XGBoost and Tweedie models deliver competitive performance, they do not outperform the decomposition-based framework in this study. This finding contrasts with the results of [17], which suggest that direct modeling of aggregate losses can yield lower prediction error.

This discrepancy highlights an important distinction between general predictive modeling and actuarial applications. Direct modeling may perform well in settings where minimizing prediction error is the sole objective. However, in insurance pricing, the frequency–severity decomposition reflects the underlying structure of claim processes and enables separate modeling of distinct risk components. This structural advantage allows for greater flexibility in capturing heterogeneous risk factors and improves risk segmentation.

In addition, direct modeling imposes a single functional relationship between covariates and aggregate loss, which may be overly restrictive when frequency and severity are influenced by different drivers. In contrast, the decomposition framework accommodates distinct covariate effects for each component, resulting in improved predictive performance in practice.

From a practical perspective, these findings have important implications for insurance pricing and risk management. If the objective is to minimize prediction error, models that incorporate machine learning for both frequency and severity are preferred. If the goal is to identify high-risk policies for underwriting or portfolio management, models that emphasize frequency modeling are more effective.

The results are also consistent with prior studies comparing gradient boosting methods with classical GLMs within the frequency–severity framework [19]. As in previous work, machine learning significantly improves claim frequency prediction. However, while earlier studies report mixed results for severity modeling, the present analysis shows that XGBoost can outperform the Gamma GLM when appropriately specified and tuned. This suggests that the relative performance of machine learning for severity is context-dependent and may improve with richer data and more flexible model configurations.

Overall, this study demonstrates that machine learning enhances predictive performance, but its effectiveness is maximized when combined with the classical frequency–severity decomposition framework. The results indicate that theoretical advantages of direct modeling do not necessarily translate into superior performance in real-world actuarial applications. These findings provide empirical support for decomposition-based pricing approaches and offer a foundation for future research on hybrid models that balance predictive accuracy and risk segmentation in complex insurance environments.