Poisson Incorporated Credibility Regression Modelling of Systematic Mortality Risk for Populations with Finite Data

This study considered the modeling of systematic mortality risk for populations with finite data using the Poisson incorporated Credibility regression model. For novelty, we have included the credibility regression approach to modelling mortality by assuming the number of annual deaths follow a Poisson distribution. Our model shows improvement in precision levels when estimating mortality risk compared to classical models used in European countries. We have illustrated that our model works optimally when using Kenyan mortality data, comparing male and female lives under the different strategies, thus making better predictions than the classical Lee–Carter (LC) and Cairns–Blake–Dowd (CBD) models. The mean absolute forecast error (MAFE), mean absolute percentage forecast error (MAPFE), root mean square error (RMSE), and root mean square forecast error (RMSFE) under the incorporated credibility regression model are much lower than the values obtained without incorporation of the Buhlmann credibility approach. The findings of this research will help insurance companies, pension firms, and government agencies in sub-Saharan countries model and forecast systematic mortality risks accurately. Finally, the results are essential in actuarial modelling and pricing, thus making life assurance products affordable for most people in low-income African countries.


Introduction
e concept of stochastic mortality risk modelling was introduced by Lee and Carter [1], making it possible to model both time and age as random variables. Moreover, mortality and longevity risk become random variables that change with age and time. is new approach led to more accurate mortality modelling, thus enabling researchers to get precise results that correspond to the actual demands of the market [2]. However, the primary source of modern risk faced by companies, pension funds, and life companies is longevity. An individual who has taken out life cover will tend to live longer than the expected age calculated by an actuary, thus posing a risk for those rms o ering actuarial services. e uncertainty posed by mortality and longevity risk has been further compounded by increased life expectancy in the 21st century.
ere has been an increase in life expectancy by over 10 years due to various factors including improved medical services, diet, and increased awareness of life choices. e increase in life expectancy has been experienced especially in developing countries, e.g., Kenya. erefore, exposure to extreme mortality risk is a crucial issue. However, life companies, governments, and pension rms are yet to formulate new approaches of dealing with or managing these risks. In addition, life companies, governments, and pension funds are yet to innovate ways of coping with or managing systematic mortality risks especially in light of the prevailing global modern nancial challenges [3]. e systematic mortality risk poses a massive challenge for the companies and governments that want to avoid insolvency from poor nancial decisions that might occur whilst they operate in the financial markets. However, many extensions of the model [1] have been made, which include [2,4], among many others.
In many African developing countries such as Kenya, with limited data points, it may be difficult to trace, analyze, and make meaningful inferences when modelling mortality and longevity risk factors. However, the concept of credibility theory, primarily used in nonlife actuarial valuations, can help model mortality and longevity risks. Many researchers use deterministic models for instance in [5], which do not consider the population's stochastic nature of death rates. e model's limitations is the difficulty in parameter estimation due to insufficient data points. e concept of credibility theory [6,7] was introduced for mortality and longevity risk modelling by Hardy and Panjer [8]. is concept has been instrumental when dealing with insufficient data in a population such as those from many developing countries, especially in Africa. is study applies a novel approach to modelling mortality using a Poisson distribution but under an incorporated credibility regression, especially when dealing with data limitations [9][10][11].
e study assumes that the established and limited populations exhibit similar longevity trends and cohort mortality levels instead of the model proposed in [12]. In our applied model, the force of mortality has inverse proportional trends as the population with limited data increases, meaning that as the credibility factor increases, this reduces our over-reliance on an established population when estimating longevity and forecasting mortality from the conventional models. e Poisson regression model is widely applied in the study of real data problems, such as mortality studies, where the aim is to investigate the number of deaths. In addition, in health insurance, where the target is to explain the number of claims made by individuals [13][14][15][16] illustrated how to use the Poisson regression model when dealing with the response variable consisting of count data. Furthermore, [14] demonstrated the almost unbiased ridge estimator in the zero-inflated Poisson regression model, which can be applied to the mortality study for an infinite number of deaths.
Most of the recent actuarial studies on modelling mortality under a framework of credibility theory were made in several studies, including [10,[17][18][19] in their published papers. e authors first applied Bühlmann's credibility theory to the mortality data of three modern developed countries: Japan, the UK, and the US. In contrast, in the second paper, they applied the concept of Bühlmann credibility to the three conventional models, namely, the model proposed in [1], the model proposed in [20], and the linear relational model proposed in [21] as a way of improving forecasting performances for the UK dataset.
We conduct our valuation using the Kenyan mortality data and compare the obtained values with the tabulated values from the same classical models for illustration purposes. Furthermore, we illustrate how our precision increases as we increase the data that affect the credibility factor, thus reducing the reliance on the established data [12].
is study helps us estimate the correct number of deaths and lives living beyond the expected life expectancy, which is vital for planning for actuaries, life assurance companies, governments, and pension firms.
Nevertheless, it has been noted that age-specific mortality displayed a downward trend with time. Additionally, suppose a case with a finite mortality data experience for a given age, but an elaborate data experience for the whole age range is provided. In that case, one can apply Poisson incorporated credibility regression techniques to capture the mortality trends [22]. e research uses the Poisson credibility regression approach as a complementary approach to Bühlmann's credibility as [10] when forecasting the mortality rates [4] for developing African countries such as Kenya, with limited data experience.
Finally, Odhiambo et al. [23] proposed using Bühlmann's credibility approach to model systematic mortality risk, especially for sub-Saharan African populations with limited mortality data, such as Kenya. In addition, Odhiambo et al. [24] proposed incorporating the deep learning technique to the Cairns-Blake-Dowd (CBD) model while modelling systematic mortality risk to enhance the accuracy of the models used in cases with limited data for systematic mortality risk modelling and forecasting. ese methods improve the accuracy and robustness of the models used in SMR forecasting, thus enhancing the valuation and pricing of actuarial products in Kenya.
e study is organized as follows: in Section 2, we briefly give reviews of the conventional models [1,20] and the random coefficients regression models. In Section 3, we propose the Poisson-incorporated credibility regression approach as a complementary approach to Bühlmann's credibility when forecasting the mortality of a given population with finite data for novelty purposes instead of conventional Poisson modelling of deaths.
In Section 4, we present the extrapolation methods used to estimate future mortality rates under the Poisson incorporated credibility regression approach. We also conduct an empirical illustration using Kenyan data, both males and females. Furthermore, we evaluate the forecasting performances of the Poisson-incorporated credibility regression approach using various methods proposed in [1,20], e.g., MAFE, MAPFE, and RMSFE measures. Finally, we use Akaike Information Criteria (AIC) and Bayesian Information Criteria (BIC) to evaluate and select the best model.
In Section 5, we apply the Bühlmann credibility approach, which is essential in modelling and forecasting the mortality data, especially for countries with limited data in many developing countries globally, mainly in African countries such as Kenya.
Finally, we conclude in Section 6, noting that using the Bühlmann credibility approach application can help improve modelling precision whenever there is a problem with mortality or life table data availability. In the appendix, we have systematic mortality risk values comparing different ages.  [8] has been effective when modelling different mortality tables for different countries regarding their demographic characteristics. e ages of a population can be from time x to time t, where t � 0, 1, 2, 3 . . ..

Definition 1.
Let D x+t denote the number of people who die with time from time x and △(x + t) that follow a conditional Poisson distribution of where the values of x, t � 1, 2, 3 . . . for all x < t. e value e(x, t) denotes the population surviving past a given age from x to x + t exposed to death and μ(x, t) is defined as the force of mortality of the established population from age x to x + t. e value will be estimated from the established models [20] using poorly established data from a country such as Kenya. It is essential to assume that the rates will change from one discrete age to another depending on an individual's age, as in [25].
We use two mortality models to model the value of μ(x, t), which will be multiplied by the number of people exposed to calculate the total number of deaths. e results are what will have been experienced within a given population with limited data sets, such as in the case of the Kenyan demographic setup.
e Poisson distribution assumes that the annual deaths are uniformly distributed throughout the study year, called the UDD assumption. ere is a constant force of mortality known as the CFM assumption between two integer ages when modelling and forecasting using the mortality models.

Mathematical Preliminaries of Stochastic Mortality Models.
is section briefly reviews the two classical mortality models, namely, the models proposed in [1,20], which are applied within the research study in the subsequent subsections.

e Lee-Carter Mortality Model.
In the original model [1], the logarithm links of the year recorded mortality rates denoted as M(x, t) � ln(m(x, t)) for all ages from x � 0, 1, 2, . . . , n − 1, and a predictor model is given by the following equation: where α x is defined as the age parameter reflecting the mean mortality at an exact age of x, k t is the period parameter indicating the overall mortality trends in year t and β x is the age parameter indicating the differences from the mean mortality at age x, with changes in the overall level of mortality index.
e error term is defined as ε(t, x) and is assumed to be normally distributed with mean zero, and a constant variance, σ 2 . In addition, it reflects a specific period as well as age effects that have not been captured by the mortality model. erefore, after assuming that the randomness of the model is independent, identical, and homoscedastic with a mean of zero, as in [1], a way of approximation is suggested using the method of SVD (singular value decomposition), which is subjected to the constraints K−1 x�0 β x � 1 and K−1 x�0 k t � 0 when obtaining the model estimates of the parameters. e estimates are as follows: Upon selection of the SVD method of estimation [1], we often apply the use of the time series models when projecting mortality using the parameters. e values of k t are obtained using the following formula:

e CBD Mortality
Model. e original model [20] links the logit transformation of the single-year probabilities of mortality or deaths denoted by M(x, t) = logit q(t, x) with model predictor as follows: where k (1) t is defined as the period parameter indicating the overall level of mortality within year t and k (2) t is the period parameter showing how mortality index affects every age, whilst x is defined as the mean age of the selected fitting age interval. ε(t, x) is the specific random effects that have not been captured on the model [20], and is expected to follow a Gaussian distribution with a mean of zero and a constant standard deviation.
Similarly, as shown in equation (1), we briefly present the model parameter estimates that can be obtained easily by regressing the logit q(t, x) on the value of (x − x) for every t value: and the values of k (2) t are given as follows: According to [25], we made an assumption that number of deaths follows a Poisson distribution with a mean of Mathematical Problems in Engineering deaths rate and E(t, x) is the number of people exposed. However, we estimate the model parameters through the maximization of the log-likelihood function of equation (5).
With an assumption of period estimates being identical and independent for each other, the use of extrapolation using a random walk having a drift parameter (ϑ i , for i � 1, 2) and then mortality forecasts for the period of h � 1, 2, 3, . . . , H are obtained by the following equation:

The Credibility Regression of Mortality Risk Models
is section deals with a mortality modelling approach within a credibility regression framework with different coefficients. Also, it describes the procedure of parameter estimation while providing a particular case with fixed coefficients within it.

e Credibility Regression Structure under Random Coefficient (RC)
Definition 2. Let D(t, x) be the number of observed deaths at a given exact age x within t years' time and E(t, x) as the mean population exposure to risk at aged x within year t.
Furthermore, we define the ratio of D(t, x)/E(t, x) as m(t, x), which denotes the age-specific mortality rates as illustrated in [26].
denote the one-year probabilities of death, which is the assumption of the constant force of mortality (CFM) within the different individual age groups.
With the response variable assumption defined by M(x, t) as the transformation of natural logarithm and logit of measure of mortality, m(t, x) and q(x, t) respectively for ages x � 0, 1, 2, 3 . . . k − 1, for integer ages t � 1, 2, 3, . . . n − 1, and for calendar years are provided in [27]. We denote B x as a random age-related risk parameter, which is defined by . . M (t n +x) ) ′ as mortality index and W x is the design matrix of the applied explanatory variables under [28]. us, the pair describing the evolution of mortality in age x is given by (B x , Mx), under these assumptions: (1) e pairs of (B x0 , Mx 0 ), (B x1 , Mx 1 ), . . . , (B xk−1 , Mx k−1 ) are independent while the pairs of (B x1 , B x2 , . . . , B xk−1 ) are independent and identically distributed.
where W x is the fixed nxq design matrix having a full rank of q( < n) and ϕ (B x ) is an undefined regression vector having a length of p.
e structure of the parameters is defined in the following way: In the above regression setting, we can estimate △x. Consequently, we estimate the values of regression coefficients using the method of generalized least squares method or GLS. is gives an individual estimator of ϕ(B x ) to be

