الطرق الحسابية البيزية لتقدير توزيع q-Weibull Bayesian computational methods for estimation of q-Weibull distribution

شارك:
المجلة: Sigma Journal of Engineering and Natural Sciences – Sigma Mühendislik ve Fen Bilimleri Dergisi
DOI: https://doi.org/10.14744/sigma.2026.2030
تاريخ النشر: 2026-01-01
المؤلف: Ibrahim Sadok
الموضوع الرئيسي: تقدير التوزيع الإحصائي وتطبيقاته

نظرة عامة

توزيع q-Weibull هو توزيع احتمالي متعدد الاستخدامات يُستخدم في مجالات مثل تحليل الموثوقية، وتحليل البقاء، ونمذجة القيم المتطرفة. تستكشف هذه الورقة تقنيات تقدير المعلمات المختلفة لتوزيع q-Weibull، مع التركيز بشكل خاص على خوارزمية متروبوليس-هاستينغز. تقيّم دراسة محاكاة شاملة أداء هذه الطرق عبر أحجام عينات متنوعة، مدعومة بتطبيقات عملية على مجموعات بيانات حقيقية، مما يُظهر مزايا نهج متروبوليس-هاستينغز في سيناريوهات معينة.

تناقش الدراسة الفجوات الكبيرة في تطبيق توزيع q-Weibull من خلال تقييم الطرق الحسابية البايزية لتقدير المعلمات، مقارنةً بتقدير الاحتمالية القصوى (MLE)، ومتروبوليس-هاستينغز المستقل (IMH)، ومتروبوليس المشي العشوائي (RWM). من خلال محاكاة مونت كارلو وتحليل البيانات الواقعية، تكشف الأبحاث أن الطرق البايزية، وخاصة RWM، تتفوق باستمرار على البدائل التقليدية مثل توزيعات فرشيه، وإيشيتا، وبيرنباوم-ساوندرز، وويبول من حيث الانحراف المعياري، ومتوسط الأخطاء المطلقة، وفترات الثقة. على الرغم من أن الطرق البايزية تتطلب حسابات مكثفة وتحتاج إلى وقت وموارد كبيرة، فإن التحسينات في الدقة والموثوقية التي تقدمها تبرر استخدامها، خاصة في مجموعات البيانات الكبيرة. تم تطوير الشيفرة المصدرية للمحاكاة والتطبيقات باستخدام برنامج Matlab.

مقدمة

تناقش مقدمة الورقة القابلية الواسعة لتطبيق توزيع ويبول في مجالات مختلفة، بما في ذلك الطب، والهندسة، والمالية، ومعالجة الصور، بشكل أساسي لنمذجة أوقات البقاء، وتوقع إجهاد المكونات الميكانيكية، وتحليل تقلبات أسعار الأسهم. يتميز توزيع ويبول الكلاسيكي بدالة كثافة الاحتمال (pdf) الخاصة به، والتي تتضمن معلمة شكل $k > 0$ ومعلمة مقياس $\lambda > 0$. تتيح هذه المرونة نمذجة سلوكيات بيانات متنوعة، خاصة في تحليل الموثوقية والبقاء. ومع ذلك، قد لا يتمكن الشكل الكلاسيكي من التقاط البيانات ذات الذيل الثقيل أو السلوكيات المعقدة بشكل كافٍ، مما يحفز تطوير أشكال عامة مثل توزيع q-Weibull، الذي يعتمد على دالة q-exponential.

يعزز توزيع q-Weibull قدرات النمذجة من خلال استيعاب السلوكيات غير الشاملة وخصائص قانون القوة، مما يجعله مناسبًا للتطبيقات في مجالات مثل علوم البيئة وتحليل الأسواق المالية. تقدم الورقة نهجًا بايزيًا لتقدير معلمات توزيع q-Weibull باستخدام تقنيات سلسلة ماركوف مونت كارلو (MCMC)، مع التركيز بشكل خاص على خوارزمية متروبوليس-هاستينغز ونسخها. يتم إجراء تحليل مقارن لتقدير الاحتمالية القصوى (MLE) والطرق البايزية (متروبوليس-هاستينغز المستقل والمشي العشوائي متروبوليس) من خلال محاكاة شاملة، مما يكشف أن نهج المشي العشوائي متروبوليس يتفوق على تنفيذات MCMC القياسية. تؤكد الدراسة على الأهمية العملية لتوزيع q-Weibull في نمذجة مجموعات البيانات المعقدة وتوفر مخططًا منظمًا للأقسام التالية من الورقة، التي ستتناول نموذج q-Weibull، وطرق التقدير، والتطبيقات الواقعية.

طرق

في هذا القسم، يقارن المؤلفون بين أربع طرق تقدير إحصائية: تقدير الاحتمالية القصوى (MLE)، ومتروبوليس-هاستينغز (MH)، ومتروبوليس-هاستينغز المستقل (IMH)، ومشي عشوائي متروبوليس (RWM). يستند المقارنة إلى مقاييس أدائها، والتي من المحتمل أن تشمل معدلات التقارب، والكفاءة الحسابية، ودقة تقديرات المعلمات.

