A new trigonometric modification of the Weibull distribution: Control chart and applications in quality control

In the most recent era, the extensions of the probability models via trigonometry methods have received great attention. This paper also offers a novel trigonometric version of the Weibull model called a type-I cosine exponentiated Weibull (for short “TICE-Weibull”) distribution. The identifiability properties for all three parameters of the TICE-Weibull model are derived. The estimators of the TICE-Weibull model are derived by implementing the maximum likelihood approach. To demonstrate the effectiveness of the TICE-Weibull model, two applications from real-world phenomena are analyzed. In addition, the proposed statistical model is established for an attribute control chart based on a time-truncated life test. The advantage of the developed charts is examined based on the average run length (ARL). The necessary tables of shift sizes and various sample sizes are offered for numerous values of the distribution parameters, as well as specified ARL and shift constants. Some numerical examples are discussed for various scheme parameters to study the performance of the new TICE-Weibull attribute control charts. According to our search and a brief study of the statistical literature, there is no published work on the development of a control chart using new probability models that are introduced using the cosine function. This is the key motivation of this work, which fills this amazing and interesting research gap.


Introduction
The two parameters (α > 0, δ > 0) Weibull model is one of the most famous and important probability models that can best describe the failure behavior of the system/components during its lifetime. It is defined on R þ and has been implemented for analyzing lifetime phenomena in engineering and other close connected fields (Pham and Lai [1], Almalki and Nadarajah [2]).
The up-gradation of the Weibull model has been done by implementing numerous approaches; see Ahmad et al. [13]. However, the literature on distribution theory has lack probability distributions that are based on trigonometric functions; most of these modifications are based on algebraic functions. The increased interest in data analysis, modeling, predicting, and directional data analysis, in particular, motivates researchers to look for the development of new approaches based on trigonometric functions. Therefore, in the recent literature about the development of new updated versions of the Weibull model, researchers have focused on implementing trigonometric functions; see Chesneau et al. [14].
The one explored by Souza et al. [15] provides a modern alternative based on a trigonometric function. They introduced a new class of distributions called a sin-G class of probability models with cumulative distribution function (CDF) F(x; ζ) given by where G(x; ζ) is CDF of the baseline random variable depending on the vector of parameters ζ. Silveira et al. [16] introduced the normal-tangent-G (NT-G) class of distributions for studying new statistical models. The CDF F(x; ζ) of the NT-G distributions is given by Another contribution work towards the trigonometric family of distributions is due to Jamal et al. [17]. They introduced a generalized version of the sin-G method called the transformed sin-G (TS-G) family of distributions. The CDF F(x; η, ζ) of the TS-G family is where η = 2 [0, 1]. Nanga et al. [18] introduced the tangent Topp-Leone (TTL) family of distributions having CDF F(x; β, ζ) given by where β > 0. From the above literature, we can see that researchers are turning their attention to introducing new probability distributions using trigonometric functions. However, based on our deep study of the literature, there is no published work about the construction of new control charts using new probability distributions that are based on the cosine function. In this paper, first, we use a trigonometric function to derive a new form of the Weibull distribution. Then, using the proposed model, we construct a new control chart. This is a key motivation for this paper.
The new form of the Weibull model is introduced by using the trigonometric function called Type-I cosine exponentiated-X (TICE-X) family. Suppose X has the TICE-X family, if its CDF is given by respectively.
In this paper, we implement the method presented in Eq (1) to define/introduce a new trigonometric-Weibull (TICE-Weibull) distribution. Section 2 offers the definition of the basic functions of the TICE-Weibull distribution. The plots for the PDF of the TICE-Weibull distribution are also presented in Section 2. Some distributional properties of the TICE-Weibull distribution are obtained in Section 3. A simulation study in Section 4 is provided to evaluate the performances of the estimators of the TICE-Weibull distribution using two well-known statistical criteria called (i) the average bias, and (ii) the average mean square error. The applicability of the TICE-Weibull distribution is shown by considering two real-life applications in Section 5. An attribute control chart based on truncated life tests is given in Section 6 and in Section 7, the final remarks are provided.

Mathematical properties
This section of the paper deals with the computation of some mathematical properties of the TICE-Weibull distribution. These properties include the quantile function (QF), identifiability property (IP), and m th moment.

The quantile function
The quantile of the TICE-Weibull distribution, say x q , is obtained by solving the following equation for x Fðx; y; ζ Þ ¼ q;

PLOS ONE
where 0 < q < 1. Using Eq (5) in Eq (7), we get On solving Eq (8), we get the QF of the TICE-Weibull distribution, which is given by

The IP of the TICE-Weibull model
In this subsection, we prove the IP of the TICE-Weibull distribution using the parameters α, δ, and θ.

The m th moment
Here, we obtain the m th moment of the TICE-Weibull model as follows Using Eq (2) in Eq (10), we have where A(x; θ, ζ) and B(x; θ, ζ) are given by and B x; y; ζ ð Þ ¼ respectively.

