Bayesian Hierarchical Scale Mixtures of Log-Normal Models for Inference in Reliability with Stochastic Constraint

: This paper develops Bayesian inference in reliability of a class of scale mixtures of log-normal failure time (SMLNFT) models with stochastic (or uncertain) constraint in their reliability measures. The class is comprehensive and includes existing failure time (FT) models (such as log-normal, log-Cauchy, and log-logistic FT models) as well as new models that are robust in terms of heavy-tailed FT observations. Since classical frequency approaches to reliability analysis based on the SMLNFT model with stochastic constraint are intractable, the Bayesian method is pursued utilizing a Markov chain Monte Carlo (MCMC) sampling based approach. This paper introduces a two-stage maximum entropy (MaxEnt) prior, which elicits a priori uncertain constraint and develops Bayesian hierarchical SMLNFT model by using the prior. The paper also proposes an MCMC method for Bayesian inference in the SMLNFT model reliability and calls attention to properties of the MaxEnt prior that are useful for method development. Finally, two data sets are used to illustrate how the proposed methodology works.


Introduction
In some practical situations, reliability analysis must engage positive data.This includes, for example, cases requiring power plant system analysis, software reliability analysis, and hydrology (flows in a river) among others (see [1,2]).Any reliability analysis must be based on precisely defined model in order to estimate the reliability measure of interest from the available data and subsequently provide a logical basis for improving the reliability of a system.The exponential model is the most fundamental parametric model used to establish reliability.We refer to [3] for excellent review of reliability analysis based on exponential model.Complex failure time events such as system failure event, appear to always produce a set of positively skewed observations and show, initially, an increase over time and, then, a decrease in hazard rate.In such cases, alternative to the exponential model, log-normal distribution (denoted by LN (µ, σ 2 )) can be used to characterize and construct a more plausible model for assessing the failure time Y of the event.The density of Y ∼ LN (µ, , y > 0, However, information on the constraint may sometimes be uncertain.In theses cases, we have to analyze reliability of log-normal model with a stochastic constraint on its parameter space, which causes restriction of reliability measures in their functional forms.For example, in a automobile system analysis, we may interest in estimating MTTF (mean time to failure) of steering, brakes, or other mechanisms based on a log-normal model.Suppose it is known beforehand that such an automotive system exhibits MTTF ∈ (a, b), but this information is uncertain.In spite of uncertain prior information, if we adopt a log-normal model whose parameter space is truncated to µ ∈ C (i.e., MTTF ∈ (a, b)) for reliability analysis, we then have two problems.First, we pass over the uncertainty of the prior information associated with our inference.Second, we may lose data information which strongly contradicts the uncertain information µ ∈ C. Furthermore, there are few valid robust models to assess poorly-distributed and fat-tailed failure time observations.A class of robust FT models can be obtained by the use of scale mixtures of log-normal models.The study about the robust model and its reliability estimation have not yet been tackled in the literature.These factors motivate the contents of the present paper.
The remainder of this paper is arranged as follows.Section 2 introduces a class of scale mixtures of log-normal failure time (SMLNFT) models by applying a scale mixture technique to log-normal model.Subsequently, various class reliability measures are obtained.Section 3 provides a two-stage MaxEnt prior of µ by applying Boltzmann's maximum entropy theorem (see, e.g., [14,15]) to the two-stages of prior hierarchy frame according to [16].We also introduce a scale (degree of prior belief) for stochastic constraint presumption µ ∈ C that is accounted for by the MaxEnt prior.Section 4 develops the Bayesian hierarchical SMLNFT model by utilizing the two-stage MaxEnt prior hierarchy involving µ and a scale mixture hierarch of the SMLNFT model.This section also explores Bayesian inference in reliability for the proposed model utilizing an MCMC sampling based approach.This section also develops MCMC sampling method based on the Gibbs sampler and the Metropolis-Hastings algorithm, and discusses Monte Carlo methods in posterior estimation of reliability measures.Section 5 illustrates the empirical performance of the proposed methodology based on a real and artificial data applications involving proposed SMLNFT models.Section 6 provides concluding remarks and discussion.

The Class of SMLNFT Models
Let failure time Y have a log-normal distribution with parameters (µ, κ(η)σ 2 ), denoted by LN (µ, κ(η)σ 2 ), where κ(η) is a suitably chosen positive weight function of a mixing variable η ∼ G(η) and G(η) denotes the cumulative distribution function (cdf) of η.Subsequently, a simple location model for the failure time Y can be formulated in terms of where so that the distribution of X is a scale mixtures of normal distributions, denoted by X ∼ SMN (µ, σ 2 , κ, G).The proposed model for the failure time Y can be defined as follows:

The density of
The SMLNFT model can vary depending on designation for the function κ(η) and the distribution of η.In the special case where the distribution of η degenerates at κ(η) = 1, the SMLNFT model produces a log-normal failure time (LNFT) model.In cases where η ∼ Gamma(ν/2, ν/2) with E[η] = 1 and κ(η) = 1/η are chosen, the SMLNFT model changes to log-t ν failure time (Lt ν FT) or log-Cauchy failure time (LCFT≡ Lt 1 FT) models, allowing for the regulation of model tail distribution by means of the degrees of freedom.The LCFT model has been particulary used for certain survival processes where significant outliers or extreme results may occur (see, e.g., [17]).We also see that the SMLNFT model approximately reduces to log-logistic failure time (LLFT) model (see, [18]), provided that the choices are κ(η) = 4η 2 and η 2 ∼ IG(2.5, 1.233), where IG(α, β) is an inverse gamma distribution with a probability function f (x) = β α x −(α+1) e β/x /Γ(α).This approximation is known to simplify the implementation of MCMC sampling.The LLFT model has been used in survival analysis as a parametric model for mortality rate, hydrology, and networked telerobots (see, e.g., [19]).Following the same procedure as in [20,21], we can construct new robust FT models through differing designations of κ(η) and the distribution of η.Particularly, in the regulation of tail distribution in the FT model, a log-slash failure time (LSFT) model can be obtained by taking κ(η) = η −1 and η ∼ Beta(1/2, 1); that is Beta distribution.Therefore, the proposed class of SMLNFT models defined by Eq.( 1) is flexible enough to include existing FT models (i.e., LNFT, LCFT, and LLFT) and also provide new robust failure time models (such as Lt ν FT and LSFT).
The densities of five SMLN(µ, σ 2 , κ, G) class members are shown in Figure 1 for different designations of the mixing distribution η and κ(η).Figure 1 shows that the tail area of the log-normal density is thinner than that of the other densities.Thus, the use of SMLN(µ, σ 2 , κ, G) models is useful for flexible and robust reliability analysis of failure time data involving fat-tailed empirical distribution.Figure 1 also demonstrates that the class of SMLN(µ, σ 2 , κ, G) models is useful to describe those situation in which early failure or occurrences dominate the distribution.As indicated in Figure 1, the hazard rate of the log-normal distribution initially increases and then decreases, approaching zero in relation to sufficiently large values of times, but this pattern is not common to all the SMLN(µ, σ 2 , κ, G) models.Reliability and hazard rate of the class of SMLNFT models are directly related to those of the SMN (µ, σ 2 , κ, G) distribution.Let Y be failure time or down time formulated by the SMLN (µ, σ 2 , κ, G) model, then respective MTTF (mean time to failure) and variance of Y be

Preprints
In addition, the reliability and hazard rates of the SMLNFT model are given by and where , and φ and Φ are the density and the cdf of a standard normal variable, respectively.Exact expressions for Eq.( 5) and Eq.( 6) based on the above five FT models are available.For example, the expressions based on the Lt ν FT model are where F ν (•) and f ν (•) denote the cdf and density function of a Student-t ν distribution, respectively.The reliable life T γ , denoting the 100(1-γ)th percentile of the failure time distribution SMLN (µ, σ 2 , κ, G) of the SMLNFT model, may also be expressed as The expectations made for the reliability measures (Eq.(4) through Eq.( 7)) are given in respect to the distribution of η.

Stochastic Constraint on Reliability Measures
Let y 1 , . . ., y n be n complete (nontruncated/noncensored) failure times generated from the SMLNFT model in Eq.( 1), and let the reliability measures, defined by Eq.(4) through Eq.( 7), have a stochastic functional constraint in terms of µ and σ 2 , where the constraint does not depend on σ 2 .Thus a reliability measure, written by B(µ, σ 2 ), needs to be located in a restricted interval B(µ, σ 2 ) ∈ (a, b) (or equivalently µ ∈ C) with a degree of prior belief α, where C ⊂ R is either an interval or a set of intervals.Nevertheless, observations from the SMLNFT model often do not provide strong evidence about information of the constraint, µ ∈ C, is true and therefore may appear to contradict the assumption of the model associated with the constraint.
In this situation, a Bayesian approach can be effectively adopted to model the stochastic (or uncertain) functional constraint on parameter µ.Bayesian reliability analysis of the model ( 1) begins with specifying prior distribution π(µ, σ 2 ), which represents information concerning the parameters µ and σ 2 that are combined with the joint probability distribution of y i 's to yield the posterior distribution: where f (y i |µ, σ 2 ) is the density given in Eq.( 3).When there are no constraints on the location parameter µ, then a common joint prior π(µ, σ 2 ) (e.g.Jeffreys prior or normal-inverse gamma prior) can be used, and posterior reliability inference can be performed without any difficulty.When we have sufficient evidence that the constraint condition on the model ( 1) is true (i.e., α = 1), then a suitable restriction on the parameter space such as using a truncated prior distribution, e.g., π(µ, σ 2 )I µ ∈ C , is expected.Here I(•) denotes an indicator function.See, e.g., [22][23][24], for various applications of truncated prior distribution in Bayesian inference.However, it is often the case that prior information about the constraint is not certain to carry out Bayesian reliability analysis.In this case, it is expected that the uncertainty, i.e., Pr µ ∈ C = α, is taken into account in eliciting a prior distribution of µ.

Two-stage MaxEnt prior
Assume that σ 2 is known and that it is possible to specify that partial information concerning the parameter µ ∈ R is of the form but nothing else about prior distribution π(µ).In this case, maximum entropy prior can be obtained by choosing π(θ) that maximizes the entropy in the presence of the partial information in the form of (8).Boltzmann's maximum entropy theorem (see, e.g., [14,15]) tells us that the density π(µ) that maximizes Ent(π), subject to the constraints E[t j (µ)] = t j , j = 1, . . ., k, takes the k-parameter exponential family form where λ 1 , λ 2 , . . ., λ k can be determined, via the k-constraints, in terms of t The theorem can be directly applied to obtain the MaxEnt prior π(µ) of the SMLNFT Model, provided that prior information on the mean and variance of unconstrained µ can be specified, i.e., 0 .Suppose we have additional information that p(µ ∈ C) = α, 0 < α < 1, we then have the following lemma.
Lemma 2. Suppose we can specify partial prior information about µ by E 0 , and µ ∈ C with degree of belief α.Then the MaxEnt prior density of µ is given by where the value of δ (δ ∈ [0, 1]) is determined by the following equation Then by setting λ 1 = 0 and λ 2 = −1/2σ 2 0 , the MaxEnt prior of µ in (9) reduces to Now the second exponential term in the right hand side of ( 12) can be determined by using the stochastic constraint E[I(µ ∈ C)] = α: Among all the possible proper prior densities of the form (12), the choice of ) yields a proper prior density.Further, this choice leads to the Eq.( 11) which satisfies the stochastic constraint From Lemma 2, we see that δ = 0 implies no functional stochastic constraint on µ of the SMLNFT model, while δ = 1 denotes that functional constraint on the model is certain so that Pr(µ ∈ C) = 1.We also see that π MaxEnt (µ) reduces to the normal density π 0 (µ) = φ(µ; θ 0 , σ 2 0 ) for the former case and the truncated normal density π 1 (µ) = φ(µ; θ 0 , σ 2 0 )/ Φ(C; θ 0 , σ 2 0 ) with the support µ ∈ C for the latter case.This indicates that δ ∈ [0, 1] is directly related to the degree α of prior belief on the constraint µ ∈ C, where α = Pr(µ ∈ C).The relationship between them is as follows.
Corollary 3. The degree α of prior belief on the constraint µ ∈ C accounted for by the MaxEnt prior , and its range is Proof.For the cases where δ = 0 and δ = 1, the degree of prior beliefs are ) by the Theorem of [25].
Corollary 3 indicates that π MaxEnt (µ) is useful for eliciting priori uncertainty concerning constraint on µ by varying the value of δ.By applying the work of [26], we can organize π MaxEnt (µ) within the frame of the two-stage prior hierarchy by [16].
Proof.From the hierarchy of the prior distributions, we see that which is equivalent to π MaxEnt (µ).
From now on, we shall call the MaxEnt prior π MaxEnt (µ), expressed by the two-stage prior hierarchy in Theorem 4, as a "two-stage MaxEnt prior" of µ.

