# Parameter Estimation and Applications of the Weibull Distribution for Strength Data of Glass Fiber

Authors: Yuxuan Wu; Hanyang Xie; Jyun-You Chiang; Gang Peng; Yan Qin
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

Glass fiber is a good substitute for metal materials. However, in the process of manufacturing, it is necessary to carry out sampling inspection on its tensile strength to infer its quality. According to previous literatures, the strength data can be well fitted by the Weibull distribution, while the poor parameter estimation method cannot obtain reliable analysis results. Therefore, a new parameter estimation method is proposed. Based on the simulation results, it is found that the proposed parameter estimation method outperforms the other competitors to obtain reliable estimates of the Weibull parameters. Finally, the proposed parameter estimation method is applied to two real data sets of glass fiber strength for illustration. The results of data analysis show that our proposed parameter estimation method is more suitable for these data sets than other estimation methods.


---

## Body

## 1. Introduction

Glass fiber is a kind of high-performance composite fiber. Compared with carbon fiber, it has lower cost and wider application in the field of construction, electronics, chemical industry, and so on. In the manufacturing process of glass fiber, it is very important to determine the failure probability distribution of several nominally identical glass fiber samples by destructive test and then statistical analysis of their strength data. Figure1 shows the typical fracture pattern for low strength fiber, which is from Feih et al. [1]. According to previous literatures, the strength data of fiber is usually described by two-parameter Weibull distribution. For example, Sakin and Ay [2] studied the bending fatigue properties of glass fiber reinforced polyester composite plates and fitted the data with two-parameter Weibull distribution. Acitas et al. [3] assumed that the strength data of glass fiber followed the Weibull distribution and estimated the parameters with maximum likelihood estimation (MLE). Zhu et al. [4] developed the technology of making high strength refractory ceramic fibers using fly ash, and the mechanical properties of a series of fly ash fibers were assessed by the Weibull distribution. Ferreira et al. [5] studied the Weibull analysis of tensile tested piassava fibers with different diameters.Figure 1
Low strength fracture surface of glass fiber.Because of the good properties of the Weibull distribution, it is currently applied in various fields. Many researches have proposed the related distribution, neutrosophic statistics, or extended distribution and discussed its properties and applications, such as Alkarni et al. [6], Almongy et al. [7], Almetwally et al. [8], Ferretti et al. [9], Zhao et al. [10], Aslam et al. [11], Fernández et al. [12], Arif and Aslam [13], Aslam et al. [14], Aslam [15], Aslam et al. [16], Aslam and Arif [17], and Zhang et al. [18].Before the reliability and failure rate of glass fiber are inferred, sample-based parameters’ estimation is required for the analysis process. There are many methods for estimating the unknown parameters and reliability analysis of the Weibull distribution. Naresh et al. [19] provided reliability analysis of tensile strength of glass/epoxy and carbon/epoxy composites based on the Weibull distribution. Malik et al. [20] evaluated reliability of a non-series-parallel system of seven components, and the Weibull distribution is used for the random variables associated with failure and operating times of the components. Jia [21] studied the reliability for the Weibull distribution with homogeneous heavily censored data. Johnson et al. [22] introduced many commonly used estimation methods for the Weibull distribution in the book, such as MLE, moment estimation (MM), percentile method (PM), method of logarithmic moment (MLM), and Bayesian estimation. Teimouri et al. [23] proposed an L-moment (LM) estimator for the Weibull distribution. Kantar [24] proposed generalized least squares and weighted least squares estimation methods for distributional parameters. Datsiou and Overend [25] investigated four estimation methods of the Weibull distribution for glass strength data. Hidekazu et al. [26] proposed an improved estimation method for MLE, which improved the estimation of Weibull shape parameter, scale parameter, and quantile. Acitas et al. [3] used particle swarm optimization to solve the MLE of the Weibull distribution and applied the results to the strength data of glass fiber. Although the above literatures on parameter estimation and reliability analysis of the Weibull distribution are valuable, there is a lack of research on the parameter estimation method of strength data fitting the Weibull distribution (Datsiou and Overend [25]). The most commonly used estimation method is MLE, while it does not perform well for the Weibull distribution, especially in small samples.The aim of this article is to propose a new parameter estimation method for the Weibull distribution. The performance of proposed method outperforms other common estimation methods for small sample size. In addition, the proposed estimation method is applied to estimate the reliability and failure rate of products to evaluate quality or reliability.The rest of this paper is organized as follows. In Section2, some existing parameter estimation methods for the Weibull distribution are addressed. In Section 3, we propose a new parameter estimation method for the Weibull distribution. Based on the proposed parameter estimation method, it is applied to the estimation of reliability and failure rate. In Section 4, a Monte Carlo simulation is conducted to evaluate the performances of the proposed parameter estimation method. Then, two real strength data sets of glass fiber are used to demonstrate the proposed parameter estimation method and reliability analysis in Section 5, and some concluding remarks are provided in Section 6.