Estimation and simulation
This section has two aims. The very first aim concerns the derivation of the maximum likelihood estimators (MLEs) ðŷ MLE ; ζ MLE Þ of the parameters (θ, ζ) of the TICE-Weibull distribution.
The second aim of this section is to demonstrate the performances ofŷ MLE and ζ MLE :

Simulation
Here, we demonstrate the performances of the MLEs ðŷ MLE ;ζ MLE Þ obtained in subsection 4.1.

The demonstration ofŷ
The above formulas are also implemented to obtain the values of Bias and MSE for the parameter vector ζ MLE : Corresponding to α = 0.7, δ = 1.0, θ = 1.3, the simulation results (SRs) are provided in Table 1 and presented graphically in Fig 2. Whereas, the SRs for α = 1.1, δ = 1.0, θ = 0.9 are presented in Table 2 and illustrated visually in

Data modeling
This section is concerned with an illustration of the usefulness of the TICE-Weibull distribution for modeling the lifetime scenarios of electronic components. Both the data sets are taken
• The W* test is one of the most widely used measures to compare the fitting capability of the estimated models. For the underlined data, a model with the smaller value of W* test has the best fit. The W* test is defined as where k and x i denote the sample size and i th observation of the underlined data, respectively.

PLOS ONE
• The A* test is also used as a comparative measure for the competing models. This test is computed as • Let G k (x) andĜðxÞ, respectively, represent the estimated CDF of the selected model (given model) and empirical CDF of the sample values. Then, the KS* is computed by sup x h G k ðxÞ ÀĜðxÞ i :

Data 1
This data represents the lifetimes of fifty-nine electronic devices. Corresponding to the first data set, the key values are provided in Table 3. Some descriptive plots of the Data 1 are presented in

Data 2
Data 2 denote the lifetimes of twenty electronic components. The observations of Data 2 with the associated summary values are displayed in Table 6. Besides the key measures, certain descriptive plots of Data 2 are also obtained; see Fig 6. Based on the statistical analysis of Data 2, the MLEs and values of the selection criteria of TICE-Weibull and other candidate models are displayed in Tables 7 and 8, respectively. Based  on Table 8, we can see that the TICE-Weibull model again turns out to be the best suitable distribution.
Some fitted plots of the TICE-Weibull model are shown in Fig 7. From these plots (see Fig  7), it can be observed that the estimated CDF and PDF of TICE-Weibull distribution closely follow the empirical CDF and the pattern of the histogram on Data 2, respectively. This behavior of the TICE-Weibull has also been confirmed by obtaining the PP, empirical SF, and QQ plots.

Attribute control chart
There are assorted valuable control charts to maintain the manufacturing process under control, comprising the lie between the lower control limit (LCL) and upper control limit (UCL). In the statistical quality control process, there are two renowned types of control charts, namely, the control chart for attributes and the control chart for variables. The variables control chart provides good information considering the process and comprises minimum sample sizes; since it applies quantitative data. The use of an attribute control chart is more flexible as compared to the variables chart due to its easiness of computation. The various attribute control chart are renowned in the literature such as the np chart, u chart, and the c chart.

PLOS ONE
The attribute control chart for time truncated life test for various distributions has been studied by various authors, please refer to Haq and Al-Omari [19], Quinino et al. [20], Rao and Al-Omari [21] and Adeoti and Rao [22]. By exploring the literature there is no work on an attribute control chart based on the TICE-Weibull distribution. In this paper, a new attribute control chart based on a truncated life test is projected using the TICE-Weibull distribution.

The proposed control chart
We project a new np control chart founded on time-truncated life testing as developed by Aslam and Jun [23]: • Step 1 Chose a simple random sample of size n from the submitted lot and examine them. The number of failures dented by φ is obtained before the experiment time x 0 = τξ 0 , where ξ 0 is

PLOS ONE
the quality consideration under the condition that the process is in-control and τ is a multiplier constant.
• Step 2 Declare that the process as out-of-control when φ > UCL and φ < LCL otherwise the process is in-control if LCL < φ < UCL.
The percentile of the T-Weibull distribution is Let Using the binomial distribution of defective products with parameters p 0 and n the proposed chart limits are, respectively, obtained as where p 0 is the probability of failed article before the testing time x 0 when the process considered as in-control, and L is the chart coefficient to be obtained. On the other hand, we can say the process is in-control once ξ = ξ 0 (i.e., α = α 0 , θ = θ 0 , δ = δ 0 ). Let us consider that the experiment time is x 0 as multiple of termination ratio B q and specified percentile life x q0 , i.e., B q x q0 in time-truncated lifetime experimentation. After simplification, the probability of failure is written as When the process is in-control, then the percentile ratio x q x q0 ¼ 1: Therefore, Eq (16) will be reduced to