Bayesian Hierarchical Model
Let (y 1 , . . ., y n ) be n complete failure times as generated from the SMLNFT model in Eq.(1), and let η i denote contribution of latent variable η in relation to each observation y i .Suppose one of reliability measures defined by Eq.(4) through Eq.( 7) has a stochastic functional constraint in terms of µ, e.g., B(µ) ∈ (a, b) (or equivalently µ ∈ C), with the degree of prior belief α, and suppose C ⊂ R is either an interval or a set of intervals.Then, from Eq.(3) and Eq.( 10), we obtain posterior distribution of µ and σ 2 given by p(µ, where p(σ 2 ) is a prior distribution of σ 2 .This is a complex function for the Bayesian inference in reliability based on the SMLNFT model.
All hyper-parameters, {θ 0 , σ 0 , δ, ν 1 , ν 2 }, are assumed to be given from the prior information of previous studies or alternate sources.In particular, the value of δ is chosen to satisfy Eq.( 11) for given values of α, θ 0 , and σ 2 0 .In cases where the prior information is not available, a convenient strategy of avoiding improper posterior distribution is to use proper priors with fixed hyper-parameters as appropriate quantity to reflect the diffuseness of the priors (i.e., limiting non-informative priors).

The Gibbs Sampler
Let η = (η 1 , . . ., η n ) be the latent variables, and let x = (x 1 , . . ., x n ) be logarithmic failure time data where x i = log(y i ).Based on the Bayesian Hierarchical SMLNFT model structure, the joint posterior distribution of (µ, µ 0 , τ, η) given the observed data x is where f (τ; ν 1 , ν 2 ) is the density of Gamma(ν 1 , ν 2 ) variate and g(η i ) is the density of the mixing variable η i .Note that the joint posterior distribution of the structure parameters, µ and τ, is p(µ, µ 0 , τ, η|x)dµ 0 dη.This is not simplified in an analytic form of known density, and is thus intractable for posterior inference.Accordingly, we treat latent observations in η as hypothetical missing data, and augment the observed data set x with η in posterior analysis.Thus we derive each conditional posterior distribution of µ, µ 0 , τ and η, for posterior inference based on Markov chain Monte Carlo (MCMC).All the full conditional posterior distributions are as follows: (1) The full conditional distribution of µ is an univariate normal given by: where ∆ = (1 − δ)σ 2 0 .(2) The full conditional distribution of µ 0 is a truncated normal given by: where ∆ 0 = δσ 2 0 .(3) The full conditional distribution of τ is a Gamma distribution:

Preprints
(4) The full conditional distributions of η i 's are independent and their densities are where g(η i ) is the density of a mixing variable.
The conditionals in Eq.( 15) through Eq.( 18) define the Gibbs sampler, provided that Eq.( 18) is a known density of standard form.Otherwise, we can define a Metropolis-within-Gibbs sampler that uses M-H (Metropolis-Hastings) step to sample from the conditional in Eq.( 18).See, [27], for the Metropolis-within-Gibbs sampler.

Markov chain Monte Carlo Method
An MCMC method that works with the posterior Eq.( 14) is not complicated since Gibbs sampling of (µ, µ 0 , τ) is routinely implemented based on each of their full conditionals outlined in Subsection 4.2.However, in posterior sampling of η, an M-H sampling algorithm may be used when the conditional posterior density Eq.( 18) does not have explicit form of known distribution.For the MCMC method, one should note the following points.note 1 : With given initial values of (µ, µ 0 , τ, η), implementation of the Gibbs (or Metropolis-within-Gibbs) sampling algorithm consists of drawing repeatedly from distributions Eq.( 15) through Eq.( 18).The R package tmvtnorm and the R package mvtnorm can be used to sample from the conditionals and to calculate δ for a given α from Eq. (11).note 2 : In cases using the LNFT model, η i degenerates at κ(η i ) = 1.This means that the conditional distribution Eq.( 18) can be eliminated from the Gibbs sampler by setting κ(η i ) = 1 for the conditionals of µ and τ. note 3 : When the Lt ν FT model is used for reliability analysis, the last stage of the Bayesian hierarchical model in Eq.( 13) becomes η i ∼ Gamma(ν/2, ν/2) with κ(η i ) = η −1 i .Thus, the conditional distribution in Eq.( 18) yields where η \η i denotes an (n − 1) × 1 vector whose elements are those of η except for η i .Note that LCFT model is a special case of the Lt ν FT model with ν = 1.To allow ν to be determined within the model, one can specify one more prior stage for the Bayesian hierarchy in Eq.( 13).As suggested by [32], a uniform prior on 1/ν (0 < 1/ν < 1) can be considered.To limit model complexity, we consider only fixed ν so that the investigation of different Lt ν FT models is possible.note 4 : The conditional distribution Eq.( 18) for LSFT model is a truncated gamma whose support is η i ∈ (0, 1).note 5 : It is easily seen that the approximate conditional distribution η 2 i for LLFT model is Then, a move to the proposal point η i and postrior sampling of η i 's can be made by using the M-H algorithm suggested by [18].See, e.g., [28], a general procedure for M-H algorithm.note 6 : The convergence of an MCMC algorithm is an important issue for the correct estimation of the posterior distribution of interest.See [29] for an example of multiple convergence diagnosis and output analysis.

Preprints
When the Markov chain is converged, Rao-Blackwellization yields good estimates of µ and τ. note 7 : As measures of model comparison among SMLNFT models, deviance information criterion (DIC) can be used.This measure can be calculated based on extensions of the MCMC method.See [29,30] and references therein for a review and comparisons of such extensions.
When we have exact expression of u(µ, σ 2 ), then the Monte Calro integration, using generated η (k) 's, vanishes from Eq. (19).For example, the estimate of posterior predictive reliability of the Lt ν FT model for specified y is , where (µ (r) , σ 2(r) ), r = 1, . . ., K, denote the rth Gibbs sample obtained from the fitted Lt ν FT model.Now (1 − α) × 100% credible interval of the reliability at time y can be calculated by where R q denotes the qth quantile value of R (r) 's with . When the LNFT model is assumed, Bayesian estimates of reliability measures can be simply obtained by setting κ(η (k) ) = 1 to Eq. (18).That is, the reliability at time y is .

Numerical Illustrations
We provide two data applications to illustrate and demonstrate the performance of the proposed methodology.The first real data application demonstrates the performance of the proposed methodology (MCMC method based on Bayesian hierarchical SMLNFT model) for inference in reliability with stochastic constraint.Proposed SMLNFT models were fitted to the second data set and compared in terms of their data-fitting performances, estimating parameters, and robustness to fat-tailed empirical distribution.For numerical implementations, we developed our program written in R, which is available from the author upon request.

Equipment Failure Data Example
This example considers failure time data (assessed in hours) from a reliability study of assembly line equipment.The data, available in [31], consists of 43 observations of the failure time (Y).It provides summary statistics as listed in Table 1.The Shapiro-Wilks (S-W) and Kolmogorov-Simirnov (K-S) tests were also implemented to determine the log-normality of Y, that is the normality of X = log(Y).The test statistic values and corresponding p-values are also in Table 1.As seen in the table, the formal tests do not reject the normality of X = log(Y), i.e., Y ∼ LN (µ, σ 2 ).Our objective of data analysis is to estimate reliability of the assembly line equipment subject to Pr(µ ∈ C) = α, where C = (c, ∞), c ≥ 7.5 and µ is the mean of log-failure time.We may have priori information about the mean µ of the log-observation model from past studies or an assembly line quality control report.The stochastic constraint on µ indicates the belief that the reliability of the LNFT model decreases slower than the empirical implication, because x = 7.143 < c(c ≥ 7.5).Thus, it is possible to use the Bayesian hierarchical LNFT model can be used to carry out reliability analysis of the assembly line.Following the proposal of this paper, the two-stage MaxEnt prior π Maxent (µ) can be used to elicit the stochastic constraint on µ.
To see the adequacy of the Bayesian hierarchical LNFT model, we fitted the failure time data to five models.The models were Bayesian hierarchical LNFT, LLFT, LCFT, Lt 5 FT, and LSFT models.They were compared in terms of their deviance information criterions (DIC) which a measure of model comparison and adequacy.Based on each model, we ran a MCMC sampling algorithm (using Gibbs sampler) on this data set and generated 100,000 random samples from the conditional posteriors.For the MCMC algorithm, we used a burn-in period of 1000, a thin interval of 100, and the following choice of hyper-parameter values: θ 0 = 7, σ 0 = 1, and ν 1 = ν 2 = 0.001.DIC values of the five models were calculated and they were DIC LNFT = 132.421,DIC LLFT = 135.704,DIC LCFT = 150.209,DIC Lt 5 FT = 135.698,and DIC LSFT = 145.971.Thus, Bayesian hierarchical LNFT model achieves minimum DIC value among the five models.This, along with the formal normality tests in Table 1, provides support in favor of the LNFT model for fitting the data.
Table 2 shows some posterior summary statistics obtained from the MCMC posterior samples based on the LNFT model with various values of C and δ.The statistics include posterior estimate, posterior standard deviation (S.D.), 0.025th quantile, 0.975th quantile, and MC error calculated from the samples.The values in the last four columns of the table were obtained by setting δ = 0.9.In estimating the Mote Carlo (MC) error, we used the batch mean method with 50 batches.See, e.g., [29, pp.39-40] for the batch mean method.The degree α of the prior belief in the constraint (µ ∈ C) were calculated by use of Eq. (11), and their values are also listed in Table 2.The small MC error values listed in Table 2 convince us of the convergence of the MCMC algorithm.The table shows that π MaxEnt (µ) induces a shrinkage effect in the Bayesian estimation of µ with the uncertain interval constraint, i.e., Pr(µ ∈ C) = α.This effect can be seen from the comparison of the posterior estimates of µ obtained from the case of δ = 0 with those from δ ∈ (0, 1].This comparison indicates the followings: (i) The Bayesian estimate of µ, based on the two-stage prior π MaxEnt (µ), shrinks toward the interval C = (c, ∞).The magnitude of shrinkage effect induced by using the proposed prior becomes more evident as value of δ (or the degree α of belief in the interval constraint µ ∈ C) gets larger.Especially, in case of C = (9, ∞) and δ = 0.9, we see that the 95% credible interval for µ does not include μBayes = 7.141 (estimate of µ with no constraint).This highlights the shrinkage effect induced by using π MaxEnt (µ).(ii) When the priori constraint is certain, i.e., δ = α = 1, Table 2 shows that estimate of µ locates in the given constraint C.These advocates the use of Bayesian hierarchical LNFT model for reflecting the uncertain parameter constraint in reliability analysis.Posterior predictive reliability measures of the LNFT model were also estimated by using the failure time data.Figure 2 depicts estimates of posterior predictive reliability and hazard rate functions for the each case where the stochastic constraint is µ ∈ C = (c, ∞), c = 0.5, 0.7, 0.9 for given δ = 0.9.The figure also compares the functions with those obtained by setting c = −∞ which denotes the case where there is no priori constraint on µ, i.e., the case of δ = 0.. Figure 2 clearly indicates that the two-stage MaxEnt prior appropriately reflects the priori stochastic constraint µ ∈ C = (c, ∞) in estimating the posterior predictive predictive measures.The left panel (a) shows that reliability of the constrained LNFT model decreases slower than that of the unconstrained model (with c = −∞), and this phenomenon is more evident as the constrained interval C locates far from the empirical mean value x = 7.143.This in turn affects shape of hazard rate functions in the right panel (b).For c > x, priori constrained LNFT models tend to produce lower hazard rates than the unconstrained LNFT model at all time point, but the pattern of their hazard rate functions is similar.We also see that the hazard rate becomes lower as the value of c gets larger.This coincides with the implication of the reliability functions in the panel (a).

Artificial Data Example
A study using an artificial data is done to evaluate the performance of proposed Bayesian hierarchical estimation methodology based on the class of SMLNFT models.For the study, we have considered Lt 5 FT model with location µ = 2 and scale parameter σ 2 = 0.25.By use of the model, we simulated a sample of size 300 complete failure times.Five models (LNFT, LCFT, Lt 5 FT, LSFT, and LLFT models) were fitted to the same simulated data.Using the data, we ran MCMC sampling algorithm based on each Bayesian hierarchical model with no constraint (i.e.δ = 0) and generated 50,000 random samples from the conditional posteriors by using burn-in period of 5,000 and thin interval of 10.For constructing each hierarchical model, we used the following choice of hyper-parameter values, θ 0 = 0, σ 0 = 3, and ν 1 = ν 2 = 0.001, which reflect the diffuseness of the priors (i.e., limiting non-informative priors) of µ and τ.
Table 3 provides posterior summaries from the posterior samples generated by using the MCMC algorithm developed for each of five models.The small MC error values listed in Table 3 indicate that we have calculated the parameter estimates with precision and the MCMC algorithm has converged to its target distributions as well.DIC values of all models are also presented in the table.We observe that the Lt 5 FT model has the lower DIC value, while the other models results in considerably higher DIC values than the corresponding one under the Lt 5 FT model.This demonstrates model selection function of the methodology.Concerning the posterior estimates, the estimate of σ 2 under the Lt 5 FT model is accurate in comparison to that under the other models.To be more specific, the scale parameter σ 2 ranges from 0.025 to 0.365 with 95% posterior probability for the Lt 5 FT model, while for the other models, none of the credible intervals include true value σ 2 = 0.25 with the same posterior probability.On the other hand, Table 3 shows that the MCMC algorithm applied to each of five models correctly estimates true location parameter value µ = 2.This implies that the proposed Lt 5 FT model is robust to fat tailed failure time observations affecting the scale parameter.

Conclusion
This paper have provided a methodology for Bayesian inference in reliability of SMLNFT models in cases where the interval constraint µ ∈ C must be incorporated due to uncertainty.For failure time modeling with stochastic restriction, we proposed a Bayesian hierarchical model involving a two-stage MaxEnt prior distribution of µ based on Boltzmann's maximum entropy theorem.The two-stage MaxEnt prior reflects uncertainty about the prior constraint, that is, the Bayesian hierarchical SMLNFT models are subject to reliability restriction with uncertainty.Furthermore, we provide an MCMC method for assessing inference in reliability of SMLNFT models based on their Bayesian hierarchical models.The effectiveness of our methodology is demonstrated by conducting two data applications.
We find that the proposed class of SMLNFT models is flexible enough to account for behavior of reliability, unlike the log-normal FT model.We also find a connection between the degree of constraint presumption and the hyper-parameters of the two-stage MaxEnt prior, indicating adequacy of the MaxEnt prior in eliciting the priori uncertain constraint (see, Corollary 3).Finally, we find that the proposed Bayesian hierarchical methodology enables the development of a simple MCMC algorithm to assess the posterior inference in reliability of the SMLNFT models with a stochastic constraint.Our proposed methodology can be easily extended to cases where the failure time data is incomplete, as a result, for example, of truncated or censored data, as well as to other FT models, including exponential, Weibull, and Gamma FT models.In addition, the class of SMLNFT models could accommodate multivariate FT models at higher dimension without much difficulty.The methodological process proposed here in relation to Bayesian estimation can be extended to the multivariate models with a multivariate stochastic constraint.

Definition 1 .
Let the model (1) construct a plausible model for a random failure time Y.Then the model for Y can then be referred to as an SMLNFT (scale mixtures of log-normal failure time) model.The distribution law of Y is a scale mixtures of log-normal (SMLN) distributions with a weight function κ(η) and the cdf G of a mixing variable η.This is written by Y ∼ SMLN (µ, σ 2 , κ, G).

Figure 2 .
Figure 2.Estimated posterior predictive reliability measures for δ = 0.9: (a) Reliability function; (b) Hazard rate function.The value for tick marks in the y-axis of (b) is multiplied by 10 4 .

Table 1 .
Summary statistics for the failure time data.

Table 2 .
Posterior summaries for each model parameters and degree α of prior belief.

Table 3 .
Posterior summaries for parameters of five models.