## 2. Review of Parameter Estimation Methods

Suppose the random variableT has a Weibull distribution with shape parameter α>0 and scale parameter β>0, denoted as T∼WDα,β. The probability density function (PDF), cumulative distribution function (CDF), reliability function (RF), and failure rate function (FRF) are, respectively, given by(1)fTt;α,β=αβtβα−1e−t/βα,t>0,(2)FTt;α,β=1−e−t/βα,t>0,(3)Rt=PX>t=e−t/βα,t>0,(4)λt=fXt;α,βRt=αβtβα−1,t>0.We assumeT1,…,Tn is a random sample from the two-parameter Weibull distribution. Five common parameter estimation methods are reviewed as follows.


### 2.1. Method of Moments (MM)

The MM is a traditional parameter estimation method, which is widely used for parameter estimation of various distributions. This estimation method assumes that the firstk population moments are equivalent to the first k sample moments and estimates the unknown parameters by solving the equations derived from the equivalence relation. Therefore, it can be used for two-parameter Weibull distribution, whose parameter estimation satisfies the following two equations:(5)Γ2/α^MM+1Γ21/α^MM+1+M^2−M^12M^1−1=0,(6)β^MM=M^1Γ1/α^MM+1,where α^MM and β^MM are the estimates of α and β, respectively, Γ. is gamma function, and M^1 and M^2 are the first and second noncentral sample moments, respectively.


### 2.2. Method of L-Moment (LM)

According to Teimouri et al. [23], they first proposed a L-moment estimator for two-parameter Weibull distribution. The LM estimator of parameters can be obtained by the following equations:(7)α^LM=−ln2ln1−2/nn−1∑i=1ni−1Ti:n−T¯/T¯,(8)β^LM=X¯Γ1/α^LM+1,where Ti:n denotes the ith order statistic in the sample size n, and T¯ is the average of sample Ti, i = 1, ..., n.


### 2.3. Maximum Likelihood Estimator (MLE)

The MLE method is one of the most commonly used parameter estimation methods. The MLE of two-parameter Weibull distribution can be obtained by the following equations:(9)β^MLE=1n∑i=1ntiαα,(10)1α^MLE=∑i=1ntiαlnti∑i=1ntiα−1n∑i=1nlnti.Since the close-form solution of the MLE is not available, some numerical optimization methods have been proposed to overcome this issue.


### 2.4. Method of Percentile (PM)

According to Johnson et al. [22], the PM estimates of the two-parameter Weibull distribution are(11)β^PM=x^0.632,(12)α^PM=ln−ln1−p2lnt^p−lnt^0.632,where t^p and t^0.632are the estimators of pth and 0.632th percentile, respectively. Former studies show that PM has a good estimation effect when p=0.15or0.31.


## 2.1. Method of Moments (MM)

The MM is a traditional parameter estimation method, which is widely used for parameter estimation of various distributions. This estimation method assumes that the firstk population moments are equivalent to the first k sample moments and estimates the unknown parameters by solving the equations derived from the equivalence relation. Therefore, it can be used for two-parameter Weibull distribution, whose parameter estimation satisfies the following two equations:(5)Γ2/α^MM+1Γ21/α^MM+1+M^2−M^12M^1−1=0,(6)β^MM=M^1Γ1/α^MM+1,where α^MM and β^MM are the estimates of α and β, respectively, Γ. is gamma function, and M^1 and M^2 are the first and second noncentral sample moments, respectively.


## 2.2. Method of L-Moment (LM)

According to Teimouri et al. [23], they first proposed a L-moment estimator for two-parameter Weibull distribution. The LM estimator of parameters can be obtained by the following equations:(7)α^LM=−ln2ln1−2/nn−1∑i=1ni−1Ti:n−T¯/T¯,(8)β^LM=X¯Γ1/α^LM+1,where Ti:n denotes the ith order statistic in the sample size n, and T¯ is the average of sample Ti, i = 1, ..., n.


## 2.3. Maximum Likelihood Estimator (MLE)

The MLE method is one of the most commonly used parameter estimation methods. The MLE of two-parameter Weibull distribution can be obtained by the following equations:(9)β^MLE=1n∑i=1ntiαα,(10)1α^MLE=∑i=1ntiαlnti∑i=1ntiα−1n∑i=1nlnti.Since the close-form solution of the MLE is not available, some numerical optimization methods have been proposed to overcome this issue.


## 2.4. Method of Percentile (PM)

According to Johnson et al. [22], the PM estimates of the two-parameter Weibull distribution are(11)β^PM=x^0.632,(12)α^PM=ln−ln1−p2lnt^p−lnt^0.632,where t^p and t^0.632are the estimators of pth and 0.632th percentile, respectively. Former studies show that PM has a good estimation effect when p=0.15or0.31.