PLOS ONE
Now, we consider the percentile ratio x q /x q0 = c = 1.0, 1.05, 1.10, . . ., 4.0 Then, the probability in Eq (17) becomes Let � φ denote the average of failures for the subgroups over the sample. If the value of p 0 is not known, the chart limits for realistic purposes can be used in the following expressions ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi � φð1 À � φ=nÞ p ; and ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi � φð1 À � φ=nÞ p g: The chance of declaring as the process is in-control for the developed control chart is given by The accomplishment of the developed control chart can be examined by its average control length (ARL) and when the process as in-control state it is expressed as follows: To investigate the performance of the proposed control chart, the study of out-of-control is needed to investigate. When the process is out-of-control, assume that p 1 as the probability of an unsuccessful item earlier than the experiment time x 0 Hence, the chance that the process is evident to be in control, whilst the declared time ratio is changed to c is given by The following step-by-step procedure can be used to acquire the tables of the developed control chart • Find out the ARL value, say r 0 , and known parametric values α = α 0 , θ = θ 0 , δ = δ 0 , α 0 .
• Determine the chart constants L, τ, and n such that the ARL 0 value is almost equal to r 0 , i.e., ARL 0 > r 0 .
• Subsequent to receiving the values in the above step, determine the ARL 1 according to shift constant c based on Eq (18).
We determined the control chart parameters and ARL 1 for various values of α = α 0 , θ = θ 0 , r 0 , and n, given in Tables 9-12 for shift values. Based on computed tables, we observed the following conclusions • It is observed that ARL 1 value in decreasing tendency as the shift value c increases.
• From Tables 9 and 10, it is evident that for fixed θ value ARL1 values are decreases with the increase of parameter α.
• Based on Tables 9-12, it is reasonable that when the sample size n increases from 20 to 40, the ARL 1 values are shows decreasing tendency.

Illustration of the suggested control chart
The demonstration of the developed control chart is as follows: let us assume that industrial output persists the TICE-Weibull distribution with parameters α = 0.7, θ = 1.3. Suppose the average target lifetime of the product is ξ 0 = 1000 hours and r 0 = 370. Using the Eq (17) the value of p 0 is 0.48866. Also, from Table 13, the chart parameters are n = 20, τ = 0.968, L = 3.029, LCL = 3, and UCL = 16. Hence, the experiment time x 0 is 968 hours. Thus, the proposed control chart is carried out in the following steps: • Step I Draw a simple random sample of size 20 from every subgroup and put them for the life testing assessment during 968 hours. Determine the number of failed units say D for the duration of the experiment time.

PLOS ONE
• Step II Pronounce the production process as under control if 3 � D � 16; else, the production process can be regarded as out-of-control.

An industrial application
Using the second data set, the estimated parameters of the TICE-Weibull distribution arê y ¼ 0:3651;â ¼ 2:2549; andd ¼ 0:02901: Table 13 provides the ARLs of the proposed control chart for the second data set. Here, n ¼ 20;ŷ ¼ 0:3651;â ¼ 2:2549; t ¼ 1:064; and r 0 = 300. The value of p 0 is 0.5281 using Eq (17). The value of ξ by Eq (15) Fig 8. From Fig 8, it is noticed that the proposed chart shows lifetimes of electronic components data are under control. Hence, the proposed chart is suitable to monitor the quality of the electronic components.

Comparison
The ARL values of the proposed control chart and the existing time truncated life testing attributed control charts for the Weibull distribution given in Adeoti and Rao (2021) are compared. The results of comparisons between the TICE-Weibull and Weibull distributions are displayed in Table 14. To compare the two types of control charts with respect to ARL values at different shift values. It is important to note that a chart having lesser out-of-control ARLs would be considered the better control chart. We noticed based on

Final remarks
In this study, a new trigonometric version of Weibull distribution with the implementation of a trigonometric function was investigated. The TICE-Weibull model can be considered a

PLOS ONE
useful updated version of the Weibull model. Certain statistical properties including shapes of PDF along with the IP, QF, and m th moment have been studied. The MLEs of the parameters involved in the TICE-Weibull distribution are obtained. Two real datasets were considered for illustration of the TICE-Weibull distribution. Moreover, an attribute control chart is proposed for TICE-Weibull distribution. The proposed control chart coefficient and out-of-control ARLs for various sample sizes and parametric values are given in Tables 9-12. The results show that the proposed control chart shows its efficiency at different sample sizes as compared with the existing Weibull distribution. Finally, the proposed control chart was also demonstrated with the help of real data for industrial application. Based on our findings, it is observed that the TICE-Weibull distribution control chart is better than the Weibull distribution control chart for monitoring the lifetime of electronic components.
In the future, we are motivated to apply the proposed distribution to the dynamical systems via data-driven approaches such as the non-Gaussian system control. We are also intended to obtain the bivariate extension of the proposed model for analyzing the bivariate data sets. Furthermore, the proposed control chart can also be implemented to monitor other industrial processes like chemical, electrical or mechanical energy in the manufacturing of industrial products or pharmaceutical manufacturing processes using other models like repetitive, multiple dependent state sampling etc.