Credibility Regression Using Fixed Coefficients (FC).
When we fix the coefficients of regression and an assumption of some weights [29], it may appear on the regression line. For example, taking the number of people who are exposed as E(t, x), for values oft � 0, 1, 2, 3 . . . n − 1 that is then applied. We apply assumptions (1) and (2) outlined in the previous subsection as they are but with a modification of the assumption (3), where Cov(Mx0/B x ) � σ 2 (Mx)P x and P x is a square m × m positive diagonal matrix having the weights, [27]. Its structural parameters are then defined as follows: whereas the coefficients vector of ordinary least squares estimator of the ϕ(B x ) is defined as follows: From the above assumptions, it easy to note that the credibility estimator of the model for the fixed coefficients (FC) model is provided as follows: 4 Mathematical Problems in Engineering where known as the estimated credibility factor [29]. Also, we apply the following estimators: Besides, the values of the estimated parameters of a and υ have been determined through the process of iteration, while imposing a constraint of υ + υ ′ /2 � 0 after every process of iteration.

Remark 1.
We fitted for the random coefficient (RC) as well as the fixed coefficient (FC) of the credibility regression models when modelling mortality rates for a given set of age x � 0, 1, 2, . . . , k − 1 for the year t � 0, 1, 2, . . . , n − 1.
For both models [1,20], the fitted mortality rates go up to year n − 1 and can be wholly written as an estimate of Mx � ϕ x W x .