## 3. Proposed Parameter Estimation Method

According to the distribution property, ifY = FTt;α,β, then Y becomes a uniform distribution from 0 to 1. Moreover, the mean and variance of Y are 1/2 and 1/12, respectively, which can be expressed as(13)μY=12,(14)σY2=112.Supposet1,t2,…,tn is a random sample with sample size n from WDα,β. Theoretically, yi=FTti;α,β, i=1,…,n, also follows a uniform distribution from 0 to 1. However, the sample has inevitable random error. Therefore, the idea of the proposed parameter estimation method is to find the estimators of α and β, which can make yi approximate the uniform distribution as much as possible. In other words, the estimators of α and β can satisfy as much as possible that the mean and variance of yi are 1/2 and 1/12, respectively. Because the sample mean y¯=∑i=1nyi/n and sample variance Sy2=∑i=1nyi−y¯2/n−1 are unbiased estimators of population mean and population variance, respectively, the proposed parameter estimators can be obtained by solving the following equations:(15)y¯≅12,Sy2≅112.Equation (15) can be rewritten as follows:(16)∑i=1nexp−tiβα≅n2,∑i=1nexp−tiβα2≅4n−112.Therefore, the estimators ofα and β can be obtained from equation (16) and denoted as α^CDF and β^CDF, respectively. Since the proposed parameter estimation method is based on CDF, we name this method CDFM (cumulative distribution function method). Then, α and β in equations (3) and (4) are replaced by α^CDF and β^CDF, respectively, so the RF and FRF based on CDFM can be obtained and expressed as(17)R^t=e−t/β^CDFα^CDF,t>0,(18)λ^t=α^CDFβ^CDFtβ^CDFα^CDF−1,t>0.Since the close-form solution of the proposed estimation is not available, the estimators can be obtained by using numerical optimization methods such as the quasi-Newton algorithm [27]. Figure 2 shows the flow chart of the parameter estimation process. When we obtain a sample data T, the estimators of the parameters, α^CDF and β^CDF, can be obtained according to the proposed parameter estimation method. Then, the goodness of fit test is used to test whether the data comes from the Weibull distribution with parameters α^CDF and β^CDF. Common goodness of fit tests include Kolmogorov-Smirnov test (K-S test) and Anderson-Darling test. When the test results show that the data does not follow the Weibull distribution, we must find the exact distribution of the data before inference. When the test results show that the data may come from the Weibull distribution, we can use equations (17) and (18) to estimate RF and FRF.Figure 2
Flowchart of proposed parameter estimation.


## 4. Monte Carlo Simulation Study