تشير النتائج إلى أنه بينما يوفر MLE تقديرات فعالة في ظل ظروف معينة، فإن طرق سلسلة ماركوف مونت كارلو (MCMC) (MH، IMH، وRWM) تقدم مزايا في السيناريوهات التي تكون فيها دالة الاحتمالية معقدة أو غير قابلة للحل. على وجه الخصوص، تُظهر طريقة RWM متانة في استكشاف فضاء المعلمات، بينما تُظهر IMH كفاءة محسّنة في أخذ العينات من خلال تقليل الارتباط الذاتي في العينات الناتجة. بشكل عام، تؤكد النتائج على أهمية اختيار الطريقة بناءً على الخصائص المحددة للبيانات والنموذج الأساسي.

مناقشة

في مناقشة توزيع q-Weibull، تؤكد الورقة على التحديات المرتبطة بتقدير المعلمات، وخاصة من خلال تقدير الاحتمالية القصوى (MLE). الطبيعة غير الخطية لمشتقات دالة اللوغاريتم الاحتمالي تعقد اشتقاق الحلول التحليلية للمعلمات \( q \)، \( k \)، و\( \lambda \). لمعالجة هذه التحديات، يقترح المؤلفون نهج تحسين مقيد، باستخدام طريقة نيوتن-رافسون لتعزيز الدقة في التطبيقات الهندسية. بالإضافة إلى ذلك، يتم بناء فترات الثقة للمعلمات بناءً على الطبيعية الأسيمتوتية لـ MLEs، مستفيدين من مصفوفة معلومات فيشر لاشتقاق مصفوفة التباين، مما يساعد في تقدير دقة تقديرات المعلمات.

تستكشف الورقة أيضًا الطرق البايزية لتقدير معلمات q-Weibull، مستخدمة توزيعات مسبقة مستقلة لـ \( q \)، \( k \)، و\( \lambda \). يتم تبرير اختيار التوزيعات المسبقة، مثل بيتا لـ \( q \) وغاما لـ \( k \) و\( \lambda \)، بناءً على ملاءمتها لنطاقات المعلمات. يستخدم المؤلفون تقنيات سلسلة ماركوف مونت كارلو (MCMC)، وتحديدًا خوارزمية متروبوليس-هاستينغز، لأخذ عينات من التوزيع البعدي، مما يبرز فعالية الخوارزمية في التنقل عبر فضاءات المعلمات المعقدة. تكشف دراسات المحاكاة أن الطرق البايزية، وخاصة المشي العشوائي متروبوليس (RWM)، تتفوق على MLE من حيث الدقة والموثوقية، خاصة في التطبيقات الواقعية التي تتضمن بيانات عمر التعب. تدعو النتائج إلى قابلية نموذج q-Weibull وميزات الطرق البايزية في التقاط الديناميات المعقدة للبقاء.

القيود

يناقش قسم القيود القيود والتحديات المرتبطة بالطرق المستخدمة في البحث. بينما تُظهر الطرق مزايا أداء ملحوظة، مثل زيادة الدقة أو الكفاءة، إلا أنها تُظهر أيضًا قيودًا معينة قد تؤثر على قابليتها للتطبيق. قد تشمل هذه القيود قضايا تتعلق بالتحجيم، أو التعقيد الحسابي، أو الحساسية لإعدادات المعلمات، مما قد يعيق فعالية الطرق في سيناريوهات العالم الحقيقي المتنوعة.

علاوة على ذلك، يؤكد القسم على الحاجة إلى النظر بعناية في هذه القيود عند تفسير النتائج. يقترح أنه بينما قد تتفوق الطرق في ظل ظروف معينة، قد تتدهور أداؤها في ظروف أخرى، مما يستلزم مزيدًا من البحث لمعالجة هذه النقائص وتعزيز موثوقيتها. بشكل عام، تعتبر القيود الموضحة اعتبارات حاسمة للعمل المستقبلي في هذا المجال.

Journal: Sigma Journal of Engineering and Natural Sciences – Sigma Mühendislik ve Fen Bilimleri Dergisi
DOI: https://doi.org/10.14744/sigma.2026.2030
Publication Date: 2026-01-01
Author(s): Ibrahim Sadok
Primary Topic: Statistical Distribution Estimation and Applications

Overview

The q-Weibull distribution is a versatile probability distribution utilized in fields such as reliability analysis, survival analysis, and extreme value modeling. This paper explores various parameter estimation techniques for the q-Weibull distribution, focusing particularly on the Metropolis-Hastings algorithm. An extensive simulation study assesses the performance of these methods across varying sample sizes, complemented by practical applications to real datasets, demonstrating the advantages of the Metropolis-Hastings approach in specific scenarios.

