BAYESIAN PARAMETER ESTIMATION FOR TRANSMUTED WEIBULL DISTRIBUTION WITH CENSORING RATES AND VARIOUS LOSS FUNCTIONS
1Jeslin J, 2Radhika A*, 3Haripriya M
1Research Scholar, Department of Statistics, Periyar University, Salem - 11 [email protected] 2Assistant Professor, Department of Statistics, Periyar University, Salem - 11 [email protected] 3Research Scholar, Department of Statistics, Periyar University, Salem - 11 [email protected]
Abstract
Statistical distributions are essential tools for describing and predicting real-world phenomena, though recent advancements in data collection have made it challenging to fit existing probability models to many practical datasets. While non-parametric models are sometimes recommended, parametric models retain substantial popularity due to their interpretability and flexibility. The quadratic rank transmutation map (QRTM) technique has been used to create new families of non-Gaussian distributions, known as transmuted distributions, which allow for modifications in moments, skewness, and kurtosis, thus increasing flexibility. The transmuted Weibull distribution (TWD) has gained attention for applications in reliability, survival analysis, and lifetime data analysis. This article focuses on a Bayesian analysis of the transmuted Weibull distribution, a generalization of the traditional Weibull model that addresses its limitations, particularly for datasets exhibiting non-monotonic failure rates. Bayesian parameter estimation is performed using a Markov Chain Monte Carlo (MCMC) algorithm, with both non-informative and informative priors. We calculate Bayes estimators (BEs) and posterior risks (PRs) under different loss functions, including the Absolute Error Loss Function (AELF), precautionary loss function (PLF), and quadratic loss function (QLF). Simulation studies evaluate the Bayes estimators' performance, investigating the effects of various priors, sample sizes, and censoring rates on estimation accuracy and credible interval width. Real-world data applications highlight the practical utility of the Bayesian approach for the TWD, showing consistent results with increasing sample sizes and underscoring the robustness of the MCMC algorithm for parameter estimation. The article is structured as follows: the TWD's parameters, including scale, shape, and transmutation, are estimated under different loss functions and priors. Bayesian credible intervals (BCIs) are also computed. Both uncensored and censored data environments are considered, with varying sample sizes and censoring rates. Posterior risks for each estimator are analyzed to assess performance, and two real datasets are used to illustrate the flexibility and applicability of the proposed distribution. This study lays a foundation for future research, such as exploring mixtures of transmuted Weibull distributions or conducting Bayesian analyses for record values.
Keywords: Transmuted Weibull distribution, Markov Chain Monte Carlo, Bayesian credible intervals, Bayes estimators, posterior risks, absolute error loss function, precautionary loss function, quadratic loss function, censoring.
I. Introduction
Statistical distributions are essential for describing and predicting real-world phenomena. However, advancements in data collection methods have led to challenges in fitting existing probability models to many significant and practical datasets [1]. In such cases, while non-parametric models may be recommended, parametric models continue to enjoy substantial popularity. The quadratic rank transmutation map (QRTM) technique has been employed to create new families of non-Gaussian distributions [2]. This technique modifies the moments, skewness, and kurtosis of a baseline distribution, resulting in what is known as a transmuted distribution. This family of distributions has attracted considerable attention from researchers, leading to the development and exploration of various new flexible distributions over the past decade. For instance, A1 Sobhi [3] introduced the transmuted modified Weibull distribution, while others presented the exponentiated transmuted Weibull distribution. More recently, a new lifetime distribution called the transmuted cubic Weibull distribution was constructed, and a novel weighted distribution known as the size-biased weighted transmuted Weibull distribution was introduced. The method of least squares and the method of moments have been utilized to estimate parameters for the transmuted Weibull distribution, with comparisons made through simulation studies using statistical measures like mean squared error (MSE) [4]. Researchers have also explored various structural properties of the transmuted Weibull distribution, including its mean, harmonic mean, standard deviation, moment generating function (MGF), skewness, and kurtosis [5][6]. Currently, transmuted distributions find applications in numerous fields, including reliability studies, lifetime analysis, engineering, economics, insurance, and environmental sciences.
The Weibull distribution is a widely recognized lifetime probability distribution, commonly used in various domains of reliability and survival analysis [7]. Its attractiveness largely arises from the variety of shapes it can assume by adjusting its shape parameter. Additionally, the Weibull distribution is versatile, encompassing many other distributions, such as the exponential and Rayleigh distributions, as special cases. Despite its widespread use and applicability, the traditional Weibull distribution does not fully capture the range of lifetime phenomena. For instance, it is not suitable for datasets exhibiting a non-monotonic failure rate, prompting investigations into generalizations of the Weibull distribution. A notable generalization applicable to survival data analysis is the transmuted generalized inverse Weibull distribution, which discusses its mathematical properties and employs maximum likelihood methods for parameter estimation. Similarly, the transmuted generalized Weibull distribution has been developed, exploring its mathematical properties, including the quantile function, moments, entropies, mean deviation, Bonferroni and Lorenz curves, and the moments of order statistics, also using maximum likelihood for parameter estimation. Furthermore, the generalized transmuted Weibull distribution has been proposed, with its properties derived. This article focuses on the Bayesian analysis of the transmuted Weibull distribution, which serves as a generalization of the Weibull probability distribution. We emphasize Bayesian analysis because maximum likelihood and moment estimators have been used for parameter estimation of the transmuted Weibull distribution [8]. To facilitate this analysis, we employ a Markov Chain Monte Carlo (MCMC) algorithm to compute posterior summaries for the unknown parameters of the distribution, comparing the results across different priors, loss functions, sample sizes, and parameter sets [9].
The objective of this paper is to define the Transmuted Weibull distribution and introduce its likelihood function, followed by the derivation of posterior distribution expressions utilizing both non-informative and informative priors, as well as marginal posterior densities for both censored
and uncensored data. The study aims to explore Bayesian estimators (BEs) and their associated posterior risks (PRs) under various loss functions. Additionally, the paper seeks to detail the estimation of unknown parameters of the proposed distribution through the MCMC algorithm for posterior summaries, encompassing different loss functions and prior types. The work will also provide a mathematical and numerical discussion of Bayesian credible intervals (BCIs) and conclude with an analysis of a real-life dataset.
II. Transmuted Weibull Distribution
As introduced by Alzaatreh et al. in 2013 [10], a random variable X follows a transmuted probability distribution if its probability density function (pdf) and cumulative distribution function (CDF) are given by:
f(x)=g(x)(l+y-y.2G(x)) Fix) = Gix)[il+Y)-Y-Gix)]
where x > 0 and the transmutation parameter y satisfies \y\ < 1. Here, G(x) is the CDF of the baseline distribution, and the functions fix) and Fix) represents the pdf and CDF of the transmuted distribution, respectively.
A random variable X is defined to follow a Weibull distribution characterized by parameters a > 0 and p > 0 if its probability density function (PDF) is given by:
a _i / xa\
gix-,a,p) =-pi*)a ^ 0
The corresponding cumulative distribution function (CDF) for this Weibull distribution can be
expressed as:
Gix) = l-exp (-y)
To find the CDF of the transmuted Weibull distribution, we substitute G (x) into the following formula:
Fix) = (1 + y) - exp (- - y - exp (-
Through algebraic manipulation, we derive the CDF for the transmuted Weibull distribution
Fix; G) = exp (- ) - Y + Y exp jj
where 9 = ia,p,y). To determine the PDF of the transmuted Weibull distribution, we differentiate this CDF with respect to x and simplify the result. The resulting PDF is then
expressed as:
= i(- ffi) ■(1-r + 2reXP(-$1)
x > 0,a,p > 0 and \y\ < 1
Jeslin J, Radhika A, Haripriya M
BAYESIAN PARAMETER ESTIMATION FOR CENSORED RT&A, No 1 (82) TRANSMUTED WEIBULL DISTRIBUTION_Volume 20, March 2025
Special Cases
• When y = 0, the Transmuted Weibull distribution simplifies to the standard Weibull distribution.
• Ifa = 1, the result is the transmuted exponential distribution. Furthermore, with y = 0, it becomes the standard exponential distribution.
• Setting both a and |3 to 1 yields the transmuted standard exponential distribution.
• When a = 2, we obtain the transmuted Rayleigh distribution.
• Ify = 0, this corresponds to the traditional Rayleigh distribution.
III. Likelihood Functions for various Sampling Schemes
Consider a complete random sample X1,X2, —Xn of size n taken from the transmuted Weibull distribution. Then, the likelihood function for the complete data set is given by:
j(a - 1) ]T logXi J exp j- £ ^J [] [l - y + 27 exp
L(x; 8) = anexp
where 9 = (a,p,y) and x = x1,x2, — xn. In many life-testing experiments, it's not possible to collect complete failure time data due to time and cost constraints. As a result, censoring plays an essential role in lifetime data analysis. LetX = X1,X2,— Xr be a type-I censored sample of size r from n items, where the lifetimes follow a transmuted Weibull distribution with parameters a,p & y. Consider a distribution with specific parameters. In the context of Type I censoring, it's important to note that the censoring time is predetermined, while the number of observed failures is random. Suppose we have n items under life testing, and we observe r failures at times t1( t2,... tr. Here, r is an integer between 0 and n, and (n — r) represents the number of items that survive or remain uncensored. According to Mendenhall and Hader (1958) [9], the likelihood function for censored data is given by:
« ^/(*,•). [l-FCnr
7 = 1
where T represents the time, r denotes the number of failures (or uncensored observations), and (n — r) are the censored observations. For a transmuted Weibull distribution applied to censored data, the likelihood function can be expressed as:
(x; 0) « aTexp | j' j x exV ^ 7 + exp
exp
7 = 1
(n — r)log — exp -—j x|l— y + yexp ——
7 = 1
Next, we examine the posterior distribution using Bayes' theorem. The posterior distribution, g(a|x) is given by:
L(x; a)n(a)
g{a \x) =■
fa L(x; a)n(a)da
where n(a) denotes the joint prior distribution of the parameters a = a1,a2,... ak,L(x; a)
Jeslin J, Radhika A, Haripriya M
BAYESIAN PARAMETER ESTIMATION FOR CENSORED RT&A, No 1 (82) TRANSMUTED WEIBULL DISTRIBUTION_Volume 20, March 2025
represents the likelihood function, and g(a I x) is the joint posterior distribution.
IV. Posterior Distribution using Uniform Prior (UP)
In Bayesian estimation, to determine unknown parameters, we specify a prior for each parameter that isn't explicitly defined by the model itself. Unlike the frequentist approach, the Bayesian method incorporates both prior knowledge about the parameters and the observed data. When prior information about the parameters is lacking, a non-informative prior can be used in Bayesian analysis. This type of prior conveys minimal information about the parameters, reflecting a lack of strong prior beliefs.
As introduced by Yousaf et al. in 2020 [11], to estimate the unknown parameters of the transmuted Weibull distribution, we assume the following prior distributions: a « l,Va £ (0, ra),^ k 1,V/? £ (0, «>) and y K l,Vl£ (-1,1). With the assumption that these parameters are independent, the joint prior distribution for a,p and lis: ni{a,p,y) k 1, where a,p > 0 and \y\ < 1. Using Bayes' theorem, the joint posterior distribution of parameters a,p and y , given data x, with a uniform prior is:
g(6\x) = „ „ [1L(X'9ln(9l-— where L(x-,9) is the likelihood function and n(6) represents the
J0 Jo 8)dydpda
uniform prior over the parameters a,p,y.
g(0\x) =
C Jo" J.1! aAoi_1expC-a^11) ^exp ^l) exp(A31)dydpda
where A01 = l + n,A11=Yf=1log-^:/ ¿21=Z?=1*f and A31=2f=1 log {l - y + 2y exp (-^f-)}. Likewise, for censored data, the posterior distribution is given by:
aBoi"1exp(-aB11)-!rexp(-^|l)exp(B31)
g{0\X) ~ jo-jo-j_iiaBo1-iexp(-«B11)7iFexp(-5li)exp(B31)drd^da ()
where B01 = 1 + r, B11=Yj=1 log B21=T1rj=1 x? and B31=Yj=1 log (l - y + 2y exp ( - +
(n - r)log [l - exp (- y) x|l- y + yexp
Since the posterior distributions for both censored and uncensored data are not available in closed form, the marginal posterior densities of the parameters a, p and y for both censored and uncensored data are obtained by integrating out the nuisance parameters, i.e., g(a|x) = jp jxg(a,p,y\x) dpdy and vice versa. Therefore, we use the MCMC technique to obtain the posterior summaries.
V. Posterior Distribution using Informative Prior (IP)
An informative prior offers specific, well-defined information about parameters through a probability distribution. In this study, we assume that the prior distributions of a, p and y are independent. Specifically, we assume Gamma(a,b) for a, JnvGamma(c,d) for p and Uniform(l1, l2) for y . The joint prior distribution of parameters a,p and y is:
g(9) « aa_1e~fca-——e P
l2~li
The joint posterior distribution of the parameters a,p &y and assuming an informative prior (IP) for the complete data, is:
gcii~1exp(-g011) *i_1expi-^|i)exp(031)
= f»f»fi C11 i , n f *-( fn , a'P>0 and \Y\ ^ 1
J0 Jo Li«Cl1 1exp(-g011)^C2i_1exp(--^jexp(031)dydgdiS
1 if Where C11 = a + n,D11 = b + log —, C21 = n + c,D21 = d + log — and D31=Yi=± log |1 —
y + 2yexp For censored data, the joint posterior distribution of a,p and y given the data,
gc" 1exp(-g012)nC71 expi-^|^)exp(032)
9^e\x)= C11 t [ "--4J————aip> 0 and |y|<l (2)
Where C12 = a + r,D12=b+l log ± C22 = c + r,D22 = d + £ xf and D32=E"=i log [l-y +
2yexp + (n — r)log |l — exp (l ~Y + yexp The marginal posterior densities
of the parameters a,p and y for both uncensored and censored data are obtained by integrating out the nuisance parameters, that is, g(a\x) = J^ j^g(a, p,y\x) dpdy and vice versa.
VI. Bayes Estimators and Posterior Risks for different Loss Functions
To estimate an unknown parameter in Bayesian analysis, a loss function must be specified. The choice depends on the specific problem, though there are no strict rules for selecting one. Loss functions can be symmetric (equal weighting to over- and underestimation) or asymmetric. For a decision d, a loss function L{p,d) > 0 represents the incurred loss when estimating unknown parameters a,p and y, and by decisions di,d2 and d2. The expected loss, or posterior risk, denoted by R(d), is given by:
R(d) = Ee]x{L{p,d)} = j L0M)p№) dp
Bayes estimators and their respective posterior risks are computed under the Absolute Error Loss Function (AELF), precautionary loss function (PLF), and quadratic loss function (QLF). Table 1 presents the expressions of Bayes estimators under various loss functions, along with their corresponding posterior risks.
Table 1: BEs and PRs for Various Loss Functions
Loss AELF PLF QLF
Function
Expression L(ß,d) = (ß-d)2 1(M>-V2 w^v'2
BEs d = Eß]x(ß) - ,- - E(P~1\x) d = ä = E(ß_2 w
PRs R(ä) = E{ß2\x)-{E{ß\x)f R(d) - 2 [jE{p2\x) - E(P\x)] R(d) - ^ ^
VII. Posterior Estimates using Markov Chain Monte Carlo
From Equation (2), we observe that the posterior density is in an intractable form, requiring a numerical technique to solve it. Since the posterior summaries are challenging to obtain directly, a Markov Chain Monte Carlo (MCMC) technique is applied, as demonstrated by Carrera B and Papaioannou I, 2024 [12]. To implement MCMC, the posterior densities using both uniform and informative priors are expressed as:
gUP(91*) « fa (n + 1, ^ log i) (n
+ 1. fr (exp ^ log jl - 7 +
glP(ß\x) « fa (n + a,^ log^)fß\a(n + b,xf) /y ie
27 exp (y
+ exp ^^ _ logy, -y + 2y exp
Here, and fa and /misrepresent the probability density functions of the gamma and inverse gamma distributions, respectively, fy while denotes the probability density function of the transmuted parameter. To obtain Bayes estimates, corresponding posterior risks, and Bayesian Credible Intervals (BCI), we proceed as follows: First, a random sample is generated from the transmuted Weibull distribution using the inverse integral transformation, i.e.,
ui = [l-e
1 —y + ye
. After simplification, we obtain:
X; =
-ßlnll
1+7- V(1 + Y)2 ~ 4iijy>
2 y )
where Uj~U(0,l) and i = 1,2, ...,n. By specifying parameter values, a desired random sample can be generated. To produce censored data, a censoring time T is set, and units with values less than or equal to T are recorded. Units with values greater than T are considered censored observations. To implement the MCMC for obtaining posterior summaries, we proceed with the Gibbs sampling steps combined with a Metropolis-Hastings step, as stated by Faucett et al. [13].
VIII. Implementation using Real Life data
The dataset contains remission times for 116 patients diagnosed with acute leukemia. The remission durations (in months) are as follows: 1.08, 0.09, 1.48, 3.87, 13.94, 8.66, 6.11, 23.63, 0.20, 2.23, 3.52, 4.98, 6.97, 9.02, 13.29, 9.40, 2.26, 3.57, 5.06, 7.09, 9.22, 13.80, 25.74, 0.50, 2.46, 3.64, 5.09, 7.26, 9.47, 14.24, 25.82, 0.51, 2.54, 3.70, 5.17, 7.28, 9.74, 14.76, 26.31, 3.81, 0.62, 2.82, 5.32, 7.32, 14.06, 10.77, 32.15, 2.64, 3.88, 5.32, 7.39,10.34,14.83, 34.26, 0.90, 2.69, 4.78, 5.34, 7.59, 0.66,15.96, 36.66,1.05, 2.69, 4.23, 5.41, 7.62, 10.75, 43.01, 1.19, 2.75, 4.26, 5.41, 7.13, 17.12, 46.12, 1.26, 2.83, 4.33, 5.49, 7.66, 11.25, 17.14, 16.62, 79.05, 1.35, 2.87, 5.62, 7.87, 11.64, 17.36, 1.45, 3.02, 4.34, 5.71, 11.93, 7.79, 18.10, 1.46, 4.40, 5.85, 8.26, 11.98, 19.13, 1.76, 3.25, 4.50, 6.25, 8.37, 2.02, 12.02, 3.31, 4.51, 6.54, 8.53 and 22.69.
To estimate the unknown parameters, we applied the methodology from previous sections, utilizing different loss functions and prior distributions. A chi-square test was conducted to verify if the data follow the transmuted Weibull distribution, yielding a p-value of 0.226, indicating a
good fit at the 5% significance level [14]. The Bayesian estimates (BEs), posterior risks (PRs) and Bayesian credible intervals for the parameters a,ß andy of the transmuted Weibull distribution were calculated using uninformative (UP) and informative (IP) priors under AELF, PLF, and QLF loss functions [15]. These results are presented in Tables 2 and 3.
Table 2: BEs and PRs of TWD with hyper parameters a=0.6, b=1.2, c=1.2 & d=1.5
Data condition Loss function Size n UP IP
Complete AELF a = 2.1 1.653 (0.0021) 1.584 (0.0023)
ß = 1.3 1.209 (0.0045) 1.278 (0.0047)
y = 0.7 0.715 (0.0298) 0.703 (0.0312)
PLF a = 2.1 1.664 (0.0052) 1.589 (0.0055)
ß = 1.3 1.229 (0.0038) 1.298 (0.0039)
Y = 0.7 0.733 (0.0755) 0.719 (0.0786)
QLF a = 2.1 1.645 (0.0076) 1.593 (0.0079)
ß = 1.3 1.215 (0.0029) 1.276 (0.0031)
Y = 0.7 0.710 (0.0662) 0.695 (0.0693)
20% Censoring AELF a = 2.1 1.719 (0.0061) 1.623 (0.0064)
ß = 1.3 1.107 (0.0028) 1.216 (0.0030)
Y = 0.7 0.641 (0.0516) 0.624 (0.0549)
PLF a = 2.1 1.732 (0.0075) 1.629 (0.0078)
ß = 1.3 1.139 (0.0042) 1.226 (0.0045)
Y = 0.7 0.629 (0.0843) 0.615 (0.0884)
QLF a = 2.1 1.708 (0.0091) 1.621 (0.0094)
ß = 1.3 1.105 (0.0026) 1.222 (0.0028)
Y = 0.7 0.598 (0.0694) 0.581 (0.0727)
40% Censoring AELF a = 2.1 1.911 (0.0109) 1.823 (0.0113)
ß = 1.3 1.095 (0.0082) 1.200 (0.0085)
Y = 0.7 0.594 (0.0581) 0.573 (0.0618)
PLF a = 2.1 1.928 (0.0121) 1.837 (0.0125)
ß = 1.3 1.098 (0.0065) 1.215 (0.0068)
Y = 0.7 0.621 (0.0821) 0.605 (0.0863)
QLF a = 2.1 1.898 (0.0144) 1.810 (0.0149)
ß = 1.3 1.082 (0.0038) 1.204 (0.0042)
Y = 0.7 0.576(0.0741) 0.556(0.0782)
From Table 2, it is evident that the BEs under both UP and IP have lower posterior risks for uncensored data compared to censored data, due to the information loss associated with censoring. Additionally, the credible intervals for uncensored data were narrower than those for censored
Table 3: 95% Bayesian Credible Intervals of TWD using UP and IP with Hyperparameters parameters a=0.6, b=1.2,
c=1.2 & d=0.5
Data Parameters UP Lower Limit UP Upper Limit IP Lower Limit IP Upper Limit
Complete a = 2.1 0.6234 2.9483 0.6257 2.9461
ß = 1.3 0.1402 1.3871 0.1415 1.3381
Y = 0.7 0.6543 1.1102 0.6401 1.1655
20% Censoring a = 2.1 0.7483 P = 1.3 0.1521 _y = 0.7_0.6798
3.1879 1.4632 1.1387
0.7491 0.1578 0.5964
3.1904 1.3463 1.2476
40% Censoring a = 2.1 0.9145 P = 1.3 0.2283 _Y = 0.7_0.7205
3.5529 1.5874 1.1618
0.9170 0.2305 0.5960
3.5224 1.3812 1.2393
The table 3 presents 95% Bayesian credible intervals for the parameters a = 2.1,p = 1.3 and y = 0.7 of the transmuted Weibull distribution under both uninformative prior (UP) and informative prior (IP) approaches. Increased censoring rates (20% to 40%) lead to wider credible intervals, indicating greater uncertainty in parameter estimates due to loss of information. The IP generally results in narrower intervals compared to the UP, suggesting the benefit of incorporating prior information. Overall, Bayesian credible intervals provide precise estimates for uncensored data, and even with censoring, the intervals remain reasonable. This analysis underscores the effectiveness of Bayesian methods for parameter estimation in incomplete data scenarios, balancing precision and uncertainty.
This article presents a Bayesian analysis of the transmuted Weibull distribution, utilizing both uniform and informative gamma priors under the AELF, PLF, and QLF loss functions. Real-world studies were conducted to evaluate the performance of the Bayes estimators, along with strategies for selecting suitable priors and loss functions across varying sample sizes and test termination times, under both complete and censored data settings. Specifically, two censoring rates—20% and 40%—were examined. Tables 2 and 3 show that the Bayes estimates demonstrated consistency, approaching the true parameter values as sample sizes grew. Posterior risks (PRs) were higher for censored data than for uncensored data, and 95% credible intervals became narrower with larger sample sizes. These findings were consistent in practical applications, supporting the effectiveness of the proposed MCMC algorithm for Bayesian parameter estimation. Future research could expand this work by studying mixtures of transmuted Weibull distributions or applying Bayesian analysis to record values with the transmuted Weibull distribution.
References
[1] Wang, Y., Albalawi, O., Alshanbari, H. M. and Alsubaie, H. H. (2024). A modified cosine-based probability distribution: Its mathematical features with statistical modeling in sports and reliability prospects. Alexandria Engineering Journal, 109:322-333.
[2] Alrweili, H. (2024). Analysis of recent decade rainfall data with new exponential-exponential distribution: Inference and applications. Alexandria Engineering Journal, 95:306-320.
[3] A1 Sobhi, M. M. (2022). The extended Weibull distribution with its properties, estimation and modeling skewed data. Journal of King Saud University - Science, 34:101801.
[4] Abubakar, H., Misiran, M. and Sayed, A. I. A. (2024). Estimation of shifted weibull distribution parameters using optimization algorithms for optimal investment decisions making. Franklin Open, 8:100152.
[5] Alharbi, T. and Hamad, F. (2024). An Optimal Strategy for Estimating Weibull distribution Parameters by Using Maximum Likelihood Method. Statistics, Optimization & Information Computing, 12:1342-1351.
[6] Awopeju Kabiru Abidemi and Abiodun, A. A. (2023). Transmuted Modified Weibull Distribution for Modeling Skewed Lifetime Dataset; Properties and Application. Journal of
IX. Discussion
Probability and Statistical Science, 21.
[7] Omar, A., Nurettin, S., Mohamed, O., and Ibrahim, G. (2023) State-of-the-art review of occupant behavior modeling and implementation in building performance simulation. Renewable and Sustainable Energy Reviews, 185:113558
[8] Yousaf, R., Ali, S. and Aslam, M. (2021). On the Bayesian analysis of two-component mixture of transmuted Weibull distribution. Scientia Iranica, 28:1711-1735.
[9] Mendenhall, W. and Hader, R. J. (1958). Estimation of Parameters of Mixed Exponentially Distributed Failure Time Distributions from Censored Life Test Data. Biometrika, 45:504-520.
[10] Alizadeh, M., Merovci, F., and Hamedani, G. (2016). Generalized transmuted family of distributions: Properties and applications. Hacettepe University Bulletin of Natural Sciences and Engineering Seiries B: Mathematics and Statistics, 46:10.
[11] Yousaf, R., Ali, S., and Aslam, M. (2020). Bayesian estimation of transmuted Weibull distribution under different loss functions. Journal of Reliability and Statistical Studies, 13:1-40.
[12] Carrera, B. and Papaioannou, I. (2024). Covariance-based MCMC for high-dimensional Bayesian updating with Sequential Monte Carlo. Probabilistic Engineering Mechanics, 77:103667.
[13] Faucett, C., and Thomas, D. (1996). Simultaneously modelling censored survival data and repeatedly measured covariates: a Gibbs sampling approach. Statistics in medicine, 15 15:1663-85.
[14] Jian, G. U. O., Zhaojun, L. I. and Thomas, K. (2018) A Bayesian approach for integrating multilevel priors and data for aerospace system reliability assessment. Chinese Journal of Aeronautics, 31(l):41-53.
[15] Singh, R. (2021). On Bayesian estimation of loss and risk functions. Science Journal of Applied Mathematics and Statistics, 9(3):73-77.