In this section, all simulation results are obtained by R software. The R codes can be obtained from the authors upon request. In this simulation study, we compare proposed estimation method CDFM with MLE, LM, PM, and MM. For convenience, we denoted the results estimated by these methods asθ^CDFM=α^CDFM,β^CDFM, θ^MLE=α^MLE,β^MLE, θ^LM=α^LM,β^LM, θ^PM=α^PM,β^PM, and θ^MM=α^MM,β^MM, respectively. Considering the sample size n = 10, 20, 30, 40, 50, 60, 70, and data set are generated from the WDα,β, where α=1,2,4, 6 and β=1,2. Ten thousand iteration runs are used to evaluate the bias and MSEs of five estimation methods. The bias and MSE can be evaluated based on equations (19) and (20), respectively:(19)bias=110000∑i=110000η^i−η,(20)MSE=110000∑i=110000η^i−η2,where η can be α or β and the subscript i denotes the five estimation methods. All obtained biases and MSEs are reported in Tables 1 and 2. In view of Tables 1 and 2, the following results are found:(1)
With the increase of sample size, the performance of bias and MSE of all methods are getting better and better.(2)
For the bias of shape parameterα, the performance of α^CDFM outperforms other methods in most cases, followed by α^LM.(3)
For the MSE ofα, the performance of α^LM outperforms other methods in most cases, followed by α^ML and α^CDFM.(4)
For the bias of scale parameterβ, the performance of β^CDFM outperforms other methods for n = 10, except α=1. And the performance of β^LM outperforms other methods in most cases.(5)
For the MSE ofβ, the performance of β^LM outperforms other methods, while the performance of β^MLM is the worst.(6)
With the increase of the population shape parameter, the bias ofα is larger and larger, and the bias of β is smaller and smaller for all methods.(7)
Whenβ = 1, the performance of each estimation method is similar to that of β = 2.Table 1
The bias and MSE of five estimation methods forβ=1,α=1,2,4,6.
β=1α=1α=2α=4α=6nMethodα^iβ^iα^iβ^iα^iβ^iα^iβ^i10CDFMBias0.0530.0420.0970.0040.2400.0010.3130.002MSE0.1290.1230.4950.0292.1480.0074.5640.003MLBias0.1690.0160.3270.0070.7030.0061.0200.006MSE0.1520.1120.5890.0272.5340.0075.5610.003LMBias0.0560.0290.1030.0110.2940.0030.4450.002MSE0.1100.1120.4240.0281.9390.0074.3670.003PMBias0.5070.0160.9810.0151.9880.0122.8420.009MSE1.3440.1534.7620.03820.940.00941.750.004MMBias0.1670.0200.1820.0080.4200.0040.6600.003MSE0.1310.1130.4410.0282.1550.0075.0660.00320CDFMBias0.0290.0210.0450.0020.0940.0010.1370.001MSE0.0510.0610.1920.0150.8160.0041.8510.002MLBias0.0800.0100.1530.0040.3100.0020.4410.003MSE0.0500.0570.1870.0140.7940.0031.7610.002LMBias0.0290.0120.0480.0060.1210.0000.1850.001MSE0.0450.0570.1560.0140.7170.0031.6620.002PMBias0.2460.0060.4570.0080.8960.0041.4110.005MSE0.3610.0851.3730.0215.2380.00512.950.002MMBias0.0910.0120.0880.0050.1760.0010.2780.002MSE0.0550.0580.1560.0140.7470.0031.8140.00230CDFMBias0.0160.0150.0260.0010.0690.0010.0680.000MSE0.0310.0400.1220.0100.4980.0021.1210.001MLBias0.0500.0080.0960.0030.2060.0010.2790.002MSE0.0280.0380.1100.0090.4410.0020.9860.001LMBias0.0170.0070.0270.0040.0870.0000.1090.000MSE0.0270.0380.1000.0090.4300.0021.0030.001PMBias0.1410.0060.2710.0060.5700.0030.8230.003MSE0.1660.0550.6250.0142.6350.0035.8840.002MMBias0.0620.0090.0540.0040.1230.0000.1720.001MSE0.0340.0390.0990.0090.4370.0021.0750.00140CDFMBias0.0110.0150.0310.0010.0920.0010.1620.002MSE0.0230.0310.0950.0070.3770.0020.8570.001MLBias0.0410.0070.0870.0010.1880.0020.2860.002MSE0.0210.0290.0840.0070.3500.0020.7470.001LMBias0.0180.0030.0330.0020.1030.0010.1780.002MSE0.0190.0290.0770.0070.3280.0020.7630.001PMBias0.0950.0100.2270.0040.4420.0040.7480.001MSE0.1000.0440.5240.0111.6920.0033.8050.001MMBias0.0560.0100.0540.0020.1280.0010.2220.002MSE0.0260.0300.0770.0070.3320.0020.8130.00150CDFMBias0.0080.0070.0200.0010.0350.0000.0480.001MSE0.0170.0240.0700.0060.2880.0010.6220.001MLBias0.0280.0030.0600.0000.1180.0010.1680.001MSE0.0150.0220.0590.0060.2400.0010.5260.001LMBias0.0090.0060.0200.0010.0450.0010.0700.000MSE0.0150.0230.0560.0060.2450.0010.5570.001PMBias0.0810.0010.1670.0030.3290.0020.4850.002MSE0.0790.0340.3130.0091.2370.0022.7180.001MMBias0.0380.0040.0350.0010.0650.0010.1060.001MSE0.0190.0230.0550.0060.2440.0010.5910.00160CDFMBias0.0100.0150.0080.0000.0320.0010.0290.000MSE0.0140.0190.0550.0050.2250.0010.4960.001MLBias0.0280.0090.0400.0020.1000.0020.1250.000MSE0.0120.0180.0430.0050.1890.0010.4020.001LMBias0.0150.0030.0080.0030.0430.0010.0450.001MSE0.0120.0180.0430.0050.1960.0010.4330.001PMBias0.0650.0120.1240.0020.2370.0030.416−0.002MSE0.0550.0290.2540.0070.9070.0022.0440.001MMBias0.0400.0120.0200.0020.0610.0020.0750.000MSE0.0170.0190.0410.0050.1960.0010.4530.00170CDFMBias0.0060.0080.0100.0000.0520.0020.0460.000MSE0.0110.0170.0490.0040.2170.0010.4520.000MLBias0.0210.0040.0350.0010.1080.0020.1340.001MSE0.0090.0160.0390.0040.1800.0010.3740.000LMBias0.0070.0020.0080.0020.0580.0020.0620.000MSE0.0100.0160.0390.0040.1840.0010.3970.000PMBias0.0580.0030.1180.0040.2330.0030.3630.002MSE0.0480.0260.2060.0060.7870.0011.8070.001MMBias0.0280.0050.0170.0020.0720.0020.0880.000MSE0.0140.0170.0380.0040.1830.0010.4150.000Table 2
The bias and MSE of five estimation methods forβ=2, α = 1,2,4,6.
β=2α=1α=2α=4α=6nMethodα^iβ^iα^iβ^iα^iβ^iα^iβ^i10CDFMBias0.0460.0790.1070.0140.2030.0030.2850.004MSE0.1310.4820.5070.1162.1000.0294.5350.013MLBias0.1690.0360.3410.0090.6660.0140.9800.011MSE0.1540.4610.6110.1092.4110.0285.3710.012LMBias0.0550.0550.1140.0160.2640.0060.4140.004MSE0.1120.4600.4350.1111.8900.0284.3620.012PMBias0.5060.0371.0040.0252.0010.0262.8250.020MSE1.3100.6355.1460.15122.200.03946.900.017MMBias0.1670.0420.1940.0120.3870.0090.6290.007MSE0.1340.4610.4540.1102.0720.0285.0720.01220CDFMBias0.0220.0330.0420.0060.1070.0010.1310.002MSE0.0500.2440.1960.0570.8180.0141.7860.007MLBias0.0740.0090.1490.0060.3120.0060.4540.005MSE0.0480.2280.1910.0540.8040.0141.7540.006LMBias0.0250.0350.0440.0090.1290.0020.1880.002MSE0.0430.2300.1610.0550.7220.0141.6480.006PMBias0.2320.0110.4470.0130.9290.0131.3570.009MSE0.3510.3421.3200.0815.4240.02112.000.009MMBias0.0880.0140.0840.0070.1830.0040.2850.003MSE0.0530.2320.1600.0550.7510.0141.8120.00630CDFMBias0.0150.0260.0310.0050.0550.0000.0960.000MSE0.0300.1600.1230.0400.4850.0101.1260.004MLBias0.0490.0110.0990.0020.1950.0040.2980.002MSE0.0280.1490.1100.0380.4340.0091.0050.004LMBias0.0160.0190.0310.0040.0710.0020.1260.000MSE0.0270.1510.1000.0380.4140.0091.0040.004PMBias0.1380.0100.2820.0080.5610.0080.8870.005MSE0.1590.2210.6400.0562.7120.0146.0700.006MMBias0.0620.0150.0560.0030.1050.0030.1850.001MSE0.0340.1540.0990.0380.4190.0091.0770.00440CDFMBias0.0140.0040.0190.0030.0310.0060.0580.000MSE0.0200.1210.0860.0270.3180.0070.8270.003MLBias0.0400.0120.0670.0020.1370.0020.2110.002MSE0.0200.1140.0730.0260.2670.0070.7040.003LMBias0.0120.0360.0190.0030.0410.0040.0910.000MSE0.0190.1170.0700.0260.2670.0070.7540.003PMBias0.1090.0240.1850.0080.4120.0010.5860.004MSE0.0950.1720.3750.0401.5770.0103.6990.005MMBias0.0450.0140.0380.0030.0660.0040.1350.001MSE0.0250.1180.0690.0260.2650.0070.8000.00350CDFMBias0.0120.0190.0170.0040.0350.0010.0620.001MSE0.0180.0970.0710.0230.2790.0060.6190.003MLBias0.0310.0100.0580.0010.1160.0030.1750.002MSE0.0150.0910.0600.0220.2340.0060.5280.002LMBias0.0110.0080.0170.0020.0450.0020.0850.000MSE0.0150.0920.0570.0220.2380.0060.5590.002PMBias0.0860.0060.1600.0050.3140.0060.4910.005MSE0.0760.1370.3080.0331.1990.0092.7270.004MMBias0.0400.0120.0330.0010.0650.0030.1230.001MSE0.0200.0950.0560.0220.2380.0060.5960.00260CDFMBias0.0110.0130.0120.0050.0260.0030.0440.002MSE0.0150.0760.0570.0190.2280.0050.5460.002MLBias0.0250.0050.0520.0010.0930.0050.1460.003MSE0.0120.0690.0460.0180.1890.0050.4530.002LMBias0.0110.0090.0160.0000.0360.0040.0610.001MSE0.0120.0710.0450.0180.1970.0050.4730.002PMBias0.0750.0090.1260.0000.2680.0070.4200.003MSE0.0620.1140.2370.0280.9660.0072.2660.003MMBias0.0350.0080.0300.0000.0530.0040.0910.002MSE0.0160.0730.0440.0180.1970.0050.4930.00270CDFMBias0.0060.0130.0060.0040.0240.0040.0660.002MSE0.0130.0680.0460.0150.2000.0040.4610.002MLBias0.0210.0050.0350.0010.0780.0050.1460.000MSE0.0100.0640.0370.0140.1550.0040.3890.002LMBias0.0090.0070.0060.0000.0360.0040.0780.002MSE0.0110.0640.0360.0140.1700.0040.4170.002PMBias0.0590.0090.1090.0040.1940.0090.3240.000MSE0.0500.0970.1860.0220.7370.0061.5150.003MMBias0.0300.0080.0160.0000.0510.0040.1000.001MSE0.0140.0670.0360.0140.1690.0040.4410.002According to the above findings, the bias ofα^CDFM outperforms other methods. Although the bias of β^CDFM is worse than β^LM, the difference is not significant. Therefore, using the proposed estimation method for inference can get reliable results.In addition, in order to compare the performance of the RF and FRF based on different parameter estimation methods, we consider that the sample sizen = 10, 20, 30, 40, 50, 60, 70, α=1,2,4, 6, β=1, and strength t = 2. All the simulation results are shown in Figures 3 and 4. In view of Figures 3 and 4, the following results are found:(1)
The performance of biases of the RF and FRF estimated by CDFM method outperforms other methods in most cases.(2)
With the increase ofα or sample size n, the performance of biases of each estimation method is better for RF and FRF.(3)
Among all estimation methods, PM method performs the worst for biases of RF and FRF.(4)
Among all estimation methods, PM method also performs the worst for MSEs of RF and FRF.(5)
With the increase ofα or sample size n, the performance of MSEs of each estimation method is better for RF and FRF.(6)
Except PM method, the performance of MSEs of other estimation methods is very close for RF and FRF.Figure 3
Comparison of biases of the RF and FRF based on CDFM and other four different parameter estimation methods (forβ = 1, α = 1, 2, 4, 6).Figure 4
Comparison of MSEs of the RF and FRF based on CDFM and other four different parameter estimation methods (forβ = 1, α = 1, 2, 4, 6).