The study addresses significant gaps in the application of the q-Weibull distribution by evaluating Bayesian computational methods for parameter estimation, comparing Maximum Likelihood Estimation (MLE), Independent Metropolis-Hastings (IMH), and Random Walk Metropolis (RWM). Through Monte Carlo simulations and real-world data analysis, the research reveals that Bayesian methods, especially RWM, consistently outperform traditional alternatives such as Fréchet, Ishita, Birnbaum-Saunders, and Weibull distributions in terms of standard deviation, mean absolute errors, and confidence intervals. Although Bayesian methods are computationally intensive and require substantial time and resources, the improvements in accuracy and robustness they offer justify their use, particularly in large datasets. The source code for the simulations and applications was developed using Matlab software.

Introduction

The introduction of the paper discusses the extensive applicability of the Weibull distribution in various fields, including medicine, engineering, finance, and image processing, primarily for modeling survival times, predicting mechanical component fatigue, and analyzing stock price variations. The classical Weibull distribution is characterized by its probability density function (pdf), which incorporates a shape parameter $k > 0$ and a scale parameter $\lambda > 0$. This flexibility allows it to model diverse data behaviors, particularly in reliability and survival analysis. However, the classical form may not adequately capture data with heavy tails or complex behaviors, prompting the development of generalized forms such as the q-Weibull distribution, which is based on the q-exponential function.

The q-Weibull distribution enhances modeling capabilities by accommodating non-extensive behaviors and power-law characteristics, making it suitable for applications in fields like environmental science and financial market analysis. The paper introduces a Bayesian approach for parameter estimation of the q-Weibull distribution using Markov Chain Monte Carlo (MCMC) techniques, particularly focusing on the Metropolis-Hastings algorithm and its variants. A comparative analysis of maximum likelihood estimation (MLE) and Bayesian methods (Independent Metropolis-Hastings and Random Walk Metropolis) is conducted through extensive simulations, revealing that the Random Walk Metropolis approach outperforms standard MCMC implementations. The study emphasizes the practical significance of the q-Weibull distribution in modeling complex datasets and provides a structured outline for the subsequent sections of the paper, which will delve into the q-Weibull model, estimation methods, and real-world applications.

Methods

In this section, the authors compare four statistical estimation methods: Maximum Likelihood Estimation (MLE), Metropolis-Hastings (MH), Independent Metropolis-Hastings (IMH), and Random Walk Metropolis (RWM). The comparison is grounded in their performance metrics, which likely include convergence rates, computational efficiency, and accuracy of parameter estimates.

The findings indicate that while MLE provides efficient estimates under certain conditions, the Markov Chain Monte Carlo (MCMC) methods (MH, IMH, and RWM) offer advantages in scenarios where the likelihood function is complex or intractable. Specifically, the RWM method demonstrates robustness in exploring the parameter space, while IMH shows improved sampling efficiency by reducing autocorrelation in the generated samples. Overall, the results underscore the importance of method selection based on the specific characteristics of the data and the underlying model.

Discussion

In the discussion of the q-Weibull distribution, the paper emphasizes the challenges associated with parameter estimation, particularly through Maximum Likelihood Estimation (MLE). The non-linear nature of the log-likelihood function’s derivatives complicates the derivation of analytical solutions for the parameters \( q \), \( k \), and \( \lambda \). To address these challenges, the authors propose a constrained optimization approach, utilizing the Newton-Raphson method to enhance accuracy in engineering applications. Additionally, confidence intervals for the parameters are constructed based on the asymptotic normality of MLEs, leveraging the Fisher information matrix to derive the covariance matrix, which aids in estimating the precision of the parameter estimates.

The paper also explores Bayesian methods for estimating the q-Weibull parameters, employing independent prior distributions for \( q \), \( k \), and \( \lambda \). The choice of priors, such as Beta for \( q \) and Gamma for \( k \) and \( \lambda \), is justified based on their suitability for the parameter ranges. The authors utilize Markov Chain Monte Carlo (MCMC) techniques, specifically the Metropolis-Hastings algorithm, to sample from the posterior distribution, highlighting the algorithm’s effectiveness in navigating complex parameter spaces. Simulation studies reveal that Bayesian methods, particularly the Random Walk Metropolis (RWM), outperform MLE in terms of accuracy and robustness, especially in real-world applications involving fatigue life data. The findings advocate for the q-Weibull model’s adaptability and the advantages of Bayesian approaches in capturing complex survival dynamics.

Limitations

The section on limitations discusses the constraints and challenges associated with the methods employed in the research. While the methods demonstrate notable performance advantages, such as increased accuracy or efficiency, they also exhibit specific limitations that may impact their applicability. These limitations could include issues related to scalability, computational complexity, or sensitivity to parameter settings, which may hinder the methods’ effectiveness in diverse real-world scenarios.

Furthermore, the section emphasizes the need for careful consideration of these limitations when interpreting the results. It suggests that while the methods may excel under certain conditions, their performance may degrade under others, necessitating further research to address these shortcomings and enhance their robustness. Overall, the limitations outlined serve as critical considerations for future work in the field.