Discussion specifying prior distributions in reliability applications

IF 1.3 4区 数学 Q3 MATHEMATICS, INTERDISCIPLINARY APPLICATIONS Applied Stochastic Models in Business and Industry Pub Date : 2023-09-03 DOI:10.1002/asmb.2812
Alfonso Suárez-Llorens
{"title":"Discussion specifying prior distributions in reliability applications","authors":"Alfonso Suárez-Llorens","doi":"10.1002/asmb.2812","DOIUrl":null,"url":null,"abstract":"<p>Firstly, I want to congratulate the authors in Reference <span>1</span> for their practical contextualization in describing the Bayesian method in real-world problems with reliability data. Undoubtedly, one of the main strengths of this article is its highly practical approach, starting from real situations and examples, and showing why Bayesian inference is many times a nice alternative for making estimations. The authors nicely describe how, in reliability applications, there are generally few failure records and, therefore, little information available. For example, this is often the case in the study of the reliability of engineering systems in the army, such as some types of weapons. Since the specific prior is a key aspect of the Bayesian framework, they are primarily concerned with guiding readers on how to make this choice properly.</p><p>Once the parameter of interest <math>\n <semantics>\n <mrow>\n <mi>θ</mi>\n <mo>=</mo>\n <mo>(</mo>\n <mi>μ</mi>\n <mo>,</mo>\n <mi>σ</mi>\n <mo>)</mo>\n </mrow>\n <annotation>$$ \\boldsymbol{\\theta} =\\left(\\mu, \\sigma \\right) $$</annotation>\n </semantics></math> has been identified, and without losing sight of real-world applications, the authors develop their exposition based on three essential premises. Firstly, they remind us that the distribution of <math>\n <semantics>\n <mrow>\n <mi>θ</mi>\n </mrow>\n <annotation>$$ \\boldsymbol{\\theta} $$</annotation>\n </semantics></math> may not always be the main focus of our interest in practical situations. Instead, our key objective might involve estimating cumulative failure probabilities at a specific time or a failure-time distribution <math>\n <semantics>\n <mrow>\n <mi>p</mi>\n </mrow>\n <annotation>$$ p $$</annotation>\n </semantics></math>-quantile, given by the expression <math>\n <semantics>\n <mrow>\n <msub>\n <mrow>\n <mi>t</mi>\n </mrow>\n <mrow>\n <mi>p</mi>\n </mrow>\n </msub>\n <mo>=</mo>\n <mi>exp</mi>\n <mo>[</mo>\n <mi>μ</mi>\n <mo>+</mo>\n <msup>\n <mrow>\n <mi>Φ</mi>\n </mrow>\n <mrow>\n <mo>−</mo>\n <mn>1</mn>\n </mrow>\n </msup>\n <mo>(</mo>\n <mi>p</mi>\n <mo>)</mo>\n <mi>σ</mi>\n <mo>]</mo>\n </mrow>\n <annotation>$$ {t}_p=\\exp \\left[\\mu +{\\Phi}^{-1}(p)\\sigma \\right] $$</annotation>\n </semantics></math>, where <math>\n <semantics>\n <mrow>\n <mi>p</mi>\n <mo>∈</mo>\n <mo>(</mo>\n <mn>0</mn>\n <mo>,</mo>\n <mn>1</mn>\n <mo>)</mo>\n </mrow>\n <annotation>$$ p\\in \\left(0,1\\right) $$</annotation>\n </semantics></math>. Secondly, censored data underly the essence of reliability analysis. Therefore, right, interval, and left censored observations play a fundamental role in all our estimations. Lastly, the authors emphasize that certain reparameterizations of the parameter <math>\n <semantics>\n <mrow>\n <mi>θ</mi>\n </mrow>\n <annotation>$$ \\boldsymbol{\\theta} $$</annotation>\n </semantics></math> can sometimes facilitate the practical interpretation of new parameters and enable greater mathematical tractability. For instance, replacing the usual scale parameter <math>\n <semantics>\n <mrow>\n <mi>exp</mi>\n <mo>(</mo>\n <mi>μ</mi>\n <mo>)</mo>\n </mrow>\n <annotation>$$ \\exp \\left(\\mu \\right) $$</annotation>\n </semantics></math> with a specific quantile <math>\n <semantics>\n <mrow>\n <msub>\n <mrow>\n <mi>t</mi>\n </mrow>\n <mrow>\n <mi>p</mi>\n </mrow>\n </msub>\n </mrow>\n <annotation>$$ {t}_p $$</annotation>\n </semantics></math> can be useful in practice. Building upon these three considerations, the authors exhaustively describe the most commonly used techniques for eliciting a prior distribution and provide a substantial number of bibliographic citations. This fact is valuable in itself because it enables the reader to be aware of the real issues associated with their data and the various approaches available to address the estimation problem.</p><p>One of the most positive aspects of this work is the effort made by the authors to describe most of the known procedures for choosing the prior distribution for the log-location-scale family. The authors provide a summary of the state of the art concerning the elicitation of non-informative distributions, informative distributions, expert opinions, or a combination of various techniques. Specifically, several methods for choosing non-informative priors are comprehensively presented, such as a Jeffreys prior, which is proportional to the square root of the determinant of the Fisher information matrix (FIM), an independence Jeffreys (IJ) prior based on the Conditional Jeffreys (CJ) prior for each parameter, a reference prior that maximizes the Kullback–Leibler divergence between the prior and the expected posterior distribution, and an ordered reference prior that specifies the order of importance of the parameter. The authors also describe the relationships among these non-informative priors, showing when some have advantages over others depending on the nature of the data. In this regard, one of the most useful contributions made by the authors is the elaboration of Table 1 in their article, which summarizes the Jeffreys, IJ, and reference non-informative prior distributions for the log-location-scale family using different parameterizations and censoring scenarios. It is also remarkable Table 2 in their article, where the authors provide a summary of the recommended prior distributions for use with log-location-scale distributions.</p><p>In my opinion, a detailed discussion about the choice of the prior distribution in other complex models could improve the article. For example, in failure processes of heterogeneous repairable systems, which are often modeled by non-homogeneous Poisson processes (NHPP). In these processes, we observe the total number of failures that occur in an engineering system in an interval time <math>\n <semantics>\n <mrow>\n <mo>(</mo>\n <mn>0</mn>\n <mo>,</mo>\n <mi>t</mi>\n <mo>]</mo>\n </mrow>\n <annotation>$$ \\left(0,t\\right] $$</annotation>\n </semantics></math>, and the estimation of the parameters of the intensity function of the model, <math>\n <semantics>\n <mrow>\n <mi>λ</mi>\n <mo>(</mo>\n <mi>t</mi>\n <mo>)</mo>\n </mrow>\n <annotation>$$ \\lambda (t) $$</annotation>\n </semantics></math>, is crucial. It that sense, the authors explicitly mention the work described in Reference <span>2</span>, where it is argued for the use of prior information for <math>\n <semantics>\n <mrow>\n <mi>λ</mi>\n <mo>(</mo>\n <mi>t</mi>\n <mo>)</mo>\n </mrow>\n <annotation>$$ \\lambda (t) $$</annotation>\n </semantics></math>. Other interesting articles in this regard are References <span>3</span> and <span>4</span>, where the authors postulate different forms of the intensity function in a study in Bayesian reliability analysis concerning train door failures on a European underground system. The problem in this context is twofold. First, choosing the intensity function, and secondly, how to evaluate the prior information about its parameters. For example, if the Weibull distribution governs the first system failure, it leads us to the popular power law process (PLP) with intensity function <math>\n <semantics>\n <mrow>\n <mi>λ</mi>\n <mo>(</mo>\n <mi>t</mi>\n <mo>|</mo>\n <mi>θ</mi>\n <mo>)</mo>\n <mo>=</mo>\n <mi>M</mi>\n <mi>β</mi>\n <msup>\n <mrow>\n <mi>t</mi>\n </mrow>\n <mrow>\n <mi>β</mi>\n <mo>−</mo>\n <mn>1</mn>\n </mrow>\n </msup>\n </mrow>\n <annotation>$$ \\lambda \\left(t|\\boldsymbol{\\theta} \\right)= M\\beta {t}^{\\beta -1} $$</annotation>\n </semantics></math>, <math>\n <semantics>\n <mrow>\n <mi>θ</mi>\n <mo>=</mo>\n <mo>(</mo>\n <mi>M</mi>\n <mo>,</mo>\n <mi>β</mi>\n <mo>)</mo>\n <mo>∈</mo>\n <msup>\n <mrow>\n <mi>ℝ</mi>\n </mrow>\n <mrow>\n <mo>+</mo>\n </mrow>\n </msup>\n <mo>×</mo>\n <msup>\n <mrow>\n <mi>ℝ</mi>\n </mrow>\n <mrow>\n <mo>+</mo>\n </mrow>\n </msup>\n </mrow>\n <annotation>$$ \\boldsymbol{\\theta} =\\left(M,\\beta \\right)\\in {\\mathbb{R}}^{+}\\times {\\mathbb{R}}^{+} $$</annotation>\n </semantics></math>. The choice of prior distributions for <math>\n <semantics>\n <mrow>\n <mi>M</mi>\n </mrow>\n <annotation>$$ M $$</annotation>\n </semantics></math> and <math>\n <semantics>\n <mrow>\n <mi>β</mi>\n </mrow>\n <annotation>$$ \\beta $$</annotation>\n </semantics></math> and its connection with the results presented in Tables 1 and 2, in their article, deserve an in-detph study.</p><p>Regarding other interesting models, I would also like to point out about potential applications in Metrology. The role of Metrology in engineering is essential because it ensures the functionality of measuring equipment, proper calibration, and quality control. Metrology suffers from the same problem as reliability data. Due to economic constraints there is a limitation of information where it is common having a random sample of size one, see<span><sup>5</sup></span> for a recent application of Bayesian techniques to evaluate measurement data.</p><p>With respect to the sensitivity, the authors are aware that the choice of a prior distribution could have a strong influence on inferences when having limited information in the data. Therefore, they address in their Section 10 a parametric analysis of sensitivity. At this respect, I would like to pay attention to recent articles about the robustness of Bayesian methods which have potential applications in reliability. All these new methods can improve the robustness study.</p><p>To conclude my discussion, I do think the authors have nicely summarized most of the known methods to choice a specific prior distribution in the log-location-scale family of distributions. As a clear strength of the article, all methods are adjusted to the most practical realities based on censoring scheme. Finally, the exhaustive overview described by the authors opens new problems in other complex model as in the NHPP one. Additionally, the properties of the bivariate likelihood function could lead us to a more precise analysis of sensitivity.</p>","PeriodicalId":55495,"journal":{"name":"Applied Stochastic Models in Business and Industry","volume":null,"pages":null},"PeriodicalIF":1.3000,"publicationDate":"2023-09-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/asmb.2812","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Applied Stochastic Models in Business and Industry","FirstCategoryId":"100","ListUrlMain":"https://onlinelibrary.wiley.com/doi/10.1002/asmb.2812","RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q3","JCRName":"MATHEMATICS, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0

Abstract

Firstly, I want to congratulate the authors in Reference 1 for their practical contextualization in describing the Bayesian method in real-world problems with reliability data. Undoubtedly, one of the main strengths of this article is its highly practical approach, starting from real situations and examples, and showing why Bayesian inference is many times a nice alternative for making estimations. The authors nicely describe how, in reliability applications, there are generally few failure records and, therefore, little information available. For example, this is often the case in the study of the reliability of engineering systems in the army, such as some types of weapons. Since the specific prior is a key aspect of the Bayesian framework, they are primarily concerned with guiding readers on how to make this choice properly.

Once the parameter of interest θ = ( μ , σ ) $$ \boldsymbol{\theta} =\left(\mu, \sigma \right) $$ has been identified, and without losing sight of real-world applications, the authors develop their exposition based on three essential premises. Firstly, they remind us that the distribution of θ $$ \boldsymbol{\theta} $$ may not always be the main focus of our interest in practical situations. Instead, our key objective might involve estimating cumulative failure probabilities at a specific time or a failure-time distribution p $$ p $$ -quantile, given by the expression t p = exp [ μ + Φ 1 ( p ) σ ] $$ {t}_p=\exp \left[\mu +{\Phi}^{-1}(p)\sigma \right] $$ , where p ( 0 , 1 ) $$ p\in \left(0,1\right) $$ . Secondly, censored data underly the essence of reliability analysis. Therefore, right, interval, and left censored observations play a fundamental role in all our estimations. Lastly, the authors emphasize that certain reparameterizations of the parameter θ $$ \boldsymbol{\theta} $$ can sometimes facilitate the practical interpretation of new parameters and enable greater mathematical tractability. For instance, replacing the usual scale parameter exp ( μ ) $$ \exp \left(\mu \right) $$ with a specific quantile t p $$ {t}_p $$ can be useful in practice. Building upon these three considerations, the authors exhaustively describe the most commonly used techniques for eliciting a prior distribution and provide a substantial number of bibliographic citations. This fact is valuable in itself because it enables the reader to be aware of the real issues associated with their data and the various approaches available to address the estimation problem.

One of the most positive aspects of this work is the effort made by the authors to describe most of the known procedures for choosing the prior distribution for the log-location-scale family. The authors provide a summary of the state of the art concerning the elicitation of non-informative distributions, informative distributions, expert opinions, or a combination of various techniques. Specifically, several methods for choosing non-informative priors are comprehensively presented, such as a Jeffreys prior, which is proportional to the square root of the determinant of the Fisher information matrix (FIM), an independence Jeffreys (IJ) prior based on the Conditional Jeffreys (CJ) prior for each parameter, a reference prior that maximizes the Kullback–Leibler divergence between the prior and the expected posterior distribution, and an ordered reference prior that specifies the order of importance of the parameter. The authors also describe the relationships among these non-informative priors, showing when some have advantages over others depending on the nature of the data. In this regard, one of the most useful contributions made by the authors is the elaboration of Table 1 in their article, which summarizes the Jeffreys, IJ, and reference non-informative prior distributions for the log-location-scale family using different parameterizations and censoring scenarios. It is also remarkable Table 2 in their article, where the authors provide a summary of the recommended prior distributions for use with log-location-scale distributions.

In my opinion, a detailed discussion about the choice of the prior distribution in other complex models could improve the article. For example, in failure processes of heterogeneous repairable systems, which are often modeled by non-homogeneous Poisson processes (NHPP). In these processes, we observe the total number of failures that occur in an engineering system in an interval time ( 0 , t ] $$ \left(0,t\right] $$ , and the estimation of the parameters of the intensity function of the model, λ ( t ) $$ \lambda (t) $$ , is crucial. It that sense, the authors explicitly mention the work described in Reference 2, where it is argued for the use of prior information for λ ( t ) $$ \lambda (t) $$ . Other interesting articles in this regard are References 3 and 4, where the authors postulate different forms of the intensity function in a study in Bayesian reliability analysis concerning train door failures on a European underground system. The problem in this context is twofold. First, choosing the intensity function, and secondly, how to evaluate the prior information about its parameters. For example, if the Weibull distribution governs the first system failure, it leads us to the popular power law process (PLP) with intensity function λ ( t | θ ) = M β t β 1 $$ \lambda \left(t|\boldsymbol{\theta} \right)= M\beta {t}^{\beta -1} $$ , θ = ( M , β ) + × + $$ \boldsymbol{\theta} =\left(M,\beta \right)\in {\mathbb{R}}^{+}\times {\mathbb{R}}^{+} $$ . The choice of prior distributions for M $$ M $$ and β $$ \beta $$ and its connection with the results presented in Tables 1 and 2, in their article, deserve an in-detph study.

Regarding other interesting models, I would also like to point out about potential applications in Metrology. The role of Metrology in engineering is essential because it ensures the functionality of measuring equipment, proper calibration, and quality control. Metrology suffers from the same problem as reliability data. Due to economic constraints there is a limitation of information where it is common having a random sample of size one, see5 for a recent application of Bayesian techniques to evaluate measurement data.

With respect to the sensitivity, the authors are aware that the choice of a prior distribution could have a strong influence on inferences when having limited information in the data. Therefore, they address in their Section 10 a parametric analysis of sensitivity. At this respect, I would like to pay attention to recent articles about the robustness of Bayesian methods which have potential applications in reliability. All these new methods can improve the robustness study.

To conclude my discussion, I do think the authors have nicely summarized most of the known methods to choice a specific prior distribution in the log-location-scale family of distributions. As a clear strength of the article, all methods are adjusted to the most practical realities based on censoring scheme. Finally, the exhaustive overview described by the authors opens new problems in other complex model as in the NHPP one. Additionally, the properties of the bivariate likelihood function could lead us to a more precise analysis of sensitivity.

查看原文
分享 分享
微信好友 朋友圈 QQ好友 复制链接
本刊更多论文
关于在可靠性应用中指定先验分布的讨论
首先,我要向参考文献 1 的作者表示祝贺,他们结合实际情况,介绍了贝叶斯方法在实际可靠性数据问题中的应用。毫无疑问,这篇文章的主要优势之一是其高度实用的方法,从实际情况和实例出发,说明了为什么贝叶斯推理在很多时候是进行估算的最佳选择。作者很好地描述了在可靠性应用中,故障记录通常很少,因此可用信息也很少。例如,在研究军队工程系统(如某些类型的武器)的可靠性时,经常会遇到这种情况。一旦确定了感兴趣的参数 θ = ( μ , σ ) $$ \boldsymbol{\theta} =\left(\mu, \sigma \right) $$,在不忽视现实应用的前提下,作者基于三个基本前提展开论述。首先,他们提醒我们,在实际情况中,θ $$ \boldsymbol{\theta} $$ 的分布可能并不总是我们关注的重点。相反,我们的主要目标可能是估计特定时间或故障时间分布 p $ $ p $ $ -quantile 的累积故障概率,其表达式为 t p = exp [ μ + Φ - 1 ( p ) σ ]。 $$ {t}_p=exp \left[\mu +{\Phi}^{-1}(p)\sigma \right] $$ ,其中 p∈ ( 0 , 1 ) $$ p\in \left(0,1\right) $$ 。其次,有删减数据是可靠性分析的基础。因此,右侧、区间和左侧删减观测值在我们的所有估计中都起着根本性的作用。最后,作者强调,参数 θ $$ \boldsymbol{\theta} $$ 的某些重参数化有时可以促进对新参数的实际解释,并使数学可操作性更强。例如,用一个特定的量子点 t p $$ {t}_p $$ 替换通常的尺度参数 exp ( μ ) $$ \exp \left(\mu \right) $$ 在实践中可能很有用。基于这三个方面的考虑,作者详尽地描述了最常用的先验分布诱导技术,并提供了大量的文献引用。这一事实本身就很有价值,因为它能让读者意识到与他们的数据相关的实际问题,以及解决估计问题的各种可用方法。这部著作最积极的方面之一,是作者努力描述了大多数已知的对数位置尺度族先验分布选择程序。作者总结了有关诱导非信息分布、信息分布、专家意见或各种技术组合的技术现状。 具体来说,作者全面介绍了几种选择非信息先验的方法,如杰弗里斯先验(与费雪信息矩阵(FIM)行列式的平方根成比例)、基于每个参数的条件杰弗里斯先验(CJ)的独立杰弗里斯先验(IJ)、使先验与预期后验分布之间的库尔贝克-莱伯勒发散最大化的参考先验,以及指定参数重要性顺序的有序参考先验。作者还描述了这些非信息先验之间的关系,说明了根据数据的性质,某些先验比其他先验更有优势。在这方面,作者最有用的贡献之一是在文章中阐述了表 1,该表总结了使用不同参数化和普查情况下对数位置尺度族的 Jeffreys、IJ 和参考非信息先验分布。在我看来,对其他复杂模型中先验分布的选择进行详细讨论可以改进这篇文章。例如,在异构可修复系统的故障过程中,通常采用非均质泊松过程(NHPP)建模。在这些过程中,我们观察工程系统在时间间隔 ( 0 , t ] 内发生故障的总数。 $$ \left(0,t\right] $$ ,模型强度函数参数 λ ( t ) $$ \lambda (t) $$ 的估计至关重要。从这个意义上说,作者明确提到了参考文献 2 中描述的工作,其中主张使用 λ ( t ) $$ \lambda (t) $$ 的先验信息。参考文献 3 和 4 是这方面其他有趣的文章,在这两篇文章中,作者在关于欧洲地下系统列车门故障的贝叶斯可靠性分析研究中假设了不同形式的强度函数。这方面的问题有两个方面。首先是选择强度函数,其次是如何评估其参数的先验信息。例如,如果 Weibull 分布控制着第一次系统故障,那么就会导致我们使用流行的幂律过程(PLP),其强度函数为 λ ( t | θ ) = M β t β - 1 $$ \lambda \left(t|\boldsymbol\{theta} \right)= M\beta {t}^{\beta -1} $$ , θ = ( M , β ) ∈ ℝ + × ℝ + $$ \boldsymbol{\theta} =\left(M,\beta \right)\in {\mathbb{R}}^{+}\times {\mathbb{R}}^{+} $$ 。关于 M $$ M $ 和 β $$ \beta $ 的先验分布的选择及其与他们文章中表 1 和表 2 所示结果的联系,值得深入研究。计量学在工程中的作用至关重要,因为它能确保测量设备的功能、正确校准和质量控制。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
求助全文
约1分钟内获得全文 去求助
来源期刊
CiteScore
2.70
自引率
0.00%
发文量
67
审稿时长
>12 weeks
期刊介绍: ASMBI - Applied Stochastic Models in Business and Industry (formerly Applied Stochastic Models and Data Analysis) was first published in 1985, publishing contributions in the interface between stochastic modelling, data analysis and their applications in business, finance, insurance, management and production. In 2007 ASMBI became the official journal of the International Society for Business and Industrial Statistics (www.isbis.org). The main objective is to publish papers, both technical and practical, presenting new results which solve real-life problems or have great potential in doing so. Mathematical rigour, innovative stochastic modelling and sound applications are the key ingredients of papers to be published, after a very selective review process. The journal is very open to new ideas, like Data Science and Big Data stemming from problems in business and industry or uncertainty quantification in engineering, as well as more traditional ones, like reliability, quality control, design of experiments, managerial processes, supply chains and inventories, insurance, econometrics, financial modelling (provided the papers are related to real problems). The journal is interested also in papers addressing the effects of business and industrial decisions on the environment, healthcare, social life. State-of-the art computational methods are very welcome as well, when combined with sound applications and innovative models.
期刊最新文献
Issue Information Foreword to the Special Issue on Mathematical Methods in Reliability (MMR23) Limiting Behavior of Mixed Coherent Systems With Lévy-Frailty Marshall–Olkin Failure Times Pricing Cyber Insurance: A Geospatial Statistical Approach Regional Shopping Objectives in British Grocery Retail Transactions Using Segmented Topic Models
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
已复制链接
已复制链接
快去分享给好友吧!
我知道了
×
扫码分享
扫码分享
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1