## 5. Real Data Analysis

### 5.1. Example 1

The tensile strength data of glass fiber, reported by Smith and Naylor [28] and further studied by Acitas et al. [3], Morais and Barreto-Souza [29], Alizadeh et al. [30], and Jones and Faddy [31], is used to illustrate the methodologies developed in this paper. In the article of Smith and Naylor [28], two experimental strength data sets of glass fiber of two lengths, 1.5 cm and 15 cm, are provided, originally obtained by workers at the UK National Physical Laboratory. Unfortunately, the units of measurement are not given in the paper. In this example, the data of 1.5 cm length is selected for analysis. The strength data is shown in Table 3.Table 3
Sample data set-1 of glass fiber strength (1.5 cm).
0.550.931.251.361.491.521.581.611.641.681.731.812.000.741.041.271.391.491.531.591.611.661.681.761.822.010.771.111.281.421.501.541.601.621.661.691.761.842.240.811.131.291.481.501.551.611.621.661.701.771.840.841.241.301.481.511.551.611.631.671.701.781.89Two-parameter Weibull distribution is used to characterize the data in Table3. Moreover, the CDFM, MLE, LM, PM, and MM estimation methods are used to obtain the estimates of the Weibull distribution parameters. All the estimation results are reported in Table 4. The estimators of LM are consistent with those of Acitas et al. [3]. The p value of the K–S test in Table 4 indicates that the Weibull distribution can well fit this data set and we found the CDFM estimation method is the best estimation method to obtain the estimates of the Weibull distribution parameters. In addition, Figure 5 presents the histogram of sample and the fitted PDFs of the Weibull distributions by using the parameter estimation results in Table 4. Figure 6 presents empirical distribution function of sample and the fitted CDFs of the Weibull distributions by using the parameter estimation results in Table 4. As shown in Figure 5, the PDF corresponding to CDFM is closest to the sample histogram, which also shows that the estimation result of CDFM is more suitable for this data. From Figure 6, we can also see that the curve corresponding to CDFM is closest to the empirical distribution function of the sample. The results in Figure 6 also show that the CDFM is more suitable for this data.Table 4
Estimate of the parameters and theP value of K–S test for simple data set-1.
Estimate ofαEstimate ofβP value of K-S testCDFM6.63811.63900.3980ML5.78071.62810.1001LM5.65731.62960.0863PM10.30261.63190.1254MM5.35541.63470.0802Figure 5
The histogram and the fitted PDFs of the Weibull distributions for sample data set-1.Figure 6
Empirical distribution function and the fitted CDFs of the Weibull distributions for sample data set-1.Suppose we want to estimate the reliability and failure rate of glass fiber, we can substitute equations (17) and (18) according to the parameter estimator results in Table 4, and the results are shown in Figure 7. From Figure 7, the results show that when t < 1.6, the RF based on CDFM is higher than other methods except PM method; when t > 1.6, the RF based on CDFM is lower than other methods except PM method. In addition, the FRF based on CDFM is higher than other methods except PM method.Figure 7
The RF and FRF based on different parameter estimation methods for sample data set-1.