Extrapolation Methods for Estimating Future Mortality
Rates. When fitting current data of the response variable for the mortality data provided, M(x) � (M (t 0 +x) , M (t 1 +x) , M (t 2 +x) , . . . M (t n +x) ) ′ , which is the mortality rate for a given year ahead and can be estimated as follows: where c � RC or FC. According to proposals of estimation of a single-year ahead estimates [18], M t n−1 +1,x is always fixed to the current fitting span, without having to remove M t 0 ,x ; hence, the fitting year span can be extended yearly sys- Also, within each of the estimation steps, we fitted the credibility regression models directly on a continuously extended response variable to yield is method of mortality extrapolation of future mortality trends is based on both the first and the current rates, determined after each step of the estimation process [10]. It is important to note that numerical results in the Mathematical Problems in Engineering 5 following section justify the methods applied efficiently in actuarial practice, especially in designing and pricing life assurance products and annuities sold to policyholders.

Empirical Mortality Trend Demonstration.
is section fits the models proposed in [1,20] and the credibility regression models on Kenyan population mortality. After that, we evaluate the forecasting results using the two statistical rations, namely, mean absolute forecast error (MAFE) and the root mean of squared forecast error (RMSFE) measures.
Kenyan data have a limitation in historical mortality observation numbers (2010-2020) and are structured according to year, age, and sex. Moreover, in many life insurance datasets,      (3) similar limitations do currently exist. is credibility regression can proficiently get the trends of underlying mortality data, particularly in those cases where a limited mortality experience during a specific period for a specific age is experienced but an elaborate experience for the complete age range (Kenyan mortality data). It makes regression methods easy to apply in several cases when dealing with larger mortality tables of developed countries [30]. Mortality evolution for the period 2010-2020 in Kenya is illustrated in Figures 1 and 2 for log (m(t, x)) and logit (q(t, x)). Both systematic mortality risk measures show linearity for discrete-time ages x � 40, 60, and 80 of males and females (Figures 1 and 2) as well as females. Also, mean mortality decline has shown a vivid downward trend over time (right panels of Figures 1 and 2) for both males and females. Figures 1 and 2 show the increase in the observed values of logm(t, x) and logitq(t, x) for the 2010 period in Kenyan, for males, females, and combined males and females at the ages 20, 30, 40, and 60.

