
This paper describes a new class of boundary value fractional-order differential equations of the q-Hilfer and q-Caputo types, with separated boundary conditions. The presented problem is converted to an equivalent integral form, and fixed-point theorems are used to prove the existence and uniqueness of solutions. Moreover, several special cases demonstrate how the proposed problems advance beyond the existing literature. Examples are provided to support the analysis presented.
Citation: Idris Ahmed, Sotiris K. Ntouyas, Jessada Tariboon. Separated boundary value problems via quantum Hilfer and Caputo operators[J]. AIMS Mathematics, 2024, 9(7): 19473-19494. doi: 10.3934/math.2024949
[1] | Refah Alotaibi, Mazen Nassar, Zareen A. Khan, Ahmed Elshahhat . Analysis of Weibull progressively first-failure censored data with beta-binomial removals. AIMS Mathematics, 2024, 9(9): 24109-24142. doi: 10.3934/math.20241172 |
[2] | Tahani A. Abushal . Parametric inference of Akash distribution for Type-Ⅱ censoring with analyzing of relief times of patients. AIMS Mathematics, 2021, 6(10): 10789-10801. doi: 10.3934/math.2021627 |
[3] | Young Eun Jeon, Suk-Bok Kang, Jung-In Seo . Pivotal-based inference for a Pareto distribution under the adaptive progressive Type-II censoring scheme. AIMS Mathematics, 2024, 9(3): 6041-6059. doi: 10.3934/math.2024295 |
[4] | Nora Nader, Dina A. Ramadan, Hanan Haj Ahmad, M. A. El-Damcese, B. S. El-Desouky . Optimizing analgesic pain relief time analysis through Bayesian and non-Bayesian approaches to new right truncated Fréchet-inverted Weibull distribution. AIMS Mathematics, 2023, 8(12): 31217-31245. doi: 10.3934/math.20231598 |
[5] | Rashad M. EL-Sagheer, Mohamed S. Eliwa, Khaled M. Alqahtani, Mahmoud El-Morshedy . Bayesian and non-Bayesian inferential approaches under lower-recorded data with application to model COVID-19 data. AIMS Mathematics, 2022, 7(9): 15965-15981. doi: 10.3934/math.2022873 |
[6] | Hatim Solayman Migdadi, Nesreen M. Al-Olaimat, Maryam Mohiuddin, Omar Meqdadi . Statistical inference for the Power Rayleigh distribution based on adaptive progressive Type-II censored data. AIMS Mathematics, 2023, 8(10): 22553-22576. doi: 10.3934/math.20231149 |
[7] | Ahmed Elshahhat, Refah Alotaibi, Mazen Nassar . Statistical inference of the Birnbaum-Saunders model using adaptive progressively hybrid censored data and its applications. AIMS Mathematics, 2024, 9(5): 11092-11121. doi: 10.3934/math.2024544 |
[8] | A. M. Abd El-Raheem, Ehab M. Almetwally, M. S. Mohamed, E. H. Hafez . Accelerated life tests for modified Kies exponential lifetime distribution: binomial removal, transformers turn insulation application and numerical results. AIMS Mathematics, 2021, 6(5): 5222-5255. doi: 10.3934/math.2021310 |
[9] | Refah Alotaibi, Mazen Nassar, Zareen A. Khan, Ahmed Elshahhat . Statistical analysis of stress–strength in a newly inverted Chen model from adaptive progressive type-Ⅱ censoring and modelling on light-emitting diodes and pump motors. AIMS Mathematics, 2024, 9(12): 34311-34355. doi: 10.3934/math.20241635 |
[10] | Haiping Ren, Jiajie Shi . A generalized Lindley distribution:Properties, estimation and applications under progressively type-Ⅱ censored samples. AIMS Mathematics, 2025, 10(5): 10554-10590. doi: 10.3934/math.2025480 |
This paper describes a new class of boundary value fractional-order differential equations of the q-Hilfer and q-Caputo types, with separated boundary conditions. The presented problem is converted to an equivalent integral form, and fixed-point theorems are used to prove the existence and uniqueness of solutions. Moreover, several special cases demonstrate how the proposed problems advance beyond the existing literature. Examples are provided to support the analysis presented.
The Weibull distribution, originally introduced by Waloddi Weibull [1], has become an essential tool in statistical analysis across a wide range of scientific and engineering fields. This model is frequently applied to characterize failures in various components and phenomena, particularly within reliability and survival analysis. Over the years, researchers have developed multiple Weibull-related distributions that extend beyond the conventional two- and three-parameter forms often discussed in reliability and statistics [2].
Applications of the Weibull model span multiple disciplines. For instance, Vallée et al. [3] used the Weibull model to estimate the strength of adhesively bonded joints, while Lewis and Withers [4] applied it to investigate particle cracking in metal matrix composites. Beyond engineering, the Weibull distribution is valuable in environmental and medical sciences. Albassam et al. [5] demonstrated its utility in scenarios with indeterminate factors, such as the assessment of unpredictable wind speed data. In medical research, the Weibull model has significant applications in survival analysis and reliability of medical treatments, as seen in studies by Ghazal and Radwan [6].
This flexibility has driven the development of many new Weibull-based models. Examples include the generalized Weibull-modified model by Emam and Alomani [7], which improves parameter estimation, and the weighted Weibull distribution by Xavier and Nadarajah [8], which broadens its applications to complex datasets. Lai et al. [9] added a modified Weibull distribution that effectively models bathtub-shaped hazard rates, and Silva et al. [10] developed the beta-modified Weibull distribution, which is especially useful in survival data analysis for its ability to capture varying hazard functions. Finally, Cousineau [11] provided a comprehensive evaluation of approaches to the three-parameter Weibull model, further extending its potential applications.
While standard Wald-type and profile likelihood intervals are widely used for estimating the Weibull distribution, their performance can be limited, especially with small sample sizes or high parameter interdependence. In this study, we derive formulas for Wald-type confidence intervals for the shape parameter, incorporating both standard and modified profile likelihood methods to address these challenges. The primary contribution is the application of a modified profile likelihood approach, designed to improve interval estimation accuracy under restrictive conditions. The effectiveness of the proposed intervals is assessed through Monte Carlo simulations, comparing them with traditional Wald-type and profile likelihood intervals based on coverage probability and interval length.
The Weibull distribution is a continuous probability distribution characterized by two parameters: Shape (β) and scale (α). The probability density function (pdf) of the Weibull distribution is given by:
f(x;α,β)=(βα)(xα)β−1e−(xα)β,x,α,β>0. | (1) |
The distribution is denoted as X∼W(β,α). Given a random sample X1,X2,...,Xn, the log of the joint likelihood function is:
l(α,β)=log((βα)nn∏i=1(xiα)β−1e−n∑i=1(xiα)β)=nlog(β)−nlog(α)+βn∑i=1log(xi)−n∑i=1log(xi)−n(β−1)log(α)−n∑i=1(xiα)β. | (2) |
As shown in Figure 1, the shape parameter influences the distribution's form, affecting its tail and peak behavior. When β < 1, the distribution has a decreasing probability density function as x increases, indicating a high initial failure rate that decreases over time. For β = 1, the Weibull distribution reduces to the exponential distribution, which is suitable for modeling events with a constant failure rate over time. This could represent components with a random chance of failure that does not change with age. When β > 1, the distribution shows an increasing failure rate, which could model situations where the likelihood of failure increases with age or use, often seen in wear-out failure modes [12,13].
The method of moments (MM) is a classical approach for estimating parameters of statistical distributions, including the Weibull distribution. This method equates the first k theoretical moments of a distribution with the first k empirical moments from sample data to solve for the unknown parameters. For the two-parameter Weibull distribution, the MM approach yields two equations [14]:
˜α=(n−1n∑i=1Xβi)1/1ββ. |
where ˆM1 and ˆM2 are the first and second noncentral sample moments, respectively:
lPF(β)=nlog(β)−nlog(˜α(β))+βn∑i=1log(xi)−n∑i=1log(xi)−n(β−1)log(˜α(β))−n∑i=1(xi˜α(β))β. |
and Γ(⋅) denotes the gamma function.
The maximum likelihood estimator (MLE) is a predominant method for estimating the parameters of the two-parameter Weibull distribution. The estimation is derived from two equations [15]:
ˆαMLE=(n−1n∑i=1XˆβMLEi)1/1ˆβMLEˆβMLEandˆβMLE=[(n∑i=1XˆβMLEiln(Xi))(n∑i=1XˆβMLEi)−1−n−1n∑i=1ln(Xi)]−1. | (3) |
Since MLE does not yield a closed-form solution for the Weibull parameters, numerical optimization techniques are typically employed to find these estimates. The elliptical shape of the contour in Figure 2 illustrates the correlation between the shape and scale parameters. Therefore, a Wald confidence interval, which typically relies on a joint likelihood function, may not be suitable for interval construction. This is because the Wald method does not account for the dependency between the two parameters, which can result in misleading confidence levels. Instead, using the profile likelihood approach for the shape parameter is a more reliable way to estimate the interval in this case because it relies on the observed data and marginalizes over the scale parameter, which better takes into account how the shape and scale parameters are connected in the Weibull model [16,17].
Consider a random sample X1,X2,...,Xn from a Weibull distribution with shape β and scale α. If we consider β as a fixed value, the MLE of α based on the profile likelihood is given by
˜α=(n−1n∑i=1Xβi)1/1ββ. |
Plugging in α in (2) by ˜α, the log profile likelihood function becomes:
lPF(β)=nlog(β)−nlog(˜α(β))+βn∑i=1log(xi)−n∑i=1log(xi)−n(β−1)log(˜α(β))−n∑i=1(xi˜α(β))β. |
Simplifying, this expression becomes:
lPF(β)=nlog(β)−nlog(∑ni=1(xi)β)+βn∑i=1log(xi)+n(log(n)−1)−n∑i=1log(xi). | (4) |
The profile likelihood function is denoted as LPF(˜α,β)=exp[lPF(β)]. The function LPF(˜α,β) can indeed be used as a likelihood function in statistical analysis. This approach serves as a method for approximating the likelihood function by profiling out nuisance parameters.
In the Weibull distribution context, parameter orthogonalization is a technique used to mitigate the potential high correlation between the MLEs for the shape and scale parameters. Reparameterizing to orthogonal parameters can lead to asymptotic independence of the MLEs. This reparameterization involves introducing a new parameter in line with the conditions set by the expected Fisher information matrix [18]. Through a differential equation that incorporates Euler's constant, the orthogonality condition can be expressed as:
jαα∂α∂β+jβα=0, | (5) |
where jαα=(β/βαα)2,jβα=(ξ−1)/(ξ−1)αα, and ξ≈0.5772. Solving Equation (5) to derive an orthogonal nuisance parameter yields λ=αexp((1−ξ)/(1−ξ)ββ). Cox and Reid [18] described a modified profile likelihood method for estimating the shape parameter, which is adjusted for the nuisance scale parameter as follows:
lMPF(β)=lPF(β)−12logdet[Jλλ(β,˜λ(β))], | (6) |
where lPF(β) is shown in Eq (4). The Jλλ(β,˜λ(β)) is the observed information matrix, ˜λ(β) is the restricted MLE of α for the specified β [19]. Yang and Xie [19] showed that Eq (6) can be derived as:
lMPF(β)=(n−2)log(β)−nlog(˜α)+βn∑i=1log(xi)−n∑i=1log(xi)−n(β−1)log(˜α)−n∑i=1(xi˜α)β, | (7) |
and the modified profile likelihood of β is denoted as LMPF(˜α,β)=exp[lMPF(β)].
To solve for the maximum modified profile likelihood estimator for the Weibull shape parameter, the function lMPF(β) can be set to zero. This modified profile likelihood approach is straightforward yet effective, as evidenced by the Monte Carlo simulations [20].
In the field of interval estimation, the analysis of censored data has garnered significant attention from scholars. The Wald method, for example, has been applied to derive approximate confidence intervals for distribution parameters when data is subject to Type-2 censoring [21]. Mweleli et al. [22] developed approximate confidence intervals for the two-parameter Weibull distribution, focusing on small Type-2 censored samples by employing the profile likelihood approach. The intervals obtained through this method are contingent on the shape of the profile likelihood function and lack explicit formulas. Heo et al. [23] introduced methods for constructing confidence limits and intervals for the quantiles of Weibull distributions, utilizing techniques such as MM, probability-weighted moments (PWM), and MLE. Silva and Peiris [24] examined the modeling of rainfall percentiles within the Weibull distribution framework, emphasizing the coverage probability of confidence intervals—a crucial factor for accurate inference regarding rainfall patterns. Vander Wiel and Meeker [25] assessed the precision of s-confidence intervals derived from the likelihood ratio, highlighting their superiority over intervals based on asymptotic normal theory, despite increased computational requirements. Mahdi [15] addressed one-sided conditional and unconditional interval estimation for the scale and shape parameters in a two-parameter Weibull model, drawing inferences from pivotal quantities suggested by Bain and Engelhardt, along with the likelihood ratio method and the Birnbaum statistic.
Researchers have expanded these approaches: Niaki et al. [26] developed Bayesian joint confidence intervals for Weibull parameters, accommodating both complete and censored data. Jana and Bera [27] focused on stress–strength reliability in k-out-of-n systems with inverse Weibull-distributed stress and strength components, proposing asymptotic, bootstrap, and HPD credible intervals. Park [28] examined interval-censored Weibull data estimation, and Yang et al. [29] studied interval estimation for the location parameter in three-parameter Weibull models with a known shape parameter, comparing coverage probability and average length across methods through simulations and a real-world example. Somsamai and Srisuradetchai [30] investigated the coverage probability and average length of confidence intervals for the shape parameter in the Weibull distribution when the scale parameter is unknown. They specifically examined the modified profile likelihood (MPF) and standard profile likelihood (PF) methods. The following expressions present the PF and MPF intervals in terms of the normalized profile likelihood for both MPF and PF:
{β|LPF(˜α,β)maxLPF(˜α,β)⩾exp(−12χ21−α,1)} | (8) |
and
{β|LMPF(˜α,β)maxLMPF(˜α,β)⩾exp(−12χ21−α,1)}. | (9) |
Since there are no closed-form solutions for Eqs (8) and (9), numerical methods are required for their computation. The normalized profile likelihood function of a sample size of 50 from W(β=5,α=1) is shown in Figure 3, with the interval bounds determined to be (3.85, 6.02).
Let θ represent the true parameter, with (Li,Ui) as the lower and upper bounds of the confidence interval for the i-th sample. Then, the coverage probability (CP) can be defined as:
CP=1NN∑i=1I(Li⩽θ⩽Ui), |
where I(⋅) is the indicator function, equaling 1 if θ lies within the interval and 0 otherwise, and N represents the total number of simulated samples. This measure indicates the proportion of intervals that successfully capture the true parameter value [29,31].
Average length (AL) refers to the mean width of the confidence intervals across repeated samples in a simulation study [29]. A smaller AL generally indicates more precise intervals, as the range of values within each confidence interval is narrower [32,33]. The AL can be defined as:
AL=1NN∑i=1(Ui−Li). |
Building on methods discussed in the literature, in this section, we introduce closed-form solutions for Wald-type intervals constructed using the profile likelihood function and the modified profile likelihood under the Weibull distribution. The following results provide explicit formulas for confidence intervals for the shape parameter when both the shape and scale parameters are unknown. These closed-form intervals offer a practical, readily applicable solution in contrast to the PF and MPF intervals presented in Eqs (8) and (9) in the literature, which lack explicit solutions.
Theorem 3.1. Consider a random sample X1,X2,...,Xn from a Weibull distribution with shape β and scale α, both of which are unknown. The corresponding maximum likelihood estimators for these parameters have been obtained. The Wald-type interval using the profile likelihood function (WPF) has a closed form as follows:
ˆβPFMLE±z(1+γ)/(1+γ)22ˆβPFMLE√(n+1n)(11.6449n−0.4904−2ˆβPFMLElog(ˆαMLE)), | (10) |
where z(1+γ)/(1+γ)22 is the (1+γ)/(1+γ)22 quantile of the standard normal distribution, corresponding to a confidence level of γ, γ∈(0,1).
Proof. First, consider β a fixed value. By taking into account Eq (3) and the MLE derived from the profile likelihood, we have: ˜α(β)=(∑ni=1xβin)1β. Substituting α in Eq (2) with ˜α(β), we obtain the log profile likelihood function:
lPF(β)=nlog(β)−nlog(˜α(β))+βn∑i=1log(xi)−n∑i=1log(xi)−n(β−1)log(˜α(β))−n∑i=1(xi˜α(β))β. |
Simplifying further:
lPF(β)=nlog(β)−nlog(∑ni=1(xi)β)+βn∑i=1log(xi)+n(log(n)−1)−n∑i=1log(xi). | (11) |
Note that the last two terms of Eq (4) are independent of the parameter β. The maximum profile likelihood estimator of β is defined as βPFMLE=argmaxlPF(β). Next, we find the score function SPF(β):
SPF(β)=∂lPF(β)∂β=∂∂β(nlog(β)−nlog(∑ni=1(xi)β)+βn∑i=1log(xi)+n(log(n)−1)−n∑i=1log(xi))=nβ−n∑ni=1xβilog(xi)∑ni=1xβi+n∑i=1log(xi). | (12) |
The value of βPFMLE can also be determined by setting Eq (12) to zero. However, deriving an analytical solution for β is complicated, as β appears in both the exponent and outside the summation. Therefore, numerical methods must be employed to solve for β.
The observed Fisher information can be derived from the following expression:
IPF(β)=−∂SPF(β)∂β=−nβ2−n(∑ni=1xβi∂∂β(∑ni=1xβilog(xi))−∑ni=1xβilog(xi)∂∂β(∑ni=1xβi))(∑ni=1xβi)2=−nβ2−n∑ni=1xβi(log(xi))2∑ni=1xβi+n(∑ni=1xβilog(xi)∑ni=1xβi)2. | (13) |
Then, the expected Fisher information is derived as follows:
JPF(β)=E[nβ2+n∑ni=1Xβi(log(Xi))2∑ni=1Xβi−n(∑ni=1Xβilog(Xi)∑ni=1Xβi)2]=nβ2+E[n∑ni=1Xβi(log(Xi))2∑ni=1Xβi]−E[n(∑ni=1Xβilog(Xi)∑ni=1Xβi)2]=nβ2+n((nαββ2)(Γ″(2)+2βΓ′(2)log(α)+β2log2(α)))nαβ−n((nα2ββ2)(Γ″(3)+(n−1)(Γ′(2))2+2βlog(α)(Γ′(3)+(n−1)Γ′(2))+(n+1)β2log2(α)))n(n+1)α2β, |
where Γ(n)(X) is the partial derivative of the gamma function, defined as Γ(n)(X)=∞∫0tx−1e−tlogk(t)dt. For ease of calculation, we have evaluated the following terms: Γ′(2)=0.4228,Γ″(2)=0.8237, Γ′(3)=1.8456, and Γ″(3)=2.4929. Thus, the expected Fisher information will be:
JPF(β)=nβ2+(nβ2)(0.8237+2(0.4228)βlog(α)+β2log2(α))−(n(n+1)β2)(2.4929+(n−1)(0.4228)2+2βlog(α)(1.8456+(n−1)0.4228)+(n+1)β2log2(α))=(nβ2){1.8237+0.8456βlog(α)−2.4929(n+1)−0.1788(n−1)(n+1)−3.6912βlog(α)(n+1)−0.8456(n−1)(n+1)βlog(α)}=n(n+1)β2(1.6449n−0.4904−2βlog(α)). |
The inverse Fisher information of the profile likelihood of β is
JPF(β)−1=(n+1)β2n(1.6449n−0.4904−2βlog(α)). | (14) |
Therefore, the theorem is proved.
Theorem 3.2. Consider a random sample X1,X2,...,Xn from a Weibull distribution with unknown shape β and scale α parameters. The Wald-type interval constructed using the Modified Profile Likelihood (WMPF) function has a closed form as follows:
ˆβMPFMLE±z(1+γ)2ˆβMPFMLE√n+11.6449n2−2.4904n−2−2nˆβMPFMLElog(ˆαMLE), | (15) |
where z(1+γ)/(1+γ)22 is the (1+γ)/(1+γ)22 quantile of the standard normal distribution, corresponding to a confidence level of γ, γ∈(0,1).
Proof. From Eq (7), the score function, SMPF(β), can be derived as follows:
SMPF(β)=∂lMPF(β)∂β=n−2β−n∑ni=1xβilog(xi)∑ni=1xβi+n∑i=1log(xi). |
The observed Fisher information will be:
IMPF(β)=∂SMPF(β)∂β=∂2lMPF(β)∂β2=∂∂β(n−2β−n∑ni=1xβilog(xi)∑ni=1xβi+n∑i=1log(xi))=−n−2β2−n∑ni=1xβi(log(xi))2∑ni=1xβi+n(∑ni=1xβilog(xi)∑ni=1xβi)2. |
Then, the expected Fisher information is as follows:
JMPF(β)=E[n−2β2+n∑ni=1xβi(log(Xi))2∑ni=1xβi−n(∑ni=1Xβilog(Xi)∑ni=1Xβi)2]=n−2β2+E[n∑ni=1xβi(log(Xi))2∑ni=1Xβi]−E[n(∑ni=1Xβilog(Xi)∑ni=1Xβi)2]=n−2β2+(nβ2)(Γ″(2)+2βΓ′(2)log(α)+β2log2(α))−(n(n+1)β2)(Γ″(3)+(n−1)(Γ′(2))2+2βlog(α)(Γ′(3)+(n−1)Γ′(2))+(n+1)β2log2(α)), |
where the gamma function's partial derivatives are evaluated as follows: Γ′(2)=0.4228, Γ″(2)=0.8237,Γ′(3)=1.8456,Γ″(3)=2.4929. Thus, the expected Fisher information is:
JMPF(β)=n−2β2+(nβ2)(0.8237+2(0.4228)βlog(α)+β2log2(α))−(n(n+1)β2)(2.4929+(n−1)(0.4228)2+2βlog(α)(1.8456+(n−1)0.4228)+(n+1)β2log2(α))=nβ2(n+1)(−2n+1.6449n−2.4904−2βlog(α)). |
The inverse Fisher information of the modified profile likelihood of β is:
JMPF(β)−1=(n+1)β21.6449n2−2.4904n−2−2nˆβMLElog(ˆαMLE). | (16) |
Therefore, the theorem is completely proved.
In the simulation study, we evaluated five interval estimation methods—W, PF, MPF, and the proposed WPF and WMPF—across a diverse set of scenarios to assess their performance. The samples were generated from a Weibull distribution with varying shape and scale parameters and multiple sample sizes to provide a comprehensive analysis. Specifically, shape parameters of 0.5, 1, 5, and 10 were paired with scale parameters of 0.5, 1, and 5, and sample sizes of 5, 10, 20, 30, 50, 80,100, and 200 were examined. This setup resulted in 96 scenarios for each method, allowing for an extensive comparison. The results are illustrated in Figure 4.
The shape parameter has a noticeable effect on coverage probability (CP) and average length (AL) of interval estimates, though its impact varies by method. When the scale and sample size are fixed, the shape parameter's influence on CP is minimal for the W, PF, and MPF methods. However, the shape parameter has a stronger effect on CP in the WPF and WMPF methods. For example, with a shape of 0.5 and a sample size of 5, the CPs for the PF method are 0.8971, 0.8925, 0.8976, and 0.8920 for shape values of 0.5, 1, 5, and 10, respectively. By contrast, the WPF method under these same conditions yields CPs of 0.9506, 0.9260, 0.7905, and 0.6858, showing a clear trend: as the shape parameter increases, CP decreases for both WPF and WMPF. This effect diminishes as sample sizes grow beyond 30, at which point the influence of the shape parameter on CP becomes negligible. In terms of AL, an increase in the shape parameter generally results in longer intervals across all methods.
The scale parameter, on the other hand, does not significantly affect CP when the shape parameter and sample size are held constant, especially at lower shape values. However, higher scale values tend to increase AL. For instance, with a shape of 0.5 and a sample size of 5 using the WMPF method, AL values are 1.1951, 1.3590, and 1.7530 for scale values of 0.5, 1, and 5, respectively. Once the sample size exceeds 30, the scale parameter's effect on both CP and AL becomes minimal.
Sample size itself has a substantial influence on both CP and AL. As sample size increases, CP generally stabilizes near the nominal level. For instance, as the sample size grows from 5 to 200, the CP for the W method levels out at the nominal confidence level of 0.95. Similarly, AL decreases as sample size increases, leading to narrower intervals, a trend that holds true across all methods.
When comparing interval methods, performance becomes largely similar when the sample size exceeds 30. For smaller sample sizes (30 or less), the proposed WPF and WMPF methods generally provide better coverage than the others. This advantage is consistent across most scenarios, except when shape values are high (5 or 10), scale is low (0.5), and sample size remains below 30. Under these specific conditions, the MPF and W intervals achieve nominal CP.
To demonstrate the applicability of the proposed interval estimation methods, data from Santiago and Smith [34] concerning hospital-acquired urinary tract infections (UTIs) was analyzed, as shown in Figure 5 (left panel). This dataset tracks 54 male patients who acquired UTIs during their hospital stay, providing essential insights into infection frequency over time—a metric critical for monitoring and intervention within healthcare settings.
The dataset includes time intervals (in days) between UTI cases, with an average interval of approximately 0.21 days, or about 5 hours. A goodness-of-fit test, conducted to assess the fit of the data to a Weibull distribution, yielded a high sample correlation of 0.9909 and a p-value of 0.6614, indicating that the Weibull distribution is an appropriate model [35,36].
Maximizing the joint likelihood function in Eq (2) yielded an estimated shape parameter of approximately 1.0401 and a scale parameter of 0.2138. Based on these estimates, confidence intervals for the shape parameter were constructed using both the proposed methods and traditional approaches (W, PF, and MPF). The resulting interval estimates for each method are as follows: W method (0.8238, 1.2564), PF method (0.8380, 1.2631), MPF method (0.8163, 1.2394), WPF method (0.8251, 1.2551), and WMPF method (0.8226, 1.2576). The closest scenario in the simulation study to this data setup used a shape parameter of 1, a scale parameter of 0.5, and a sample size of 50. In this case, the simulation results showed that WMPF had the highest CP as well as AL, which aligns with the real data findings, where the AL for WMPF is 0.435—the largest among the intervals.
We derived explicit formulas for Wald-type intervals using modified and non-modified profile likelihoods, providing a practical approach for interval estimation. Simulations showed that the WMPF method generally performs best, particularly for small sample sizes (under 30). Additionally, as the shape parameter increases, the CP of the proposed intervals slightly decreases, while changes in the scale parameter have minimal effect. For larger sample sizes (50 or more), all interval methods exhibit similar performance, confirming the reliability of the derived formulas across various settings. Furthermore, the real data application to hospital-acquired urinary tract infections demonstrated the practical value of these intervals in healthcare, supporting early infection monitoring and response.
Future research could extend these Weibull interval estimation methods to high-dimensional data contexts, which often involve multiple predictors or response variables. While we focused on a single response variable, high-dimensionality could arise in similar healthcare settings if multiple factors, such as patient demographics, treatment types, or environmental conditions, are considered alongside infection times. Recent advances, like those by Chaipitak and Choopradit [37] in high-dimensional covariance testing, suggest that adaptive approaches—such as dimension reduction or regularization—may help extend Weibull-based interval estimation to complex, multi-variable datasets. Such extensions could broaden the applicability of Weibull methods in fields requiring robust interval estimation across many variables, particularly in healthcare and engineering reliability.
P. Srisuradetchai: Conceptualization, formal analysis, investigation, methodology, validation, software, writing-original draft, writing–review and editing; J. Somsamai: Data curation, formal analysis, methodology, software, visualization; W. Phaphan: Formal analysis, investigation, project administration, resources, supervision, validation, writing – original draft. All authors have read and approved the final version of the manuscript for publication.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
We express our sincere gratitude for the insightful comments and valuable suggestions provided by the reviewers. This research was funded by Thailand Science Research and Innovation Fund (NSRF), and King Mongkut's University of Technology North Bangkok with Contract no. KMUTNB-FF-68-B-35.
The authors declare that they have no conflicts of interest regarding the publication of this paper.
[1] | A. A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and applications of fractional differential equations, Elsevier, 2006. |
[2] | S. G. Samko, Fractional integrals and derivatives, Theory and Applications, 1993. |
[3] | I. Podlubny, Fractional differential equations, 1999. |
[4] | F. Mainardi, Fractional calculus and waves in linear viscoelasticity: An introduction to mathematical models, World Scientific, 2010. https://doi.org/10.1142/9781848163300 |
[5] |
D. S. Oliveira, E. C. de Oliveira, Hilfer-Katugampola fractional derivatives, Comput. Appl. Math., 37 (2018), 3672–3690. https://doi.org/10.1007/s40314-017-0536-8 doi: 10.1007/s40314-017-0536-8
![]() |
[6] |
T. J. Osler, The fractional derivative of a composite function, SIAM J. Math. Anal., 1 (1970), 288–293. https://doi.org/10.1137/0501026 doi: 10.1137/0501026
![]() |
[7] |
R. Almeida, A Caputo fractional derivative of a function with respect to another function, Commun. Nonlinear Sci. Numer. Simul., 44 (2017), 460–481. https://doi.org/10.1016/j.cnsns.2016.09.006 doi: 10.1016/j.cnsns.2016.09.006
![]() |
[8] |
J. V. da C. Sousa, E. C. de Oliveira, On the ψ-Hilfer fractional derivative, Commun. Nonlinear Sci. Numer. Simul., 60 (2018), 72–91. https://doi.org/10.1016/j.cnsns.2018.01.005 doi: 10.1016/j.cnsns.2018.01.005
![]() |
[9] |
Y. Y., Gambo, F. Jarad, D. Baleanu, T. Abdeljawad, On Caputo modification of the Hadamard fractional derivatives, Adv. Differ. Equ., 2014, (2014), 10. https://doi.org/10.1186/1687-1847-2014-10 doi: 10.1186/1687-1847-2014-10
![]() |
[10] |
F. Jarad, T. Abdeljawad, D. Baleanu, Caputo-type modification of the Hadamard fractional derivatives, Adv. Differ. Equ., 2012 (2012), 142. https://doi.org/10.1186/1687-1847-2012-142 doi: 10.1186/1687-1847-2012-142
![]() |
[11] |
F. Jarad, T. Abdeljawad, D. Baleanu, On the generalized fractional derivatives and their Caputo modification, J. Nonlinear Sci. Appl., 10 (2017), 2607–2619. https://doi.org/10.22436/jnsa.010.05.27 doi: 10.22436/jnsa.010.05.27
![]() |
[12] |
K. Balachandran, S. Kiruthika, J. Trujillo, Existence of solutions of nonlinear fractional pantograph equations, Acta Math. Sci., 33 (2013), 712–720. https://doi.org/10.1016/s0252-9602(13)60032-6 doi: 10.1016/s0252-9602(13)60032-6
![]() |
[13] | D. Vivek, K. Kanagarajan, S. Sivasundaram, Dynamics and stability of pantograph equations via Hilfer fractional derivative, Nonlinear Stud., 23 (2016), 685–698. |
[14] | D. Vivek, K. Kanagarajan, S. Sivasundaram, Theory and analysis of nonlinear neutral pantograph equations via Hilfer fractional derivative, Nonlinear Stud., 24 (2017), 699–712. |
[15] |
A. Anguraj, A. Vinodkumar, K. Malar, Existence and stability results for random impulsive fractional pantograph equations, Filomat, 30 (2016), 3839–3854. https://doi.org/10.2298/fil1614839a doi: 10.2298/fil1614839a
![]() |
[16] |
K. Shah, D. Vivek, K. Kanagarajan, Dynamics and stability of ψ-fractional pantograph equations with boundary conditions, Bol. Soc. Paran. Mat., 39 (2021), 43–55. https://doi.org/10.5269/bspm.41154 doi: 10.5269/bspm.41154
![]() |
[17] |
H. M. Srivastava, Fractional-order integral and derivative operators and their applications, Mathematics, 8 (2020), 1016. https://doi.org/10.3390/math8061016 doi: 10.3390/math8061016
![]() |
[18] |
A. Lachouri, M. S. Abdo, A. Ardjouni, S. Etemad, S. Rezapour, A generalized neutral-type inclusion problem in the frame of the generalized Caputo fractional derivatives, Adv. Differ. Equ., 2021 (2021), 404. https://doi.org/10.1186/s13662-021-03559-7 doi: 10.1186/s13662-021-03559-7
![]() |
[19] |
C. Ravichandran, K. Logeswari, F. Jarad, New results on existence in the framework of Atangana-Baleanu derivative for fractional integro-differential equations, Chaos Soliton. Fract., 125 (2019), 194–200. https://doi.org/10.1016/j.chaos.2019.05.014 doi: 10.1016/j.chaos.2019.05.014
![]() |
[20] |
M. S. Abdo, S. K. Panchal, Fractional integro-differential equations involving ψ-Hilfer fractional derivative, Adv. Appl. Math. Mech., 11 (2019), 338–359. https://doi.org/10.4208/aamm.oa-2018-0143 doi: 10.4208/aamm.oa-2018-0143
![]() |
[21] |
K. Karthikeyan, P. Karthikeyan, N. Patanarapeelert, T. Sitthiwirattham, Mild solutions for impulsive integro-differential equations involving Hilfer fractional derivative with almost sectorial operators, Axioms, 10 (2021), 313. https://doi.org/10.3390/axioms10040313 doi: 10.3390/axioms10040313
![]() |
[22] |
G. Wang, A. Ghanmi, S. Horrigue, S. Madian, Existence result and uniqueness for some fractional problem, Mathematics, 7 (2019), 516. https://doi.org/10.3390/math7060516 doi: 10.3390/math7060516
![]() |
[23] |
A. Morsy, C. Anusha, K. S. Nisar, C. Ravichandran, Results on generalized neutral fractional impulsive dynamic equation over time scales using nonlocal initial condition, AIMS Math., 9 (2024), 8292–8310. https://doi.org/10.3934/math.2024403 doi: 10.3934/math.2024403
![]() |
[24] |
K. Zhao, J. Liu, X. Lv, A unified approach to solvability and stability of multipoint bvps for Langevin and Sturm-Liouville equations with CH-fractional derivatives and impulses via coincidence theory, Fractal Fract., 8 (2024), 111. https://doi.org/10.3390/fractalfract8020111 doi: 10.3390/fractalfract8020111
![]() |
[25] |
H. Srivastava, A. El-Sayed, F. Gaafar, A class of nonlinear boundary value problems for an arbitrary fractional-order differential equation with the Riemann-Stieltjes functional integral and infinite-point boundary conditions, Symmetry, 10 (2018), 508. https://doi.org/10.3390/sym10100508 doi: 10.3390/sym10100508
![]() |
[26] |
M. Alam, A. Zada, T. Abdeljawad, Stability analysis of an implicit fractional integro-differential equation via integral boundary conditions, Alex. Eng. J., 87 (2024), 501–514. https://doi.org/10.1016/j.aej.2023.12.055 doi: 10.1016/j.aej.2023.12.055
![]() |
[27] | K. S. Miller, B. Ross, An introduction to the fractional calculus and fractional differential equations, New York: John Wiley, 1993. |
[28] |
B. Ahmad, S. K. Ntouyas, R. P. Agarwal, A. Alsaedi, Existence results for sequential fractional integro-differential equations with nonlocal multi-point and strip conditions, Bound. Value Probl., 2016 (2016), 205. https://doi.org/10.1186/s13661-016-0713-5 doi: 10.1186/s13661-016-0713-5
![]() |
[29] |
B. Ahmad, S. K. Ntouyas, A. Alsaedi, Sequential fractional differential equations and inclusions with semi-periodic and nonlocal integro-multipoint boundary conditions, J. King Saud Uni. Sci., 31 (2019), 184–193. https://doi.org/10.1016/j.jksus.2017.09.020 doi: 10.1016/j.jksus.2017.09.020
![]() |
[30] |
A. Alsaedi, S. K. Ntouyas, R. P. Agarwal, B. Ahmad, On Caputo type sequential fractional differential equations with nonlocal integral boundary conditions, Adv. Differ. Equ., 2015 (2015), 33. https://doi.org/10.1186/s13662-015-0379-9 doi: 10.1186/s13662-015-0379-9
![]() |
[31] |
B. Ahmad, R. Luca, Existence of solutions for sequential fractional integro-differential equations and inclusions with nonlocal boundary conditions, Appl. Math. Comput., 339 (2018), 516–534. https://doi.org/10.1016/j.amc.2018.07.025 doi: 10.1016/j.amc.2018.07.025
![]() |
[32] |
W. Saengthong, E. Thailert, S. K. Ntouyas, Existence and uniqueness of solutions for system of Hilfer-Hadamard sequential fractional differential equations with two point boundary conditions, Adv. Differ. Equ., 2019 (2019), 525. https://doi.org/10.1186/s13662-019-2459-8 doi: 10.1186/s13662-019-2459-8
![]() |
[33] |
S. Sitho, S. K. Ntouyas, A. Samadi, J. Tariboon, Boundary value problems for ψ-Hilfer type sequential fractional differential equations and inclusions with integral multi-point boundary conditions, Mathematics, 9 (2021), 1001. https://doi.org/10.3390/math9091001 doi: 10.3390/math9091001
![]() |
[34] | S. K. Ntouyas, D. Vivek, Existence and uniqueness results for sequential ψ-Hilfer fractional differential equations with multi-point boundary conditions, Acta Math. Univ. Comenianae, 90 (2021), 171–185. |
[35] |
G. Wang, X. Ren, L. Zhang, B. Ahmad, Explicit iteration and unique positive solution for a Caputo-Hadamard fractional turbulent flow model, IEEE Access, 7 (2019), 109833–109839. https://doi.org/10.1109/access.2019.2933865 doi: 10.1109/access.2019.2933865
![]() |
[36] |
P. Borisut, S. Phiangsungnoen, Existence and uniqueness of positive solutions for the fractional differential equation involving the ρ (τ)-Laplacian operator and nonlocal integral condition, Mathematics, 11 (2023), 3525. https://doi.org/10.3390/math11163525 doi: 10.3390/math11163525
![]() |
[37] |
J. Tariboon, A. Cuntavepanit, S. K. Ntouyas, W. Nithiarayaphaks, Separated boundary value problems of sequential Caputo and Hadamard fractional differential equations, J. Funct. Space., 2018 (2018), 6974046. https://doi.org/10.1155/2018/6974046 doi: 10.1155/2018/6974046
![]() |
[38] |
S. Asawasamrit, A. Kijjathanakorn, S. K. Ntouyas, J. Tariboon, Nonlocal boundary value problems for Hilfer fractional differential equations, Bull. Korean Math. Soc., 55 (2018), 1639–1657. https://doi.org/10.4134/BKMS.b170887 doi: 10.4134/BKMS.b170887
![]() |
[39] | F. H. Jackson, On q-definite integrals, Quart. J. Pure Appl. Math., 41 (1910), 193–203. |
[40] |
F. H. Jackson, XI.–On q-functions and a certain difference operator, Earth Env. Sci. T. R. So., 46 (1909), 253–281. https://doi.org/10.1017/S0080456800002751 doi: 10.1017/S0080456800002751
![]() |
[41] | G. Gasper, M. Rahman, Basic hypergeometric series, Cambridge university press, 2009. https://doi.org/10.1017/cbo9780511526251 |
[42] | V. G. Kac, P. Cheung, Quantum calculus, New York: Springer, 2002. https://doi.org/10.1007/978-1-4613-0071-7 |
[43] |
B. Ahmad, S. K. Ntouyas, I. K. Purnaras, Existence results for nonlocal boundary value problems of nonlinear fractional q-difference equations, Adv. Differ. Equ., 2012 (2012), 140. https://doi.org/10.1186/1687-1847-2012-140 doi: 10.1186/1687-1847-2012-140
![]() |
[44] |
W. A. Al-Salam, Some fractional q-integrals and q-derivatives, P. Edinburgh Math. Soc., 15 (1966), 135–140. https://doi.org/10.1017/s0013091500011469 doi: 10.1017/s0013091500011469
![]() |
[45] |
R. P. Agarwal, Certain fractional q-integrals and q-derivatives, Math. Proce. Cambridge, 66 (1969), 365–370. https://doi.org/10.1017/s0305004100045060 doi: 10.1017/s0305004100045060
![]() |
[46] |
S. Salahshour, A. Ahmadian, C. S. Chan, Successive approximation method for Caputo q-fractional IVPs, Commun. Nonlinear Sci. Numer. Simul., 24 (2015), 153–158. https://doi.org/10.1016/j.cnsns.2014.12.014 doi: 10.1016/j.cnsns.2014.12.014
![]() |
[47] |
W. X. Zhou, H. Z. Liu, Existence solutions for boundary value problem of nonlinear fractional q-difference equations, Adv. Differ. Equ., 2013 (2013), 113. https://doi.org/10.1186/1687-1847-2013-113 doi: 10.1186/1687-1847-2013-113
![]() |
[48] |
S. Abbas, M. Benchohra, N. Laledj, Y. Zhou, Existence and Ulam stability for implicit fractional q-difference equations, Adv. Differ. Equ., 2019 (2019), 480. https://doi.org/10.1186/s13662-019-2411-y doi: 10.1186/s13662-019-2411-y
![]() |
[49] |
N. Allouch, J. R. Graef, S. Hamani, Boundary value problem for fractional q-difference equations with integral conditions in banach spaces, Fractal Fract., 6 (2022), 237. https://doi.org/10.3390/fractalfract6050237 doi: 10.3390/fractalfract6050237
![]() |
[50] | P. Rajkovic, S. Marinkovic, M. Stankovic, On q-analogues of Caputo derivative and Mittag-Leffler function, Fract. Calc. Appl. Anal., 10 (2007), 359–373. |
[51] | N. Limpanukorn, Existence of solution to q-Hilfer fractional difference equation with a time-varying order of operations, Mathematical Journal by The Mathematical Association of Thailand Under The Patronage of His Majesty The King, 67 (2022), 1–11. |
[52] | M. A. Krasnosel'ski˘i, Two remarks on the method of successive approximations, Uspekhi Mat. Nauk, 10 (1955), 123–127. |
[53] | A. Granas, J. Dugundji, Fixed point theory, New York: Springer, 2003. https://doi.org/10.1007/978-0-387-21593-8 |
[54] | K. Deimling, Nonlinear functional analysis, Berlin: Springer, 1985. https://doi.org/10.1007/978-3-662-00547-7 |