### 5.2. Example 2

In this example, the data of 15 cm length reported by Smith and Naylor [28] is selected for analysis. The strength data is shown in Table 5. Two-parameter Weibull distribution is also used to characterize the data in Table 5. Moreover, the CDFM, MLE, LM, PM, and MM estimation methods are used to obtain the estimates of the Weibull distribution parameters. All the estimation results are reported in Table 6. The p value of the K-S test in Table 6 indicates that the Weibull distribution can well fit this data set and we found the CDFM estimation method is the best estimation method to obtain the estimates of the Weibull distribution parameters. In addition, Figure 8 presents the histogram of sample and the fitted PDFs of the Weibull distributions by using the parameter estimation results in Table 6. Figure 9 presents empirical distribution function of sample and the fitted CDFs of the Weibull distributions by using the parameter estimation results in Table 6. The results in Figures 8 and 9 show that the fitting performance of these estimation methods for this sample is similar.Table 5
Sample data set-2 of glass fiber strength (15 cm).
0.370.400.700.750.800.810.830.860.920.920.940.950.981.031.061.061.081.091.101.101.131.141.151.171.201.201.211.221.251.281.281.291.291.301.351.351.371.371.381.401.401.421.431.511.531.61Table 6
Estimate of the parameters and theP value of K–S test for sample data set-2.
Estimate ofαEstimate ofβP value of K-S testCDFM5.05621.24320.9792ML5.14731.22970.9302LM4.82651.23320.9354PM5.60341.26340.9183MM4.74721.23440.9011Figure 8
The histogram and the fitted PDFs of the Weibull distributions for sample data set-2.Figure 9
Empirical distribution function and the fitted CDFs of the Weibull distributions for sample data set-2.Suppose we want to estimate the reliability and failure rate of glass fiber, we can substitute equations (17) and (18) according to the parameter estimator results in Table 6, and the results are shown in Figure 10. It can be seen from Figure 10 that the RF and FRF results estimated by each method are almost the same, because the parameter estimator results of this data set for each estimation method are almost the same.Figure 10
The RF and FRF based on different parameter estimation methods for sample data set-2.