Remark 2.
ere is a general decrease in mortality between ages 25 and 30 before an upsurge of the trends in the remainder of the lifetimes of both males and females in Kenya.

Forecasting of Mortality Results.
For purposes of numerical illustration in the following section, we applied the empirical data of age-specific mortality rates of m(t, x) from 2010 to 2020 for both genders, such as males and females, from ages 20 years to 80 years. Besides, the age span choice as in similar studies [19], mainly when corresponding to a person's age up to the overall level of complete expectation in the developed countries.
To ensure robustness relative to the fitting range of data changes, we have applied two age and three-period spans when extracting forecasts for a 10-year (H � 10) forecasting horizon, as listed in Table 1. Specifically, in using the FC model, we used M y � I as proportions or weights. e illustration of credibility regression methods, the LC model and the CBD model, has been implemented in the R program to use Poisson LC and CBD models; while fitting these methods, we have used the "Life Metrics," and "StMoMo" packages in the R statistical programming language [31].
As it retains the concept of linearity over all the available fitting period, the logarithmic transform M t n−1 ,x is used for all the age-specific mortality rates and the logit for the given-year death probabilities.
When forecasting errors that were then calculated over the 10-year forecasting horizon when using MAFE and RMSFE, smaller values can indicate a much better performance when forecasting. We calculate the mean or average of MAFE and RMSFE values using the following equation: Similarly, in the case of using M t,x � logitq(t, x) as a given response variable, m(t, x) should be replaced by the q(t, x) as per equations (17) and (18), respectively.
When the forecast accuracy results at a percentage (%) scale, which are evaluated over the period 2010-2020. e values of MAFE and RMSFE values especially for fitting ages 20-80, which are using M t,x � logm(t, x), are listed in Tables 2  and 3, respectively; whereas, the corresponding values for ages 60-100 with M t,x � logm(t, x) are listed in Tables 4 and 5, respectively. It is important to note that the CBD model has been included only for comparisons when fitting ages 55-84, which has particularly been designed for the higher ages.
For both genders, males and females, the results are tabulated in Tables 2 and 3 when fitting ages 20-60 and  Tables 4 and 5 for ages 60-80, thus indicating that, for every fitting period, credibility regression models outperform LC and CBD models for both error measures.  e mean values are provided in the last rows of every measure's suitability over the whole period. To be more precise, for specific ages 15-84, the FC-MEM and FC-SEM produce, thus making the smallest mean of MAFE and RMSFE. In contrast, for ages 60-80, RC-MEM performs better in the mean on both measures, indicating that it predicts higher ages based on the recent mortality trends.