## 5.1. Example 1

The tensile strength data of glass fiber, reported by Smith and Naylor [28] and further studied by Acitas et al. [3], Morais and Barreto-Souza [29], Alizadeh et al. [30], and Jones and Faddy [31], is used to illustrate the methodologies developed in this paper. In the article of Smith and Naylor [28], two experimental strength data sets of glass fiber of two lengths, 1.5 cm and 15 cm, are provided, originally obtained by workers at the UK National Physical Laboratory. Unfortunately, the units of measurement are not given in the paper. In this example, the data of 1.5 cm length is selected for analysis. The strength data is shown in Table 3.Table 3
Sample data set-1 of glass fiber strength (1.5 cm).
0.550.931.251.361.491.521.581.611.641.681.731.812.000.741.041.271.391.491.531.591.611.661.681.761.822.010.771.111.281.421.501.541.601.621.661.691.761.842.240.811.131.291.481.501.551.611.621.661.701.771.840.841.241.301.481.511.551.611.631.671.701.781.89Two-parameter Weibull distribution is used to characterize the data in Table3. Moreover, the CDFM, MLE, LM, PM, and MM estimation methods are used to obtain the estimates of the Weibull distribution parameters. All the estimation results are reported in Table 4. The estimators of LM are consistent with those of Acitas et al. [3]. The p value of the K–S test in Table 4 indicates that the Weibull distribution can well fit this data set and we found the CDFM estimation method is the best estimation method to obtain the estimates of the Weibull distribution parameters. In addition, Figure 5 presents the histogram of sample and the fitted PDFs of the Weibull distributions by using the parameter estimation results in Table 4. Figure 6 presents empirical distribution function of sample and the fitted CDFs of the Weibull distributions by using the parameter estimation results in Table 4. As shown in Figure 5, the PDF corresponding to CDFM is closest to the sample histogram, which also shows that the estimation result of CDFM is more suitable for this data. From Figure 6, we can also see that the curve corresponding to CDFM is closest to the empirical distribution function of the sample. The results in Figure 6 also show that the CDFM is more suitable for this data.Table 4
Estimate of the parameters and theP value of K–S test for simple data set-1.
Estimate ofαEstimate ofβP value of K-S testCDFM6.63811.63900.3980ML5.78071.62810.1001LM5.65731.62960.0863PM10.30261.63190.1254MM5.35541.63470.0802Figure 5
The histogram and the fitted PDFs of the Weibull distributions for sample data set-1.Figure 6
Empirical distribution function and the fitted CDFs of the Weibull distributions for sample data set-1.Suppose we want to estimate the reliability and failure rate of glass fiber, we can substitute equations (17) and (18) according to the parameter estimator results in Table 4, and the results are shown in Figure 7. From Figure 7, the results show that when t < 1.6, the RF based on CDFM is higher than other methods except PM method; when t > 1.6, the RF based on CDFM is lower than other methods except PM method. In addition, the FRF based on CDFM is higher than other methods except PM method.Figure 7
The RF and FRF based on different parameter estimation methods for sample data set-1.


## 5.2. Example 2