Model Selection.
From the two models, we will look at the best among the two in systematic mortality risk modelling. We do model selection using the log-likelihood and the number of the effective parameters along with AIC and BIC values of the mortality models for males and females. Table 6 lists the AIC and BIC values for both models. e AIC and BIC values for the two models show that Cairns-Blake-Dowd (CBD) model is the better model Using the Bühlmann credibility approach in mortality modelling is essential when the country is experiencing limited data for modelling and predicting systematic mortality risk. Ultimately, this is helpful during the correct pricing of the life assurance products sold in the market for the consumers.

Effects of Credibility Approach on Stochastic Mortality
Models. In the previous subsection, we used the proposed credibility regression methods when estimating the actual mortality trends for a specific age through the proportion of the mortality trends for this age and the mean trend over a much wider group of ages, which provide more information: Remark 4. Figures 3 to 5 show the intercept as well as slope estimates of logitq (2000 + h, x) for h � 1, 2, . . . , 10 and ages x � 60 for males and x � 80 for females, with credibility for FC-MEM and RC-MEM and without credibility for the mortality models, making the incorporated credibility more accurate in predicting the mortality models from the two mortality risk models, namely, Lee-Carter and CBD models.
On the effect of the credibility on the trends, we have the models listed below.
e Poisson incorporated credibility regression modelling of mortality, especially for the populations with finite data, can be determined from the Poison distribution as follows: Figure 6 shows the Poisson incorporated credibility regression modelling of mortality for Kenyan males, which shows how deaths of males have been trending over the last decade based on the two mortality models, namely, the Lee-Carter and CBD models. Figure 7 shows the Poisson incorporated credibility regression modelling of mortality for Kenyan females, which shows how deaths of males have been trending over the last decade based on the two mortality models, namely, the Lee-Carter and CBD models. Figure 8 shows the Poisson incorporated credibility regression modelling of mortality for Kenyan males and females. It also shows how deaths of males have been trending over the last decade based on the two mortality models, namely, Lee-Carter and CBD models.  life assurance products sold in the Kenyan market. e deaths have decreased from the data since people live longer than expected, thus there is a massive challenge for the financial institutions making these life annuities and assurance products.

Application of the Bühlmann
Credibility Approach [32] proposed a Bühlmann credibility approach that can be used when forecasting mortality rates for both males and females in the three most developed countries: Japan, the UK, and the US. Our model improves its model by modelling it as an incorporated Poisson distribution.

Definition 3.
Let the general regression model be Z y � (1, 1, . . . 1) ′ and W y � I for y � y 0 , y 1 , . . . y k−1 , then it is easy to estimate the parameters of the regression model as follows: where U and K are mortality index and Bühlmann credibility factor estimates, respectively.       e systematic mortality risk from the Bühlmann credibility estimates becomes all for the values of y � y 0 , y 1 , . . . y k−1 and Q y and Q are for the mortality risk of population estimates for finite and infinite data respectively. Contrary to the credibility regression approaches that aims at capturing the central death rate, m(t, x) downward trend as in s over t, application of the Bühlmann credibility approach eliminates it during modelling.
For the same reason, Tsai and Lin [19] did an application of the Bühlmann credibility model on a specified time series of mortality rate that changes as opposed to flat levels of the mortality rate, such that, ereafter, they applied two strategies during the estimation of Q t+h,x , h � 2, 3, 4, . . . , H. e first strategy expands the fitting window (EW strategy) by a year, similar to the EEM regression method, as described in Section 4, and the second one moves the fitting window (MW) through a year, similar to the MEM regression method.
Similarly, we have done a comparison of the predicting performance in between the Bühlmann and the credibility regression approaches on Kenyan data [32]. To be consistent with the Bühlmann modelling structure [19]. Also, the age fitting spans between 20-100 and 60-100, which were chosen, and forecast errors were also evaluated under the mean MAPFE values denoted by MAPFE mean as follows: It is important to note that the error values for every gender were evaluated by fitting Q t,x , s on period 2010-2020 and 2020-2050.
A comparison of mean values of MAFE, RMSFE, and MAPF results between Bühlmann and credibility regression methods is given for both genders in Table 4 for ages 20-60 and Table 5 for ages 50-90. e results indicate that credibility regression methods produce the smallest MAFE, RMSFE, and MAPFE values for the most selected fitting periods for both age spans. e FC-MEM method has the optimum expected performance per the MAFE and MAPFE values for ages 20-80. In contrast, the RC-MEM method seems more excellent in capturing future mortality trends for older ages 60-80.
We note that the smallest values in expected values are produced by different regression methods, which depend on the measures used. Such inconsistencies are always expected due to the nature of MAFE, RMSFE, and MAPFE formulas. e concept was discussed, as pointed out by Tsai and Yang [21].
Credibility regression techniques have proved vital for mortality datasets of a comparatively short historical time or period of recordings (for finite data). ey have efficiently captured the underlying mortality trends for a specific age from all the information gathered from populations of many other periods.   (3) is research study has proposed the mortality modelling approach that embeds it within the framework of credibility regression. Using Kenyan data for illustration, credibility regression approaches have led to better SMR forecasts for males and females. Comparing MAFE and RMSFE measures to the conventional LC and CBD models, they fit well during the credibility approach modelling.
Finally, we used a structure proposal of credibility regression mortality when dealing with randomly varying coefficients on a particular case with fixed coefficients. us, as explained in previous section five, we have presented extrapolation methods for each credibility approach when estimating the future mortality rates. From the analysis, we have concluded that, in aggregate, credibility modelling methods performed better than the LC and CBD mortality models. For forecasting accuracy during the study, we have seen a better analysis for the whole age fitting span, with fixed coefficients credibility methods that performed much better on average when compared to the RC-MEM, which makes a better analysis.