In this example, the data of 15 cm length reported by Smith and Naylor [28] is selected for analysis. The strength data is shown in Table 5. Two-parameter Weibull distribution is also used to characterize the data in Table 5. Moreover, the CDFM, MLE, LM, PM, and MM estimation methods are used to obtain the estimates of the Weibull distribution parameters. All the estimation results are reported in Table 6. The p value of the K-S test in Table 6 indicates that the Weibull distribution can well fit this data set and we found the CDFM estimation method is the best estimation method to obtain the estimates of the Weibull distribution parameters. In addition, Figure 8 presents the histogram of sample and the fitted PDFs of the Weibull distributions by using the parameter estimation results in Table 6. Figure 9 presents empirical distribution function of sample and the fitted CDFs of the Weibull distributions by using the parameter estimation results in Table 6. The results in Figures 8 and 9 show that the fitting performance of these estimation methods for this sample is similar.Table 5
Sample data set-2 of glass fiber strength (15 cm).
0.370.400.700.750.800.810.830.860.920.920.940.950.981.031.061.061.081.091.101.101.131.141.151.171.201.201.211.221.251.281.281.291.291.301.351.351.371.371.381.401.401.421.431.511.531.61Table 6
Estimate of the parameters and theP value of K–S test for sample data set-2.
Estimate ofαEstimate ofβP value of K-S testCDFM5.05621.24320.9792ML5.14731.22970.9302LM4.82651.23320.9354PM5.60341.26340.9183MM4.74721.23440.9011Figure 8
The histogram and the fitted PDFs of the Weibull distributions for sample data set-2.Figure 9
Empirical distribution function and the fitted CDFs of the Weibull distributions for sample data set-2.Suppose we want to estimate the reliability and failure rate of glass fiber, we can substitute equations (17) and (18) according to the parameter estimator results in Table 6, and the results are shown in Figure 10. It can be seen from Figure 10 that the RF and FRF results estimated by each method are almost the same, because the parameter estimator results of this data set for each estimation method are almost the same.Figure 10
The RF and FRF based on different parameter estimation methods for sample data set-2.


## 6. Conclusions

This study aims to solve the problem of parameter estimation for the Weibull distribution and then infer the reliability of products.


### 6.1. Unique Contributions of the Study

We propose a new parameter estimation method (CDFM) and use this method to estimate the reliability of products. According to Monte Carlo simulation results, it can be found that the CDFM estimation method outperforms the other competitors to obtain reliable estimates of two-parameter Weibull parameters. Therefore, two strength data sets of glass fiber provided by Smith and Naylor [28] are estimated according to the proposed estimation method. The analysis results show that the estimates of CDFM are more suitable for these data sets. And according to the estimates of CDFM, the reliability and failure rate of glass fiber are also estimated. From the analysis results, it can be seen that the reliability and failure rate estimated by different parameter estimates will be different, so it is very important to use a reliable parameter estimate for inference. In addition, the estimation method proposed in this paper can be applied to other fibers, materials, or tests as long as the data follow the Weibull distribution.


### 6.2. Limitations

In particular, the proposed parameter estimation method can be used as long as the data follows the Weibull distribution, but it cannot guarantee the same performance when the data follows the other distributions.


### 6.3. Future Research

In the current work, we only consider the parameter estimation and application of two-parameter Weibull distribution without inferring the theoretical properties and limit distribution of the estimators. In future research, in addition to applying this estimation method to other distributions, we can also infer the limit distribution and statistical properties of these estimators. We are currently exploring possible extensions in this direction.


## 6.1. Unique Contributions of the Study

We propose a new parameter estimation method (CDFM) and use this method to estimate the reliability of products. According to Monte Carlo simulation results, it can be found that the CDFM estimation method outperforms the other competitors to obtain reliable estimates of two-parameter Weibull parameters. Therefore, two strength data sets of glass fiber provided by Smith and Naylor [28] are estimated according to the proposed estimation method. The analysis results show that the estimates of CDFM are more suitable for these data sets. And according to the estimates of CDFM, the reliability and failure rate of glass fiber are also estimated. From the analysis results, it can be seen that the reliability and failure rate estimated by different parameter estimates will be different, so it is very important to use a reliable parameter estimate for inference. In addition, the estimation method proposed in this paper can be applied to other fibers, materials, or tests as long as the data follow the Weibull distribution.


## 6.2. Limitations

In particular, the proposed parameter estimation method can be used as long as the data follows the Weibull distribution, but it cannot guarantee the same performance when the data follows the other distributions.


## 6.3. Future Research

In the current work, we only consider the parameter estimation and application of two-parameter Weibull distribution without inferring the theoretical properties and limit distribution of the estimators. In future research, in addition to applying this estimation method to other distributions, we can also infer the limit distribution and statistical properties of these estimators. We are currently exploring possible extensions in this direction.



---

*Source: 9175170-2021-09-22.xml*