Conclusion and Practical Applications
e proposed credibility regression model have resulted in better SMR forecasts based on MAFE and RMSFE measures. Credibility regression thus offers excellent forecasting performance when applied to datasets of developing countries (having a relatively small population), such as Kenyan mortality data. As a proposal for further research, the forecasting comparison between many other countries' datasets would be a good avenue of explorative data analysis.
During the comparison, the Bühlmann credibility approach was applied to the Kenyan life table in Section 5. e credibility regression methods yielded excellent SMR forecasts based on MAFE, MAPFE, and RMSFE measures. Moreover, credibility regression methods had an excellent predicting performance.
We applied the proposed model to the datasets of other countries (with a relatively small population) for a few selected fitting periods for many countries with developed structures in mortality modelling. Additional prediction comparison between datasets of other countries has been left out for future research work. Our numerical illustration yielded relevant results and proved the credibility modelling approach, including those that could contribute to future systematic mortality risk projection studies.
On policy recommendation, the government and policymakers should model the systematic mortality risks correctly since it will always determine the amount of money in terms of benefits payable to the policyholders whenever they are in financial need from the life assurance and annuities bought. Ensuring that the affordable premiums are payable makes it an avenue for many policyholders who may now purchase these products sold by the pension, government, and insurance companies in developing countries. e results of this study will help better systematic mortality risk modelling, which is important in actuarial modelling and pricing. Many pension and insurance firms will improve their estimates, which ultimately helps enhance the pricing rates of the actuarial rates sold in the Kenyan markets. By making these premiums affordable, policyholders can buy more products while reducing the burden of deaths that many Africans experience while living in many sub-Saharan countries such as Kenya.  show the SMR values of MAFE, RMSFE, and MAPFE measures of forecast errors during the different periods for ages for both males and females when modelling the mortality rates. e authors note that incorporating the Bühlmann credibility approach helps enhance the accuracy and prediction of future mortality risks from the lower levels of errors than classical models. Table 7 lists the MAFE values of forecast errors during the period 2010-2020 for ages 20-60 for the two models under the credibility approach. Table 8 lists the RMSFE values of forecast errors during the period 2010-2020 between the ages 20 and 60 for both males and females. Table 9 lists the MAPFE values of forecast errors during the period 2010-2020 for ages 20-60 for males and females. It shows how high levels of precision are achieved under the credibility approach. Table 10 lists the MAFE values of forecast errors during the period 2010-2020 for ages 20-80 for males and females. Table 11 lists the RMSFE values of forecast errors during the period 2010-2020 for ages 60-100 for both males and femalesTable 13. Table 12 lists the MAPFE values of forecast errors during the period 2010-2020 for ages 60-100 for both males and females when modelling the systematic mortality rates.

Data Availability
e data used to support the findings of this study can be obtained from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.