Inspection plan for COVID-19 patients for Weibull distribution using repetitive sampling under indeterminacy

Abstract

Background

This research work is elaborated investigation of COVID-19 data for Weibull distribution under indeterminacy using time truncated repetitive sampling plan. The proposed design parameters like sample size, acceptance sample number and rejection sample number are obtained for known indeterminacy parameter.

Methods

The plan parameters and corresponding tables are developed for specified indeterminacy parametric values. The conclusion from the outcome of the proposed design is that when indeterminacy values increase the average sample number (ASN) reduces.

Results

The proposed repetitive sampling plan methodology application is given using COVID-19 data belong to Italy. The efficiency of the proposed sampling plan is compared with the existing sampling plans.

Conclusions

Using the tables and COVID-19 data illustration, it is concluded that the proposed plan required a smaller sample size as examined with the available sampling plans in the literature.

Background

It is broadly established that a huge number of COVID-19 cases are unnoticed worldwide. A rudimentary measure of population occurrence is the small part of positive cases for a given date in any country. On the other hand, this is subject to largely found that bias since tests are normally only ordered from suggestive cases, whereas a large proportion of infected people might show little symptoms or sometimes no symptoms for more details see [1]. Most governments are applying the mechanism of test randomly selected individuals to estimate the true disease occurrence in inhabitants in a particular locality. Nevertheless, when the disease occurrence is low and difficult to acquire from each patient/person by tests, under such situations we may use an acceptance sampling plan under indeterminacy. The health practitioners are paying attention to estimate the average number of deaths or ratio of deaths to the total number of COVID-19 death cases on daily basis, for the coming days, next week or month, etc. Reference [2]. In such a case, the health practitioners are paying attention to test the null hypothesis that the average number of deaths or ratio of deaths to the total number of COVID-19 death cases on daily basis is equal to the specified average number of deaths due to COVID-19 against the alternative hypothesis that the average number of deaths due to COVID-19 varies significantly. In this situation for testing of the hypothesis, practically it is difficult to record the average number of death for the whole year, whereas it is easy to record the daily basis and the average number of deaths can be obtained from the randomly selected days. The null hypothesis may be rejected if the daily average number of deaths due to COVID-19, state acceptance number of days, is more than or equal to the specified average number of deaths due to COVID-19 throughout the given number of days.

Many researchers have done studies on the time truncated life test for various distributions. Some of them are [3] developed the acceptance sampling plan for life tests: log-logistic models. Reference [4] derived acceptance sampling based on truncated life tests for generalized Rayleigh distribution. Reference [5] developed the acceptance sampling plans based on the generalized Birnbaum-Saunders distribution. Reference [6, 7] constructed the acceptance sampling plans for Birnbaum-Saunders and Burr XII distributions. References [8, 9] constructed acceptance sampling plans for extended exponential and generalized inverted exponential distributions. The details about the acceptance sampling plans can be seen in [10, 11]. The generalization of a single acceptance sampling plan namely repetitive sampling plan, [12] derived the decision rule of the repetitive acceptance sampling plan. The method of repetitive group acceptance sampling plan (RGASP) was first proposed by [13] for an attribute. Reference [14, 15] constructed the RASP for inverse Gaussian distribution and Burr type XII. Reference [16] developed generalized inverted exponential distributions. References [17,18,19] studied the repetitive sampling plan under different situations.

More details about the neutrosophic logic, their measure of determinacy, and indeterminacy are given by [20]. Numerous authors studied the neutrosophic logic for various real problems and showed its efficiency over fuzzy logic, for more details refer [21,22,23,24,25,26]. The idea of neutrosophic statistics was given using the idea of neutrosophic logic, [27,28,29]. The neutrosophic statistics give information about the measure of determinacy and measure of indeterminacy, see [30]. The neutrosophic statistics reduce to classical statistics if no information is recorded about the measure of indeterminacy. References [31,32,33] proposed the acceptance sampling plans using neutrosophic statistics [34]. proposed the time-truncated group plans for the Weibull distribution. Reference [35] worked on neutrosophic Weibull and neutrosophic family of Weibull distribution.

The existing sampling plans based on classical statistics and fuzzy philosophies do not give information about the measure of indeterminacy. Reference [36] worked on the single sampling plan using a fuzzy approach. Reference [37] discussed the effect of sampling error on inspection using a fuzzy approach. Reference [38] proposed a single plan using fuzzy logic. Reference [39] proposed the improved sampling plan using fuzzy logic. For details, the reader may refer to [40, 41]. To the best of our knowledge, there is no work on a time-truncated sampling plan for Weibull distribution under indeterminacy. In this paper, a repetitive acceptance sampling plan for Weibull distribution under indeterminacy is developed to testing the daily average deaths. We are anticipated the proposed sampling plan shows a fewer sample size as compared with the existing sampling plans for testing the daily average deaths.

In Section 2, we present an introduction of a repetitive acceptance sampling plan for Weibull distribution under indeterminacy. In Section 3, the proposed repetitive acceptance sampling plan under indeterminacy is compared with the single sampling plan proposed by [42]. The proposed sampling plan is illustrated using COVID-19 data belong to Italy, which was recorded from 1 April to 20 July 2020 in Section 4. Finally, the conclusions and future research works are established in Section 5.

Methods

The repetitive acceptance sampling plan depends upon the truncated life test procedure is developed by [43,44,45]. The operational steps of this test are given as follows:

• Step 1: Draw a sample of size n from the lot. These samples can be put on a life test for a fixed time t₀. Specify the average μ₀ and indeterminacy parameter I_Nϵ[I_L, I_U].

• Step 2: Accept H₀ : μ_N = μ_0N if the daily average deaths in c₁ days are more than or equal to μ₀ (i.e., μ₀ ≤ c₁). If daily average deaths in c₂ days are less than to μ₀ (i.e., μ₀ >c₂) then we reject H₀ : μ_N = μ_0N and terminate the test, where c_1≤ c₂.

• Step 3: If c₁< μ₀ ≤c₂ then go to step 1 and repeat the above experiment.

The above procedure of repetitive acceptance sampling plan (RASP) mainly depends on four characteristics those are n, c₁, c₂ and I_N, where I_Nϵ[I_L, I_U] is considered as the specified parameter and set according to the uncertainty level. RASP is nothing but the generalization of an ordinary single sampling plan under uncertainty. If c₁ = c₂ in RASP, it ultimately reduces to a single sampling plan under uncertainty. Suppose that t₀ = aμ₀ be the time in days, where a is the termination ratio. The lot acceptance probability is to be determined with the help of operating characteristic (OC) function for details see [13] and it is given by

$$L(p)=\frac{P_a(p)}{P_a(p)+{P}_r(p)};0<p<1$$

(1)

Here P_a(p) is the probability of accepting H₀ : μ_N = μ_0N and P_r(p) is the probability of rejecting H₀ : μ_N = μ_0N, which are given by

$${P}_a(p)=\sum \limits_{i=0}^{c_1}\left(\begin{array}{*{20}l} n\\ {}i\end{array}\right){p}^i{\left(1-p\right)}^{n-i}$$

(2)

and

$${P}_r(p)=1-\sum \limits_{i=0}^{c_2}\left(\begin{array}{*{20}l} n\\ {}i\end{array}\right){p}^i{\left(1-p\right)}^{n-i}$$

(3)

where p is the probability of unreliability.

Therefore eq. (1) becomes

$$L(p)=\frac{\sum \limits_{i=0}^{c_1}\left(\begin{array}{*{20}l} n\\ {}i\end{array}\right){p}^i{\left(1-p\right)}^{n-i}}{\sum \limits_{i=0}^{c_1}\left(\begin{array}{*{20}l} n\\ {}i\end{array}\right){p}^i{\left(1-p\right)}^{n-i}+1-\sum \limits_{i=0}^{c_2}\left(\begin{array}{*{20}l} n\\ {}i\end{array}\right){p}^i{\left(1-p\right)}^{n-i}};0<p<1$$

(4)

The Weibull distribution under neutrosophic statistics is developed by [42] for the design of the sampling scheme plan for testing the average wind speed under an indeterminate environment.

Suppose that f(x_N) = f(x_L) + f(x_U)I_N; I_Nϵ[I_L, I_U] be a neutrosophic probability density function (npdf) having determinate part f(x_L), indeterminate part f(x_U)I_N and indeterminacy interval I_Nϵ[I_L, I_U]. Note that x_Nϵ[x_L, x_U] be a neutrosophic random variable follows the npdf. The npdf is the generalization of pdf under classical statistics. The proposed neutrosophic form of f(x_N)ϵ[f(x_L), f(x_U)] reduces to pdf under classical statistics when I_L =0. Based on this information, the npdf of the Weibull distribution is defined as follows.

$$f\left({x}_N\right)=\left\{\left(\frac{\beta }{\alpha}\right){\left(\frac{x_N}{\alpha}\right)}^{\beta -1}{e}^{-{\left(\frac{x_N}{\alpha}\right)}^{\beta }}\right\}+\left\{\left(\frac{\beta }{\alpha}\right){\left(\frac{x_N}{\alpha}\right)}^{\beta -1}{e}^{-{\left(\frac{x_N}{\alpha}\right)}^{\beta }}\right\}{I}_N;{I}_N\epsilon \left[{I}_L,{I}_U\right]$$

(5)

where α and β are scale and shape parameters, respectively. Note here that the proposed npdf of the Weibull distribution is the generalization of pdf of the Weibull distribution under classical statistics. The neutrosophic form of the npdf of the Weibull distribution reduces to the Weibull distribution when I_L =0. The neutrosophic cumulative distribution function (ncdf) of the Weibull distribution is given by

$$F\left({x}_N\right)=1-\left\{{e}^{-{\left(\frac{x_N}{\alpha}\right)}^{\beta }}\left(1+{I}_N\right)\right\}+{I}_N;{I}_N\epsilon \left[{I}_L,{I}_U\right]$$

(6)

The neutrosophic mean of the Weibull distribution is given by.

$${\mu}_N=\alpha \Gamma \left(1+1/\beta \right)\left(1+{I}_N\right);{I}_N\epsilon \left[{I}_L,{I}_U\right]$$

(7)

The null and alternative hypotheses for the daily average deaths are stated as follows:

$${H}_0:{\mu}_N={\mu}_{0N}\kern0.5em \mathrm{Vs}.{H}_1:{\mu}_N\ne {\mu}_{0N}.$$

Where μ_N is a true daily average death and μ_0N is the specified daily average deaths. Suppose that t_0N = aμ_0N be the time in days, where a is the termination ratio. The probability of the item will fail before it reaches the experiment time t_0N is defined as follows:

$${\displaystyle \begin{array}{*{20}c}{p}_N=1-\left\{{e}^{-{\left(\frac{t_{0N}}{\alpha}\right)}^{\beta }}\left(1+{I}_N\right)\right\}+{I}_N;{I}_N\epsilon \left[{I}_L,{I}_U\right]\\ {}=1-\left\{\exp \left(-{a}^{\beta }{\left({\mu}_N/{\mu}_{0N}\right)}^{-\beta }{\left(\Gamma \left(1/\beta \right)/\beta \right)}^{\beta }{\left(1+{I}_N\right)}^{\beta}\right)\left(1+{I}_N\right)\right\}+{I}_N\end{array}}$$

(8)

where μ_N/μ_0N is the ratio of true average daily wind speed to specified average daily wind speed. Suppose that $\tilde{\alpha }$ and $\tilde{\beta }$ be type-I and type-II errors. The medical practitioners are interested to apply the proposed plan for testing H₀ : μ_N = μ_0N such that the probability of accepting H₀ : μ_N = μ_0N when it is true should be larger than $1-\tilde{\alpha }$ at μ_N/μ_0N and the probability of accepting H₀ : μ_N = μ_0N when it is wrong should be smaller than $\tilde{\beta }$ at μ_N/μ_0N = 1. In order to find the design parameters n, c₁, c₂ and I_N for the proposed RASP, we consider two points on the OC function. In our approach, the quality level mainly depends on the ratio μ_N/μ_0N. This ratio is helpful for the producer to improve the lot quality. From in producer point of view, the probability of acceptance should be at least $1-\tilde{\alpha }$ at acceptable quality level (AQL), p_1N. So, the producer demands the lot should be accepted at various levels of μ_N/μ_0N. Similarly, from in consumer point of view the lot rejection probability should not be exceeded $\tilde{\beta }$ at limiting quality level (LQL), p_2N. The design parameters are determined by satisfying the following two inequalities

$$$$L\left({p}_{1N}\left|{\mu}_N/{\mu}_{0N}\right.\right)=\frac{\sum \limits_{i=0}^{c_1}\left(\begin{array}{*{20}l} n\\ {}i\end{array}\right){p_{1N}}^i{\left(1-{p}_{1N}\right)}^{n-i}}{\sum \limits_{i=0}^{c_1}\left(\begin{array}{*{20}l} n\\ {}i\end{array}\right){p_{1N}}^i{\left(1-{p}_{1N}\right)}^{n-i}+1-\sum \limits_{i=0}^{c_2}\left(\begin{array}{*{20}l} n\\ {}i\end{array}\right){p_{1N}}^i{\left(1-{p}_{1N}\right)}^{n-i}}\ge 1-\tilde{\alpha}$$

(9)

$$L\left({p}_{2N}\left|{\mu}_N/{\mu}_{0N}=1\right.\right)=\frac{\sum \limits_{i=0}^{c_1}\left(\begin{array}{*{20}l} n\\ {}i\end{array}\right){p_{2N}}^i{\left(1-{p}_{2N}\right)}^{n-i}}{\sum \limits_{i=0}^{c_1}\left(\begin{array}{*{20}l} n\\ {}i\end{array}\right){p_{2N}}^i{\left(1-{p}_{2N}\right)}^{n-i}+1-\sum \limits_{i=0}^{c_2}\left(\begin{array}{*{20}l} n\\ {}i\end{array}\right){p_{2N}}^i{\left(1-{p}_{2N}\right)}^{n-i}}\le \tilde{\beta}$$

(10)

where p_1N and p_2N are defined by

$${p}_{1N}=1-\left\{\exp \left(-{a}^{\beta }{\left(\mu /{\mu}_0\right)}^{-\beta }{\left(\Gamma \left(1/\beta \right)/\beta \right)}^{\beta }{\left(1+{I}_N\right)}^{\beta}\right)\left(1+{I}_N\right)\right\}+{I}_N$$

(11)

$${p}_{2N}=1-\left\{\exp \left(-{a}^{\beta }{\left(\Gamma \left(1/\beta \right)/\beta \right)}^{\beta }{\left(1+{I}_N\right)}^{\beta}\right)\left(1+{I}_N\right)\right\}+{I}_N$$

(12)

The estimated designed parameters of the proposed plan should be minimizing the average sample number (ASN) at an acceptable quality level. The ASN for the proposed plan with fraction defective (p) is derived to be

$$ASN=\frac{n}{P_a(p)+{P}_r(p)}$$

(13)

Therefore, the design parameters for the proposed plan with minimum sample size will be obtained by solving the below optimization technique

$${\displaystyle \begin{array}{*{20}l} \operatorname{Minimize}\ ASN\left({p}_{1N}\right)\\ {}\mathrm{subject}\ \mathrm{to}\ \\ {}L\left({p}_{1N}\right)\ge 1-\tilde{\alpha}\\ {}L\left({p}_{2N}\right)\le \tilde{\beta}\\ {}0\le {c}_1\le {c}_2\\ {}\mathrm{where}\kern0.5em n,{c}_1,{c}_2\in z\end{array}}$$

(14)

The values of the designed parameters n, c₁ and c₂ for various values of $\tilde{\beta }$ =0.25, 0.10, 0.05, 0.01; $\tilde{\alpha }=0.10$; a = 0.5 and 1.0, μ_N/μ_0N =1.1, 1.2, 1.3, 1.4, 1.5, 1.8, 2.0 and I_N =0.0, 0.02, 0.04 and 0.05 when shape parameter β = 1, 2 and 3 are given in Tables 1, 2, 3, 4, 5 and 6. Tables 1 and 2 are shown for the exponential distribution case. For exponential distribution, it can be seen that the values of ASN decrease as the values of a increases from 0.5 to 1.0. On the other hand for other the same parameters, the values of n decreases as the values of β increases. Note here that the indeterminacy parameter I_N also plays a significant role in minimizing the sample size. As indeterminacy parameter I_N increases the ASN values are decreasing.

Table 1 The plan parameter when $\tilde{\alpha }=0.10;\beta =1$ and a = 0.50

Table 2 The plan parameter when $\tilde{\alpha }=0.10;\beta =1$ and a = 1.0

Table 3 The plan parameter when $\tilde{\alpha }=0.10;\beta =2$ and a = 0.5

Table 4 The plan parameter when $\tilde{\alpha }=0.10;\beta =2$ and a = 1.00

Table 5 The plan parameter when $\tilde{\alpha }=0.10;\beta =3$ and a = 0.5

Table 6 The plan parameter when $\tilde{\alpha }=0.10;\beta =3$ and a = 1.0

Results

A comparative study is carried out between the proposed sampling plans with the existing sampling plans available in the literature with respect to the sample size in this section. We know the cost of the study is always directly proportional to the sample size, a plan is said to be economical if it requires a smaller number of samples for testing the hypothesis about the daily new deaths from COVID-19. The proposed repetitive sampling plan under uncertainty/indeterminacy for Weibull distribution is the generalization of the testing average wind speed using sampling plan for Weibull distribution under indeterminacy plan developed by [42]. The comparison for the proposed and the existing sampling plan for Weibull distribution under indeterminacy plan developed by [42] are displayed in Tables 7 and 8 for $\tilde{\alpha }=0.10;\beta =2$ at a = 0.5 and 1.0. The developed sampling plan reduces to the existing sampling plan when c₁ = c₂ = c. From Tables 7 and 8, it is noticed that the values of the sample size required for testing H₀ : μ_N = μ_0N smaller for the proposed sampling plan as compared with the existing sampling plan developed by [42]. For example, when μ_N/μ_0N =1.1 and a =0.5 from Table 7, it can be seen that ASN = 491.58 from the plan proposed sampling plan whereas existing sampling plan sample size n = 617 when I_N =0.02, β = 2 and a = 0.5. Hence, the proposed sampling plan is more economical than the existing sampling plan.

Table 7 Sample size comparison between the proposed plan and existing plan for $\tilde{\alpha }=0.10;\beta =2$ and a = 0.5

Table 8 Sample size comparison between the proposed plan and existing plan for $\tilde{\alpha }=0.10;\beta =2$ and a = 1.0

Discussions

At this juncture, application of the proposed methodology will be illustrated using COVID-19 data belong to Italy of 111 days that are recorded from 1 April to 20 July 2020. The data are available at https://covid19.who.int/. This data is made up of the ratio of daily new deaths (i.e. daily number of deaths over new cases). The data is reported in Table 9. We have taken this data from [46] and they studied applications of COVID-19 data for Kumaraswamy inverted Topp-Leone distribution. Coronavirus disease (COVID-19) is an infectious disease caused by a newly discovered coronavirus. A large number of people affected by the COVID-19 virus and it are infected at random and uncertain, the COVID-19 data follows a certain statistical distribution under neutrosophic statistics. The World health organization and different countries’ health administrators are involved to check the daily affected cases, recovered cases and deaths under indeterminacy. It is found that the COVID-19 data follows the Weibull distribution with shape parameter $\hat{\beta}=$ 2.2222 with the standard error (SE) as 0.1596 and scale parameter $\hat{\alpha}=0.1880$ with SE value as 0.00845. The Kolmogorov-Smirnov test and it p value are D = 0.0684 and p = 0.6766. The goodness of fit of the Weibull distribution is highlight by depicts the histogram and quantile-quantile (Q-Q) plot in Fig. 1. We also applied various life distributions to fit the COVID-19 data set for the intention of comparative study. We have considered here the existing three models like odds Weibull distribution (OWD), Nadarajah-Haghighi distribution (NHD) and Exponentiated Nadarajah-Haghighi distribution (ENHD) for the same data. For more details please refer to [47].

Table 9 COVID-19 data belong to Italy from 1 April to 20 July 2020

Fig. 1

The empirical and theoretical pdfs, empirical and theoretical cdfs, Q-Q plots and p-p plot for the WD for the daily new deaths

Pdf and cdf of Weibull distribution are respectively

$$f(x)=\left(\frac{\beta }{\alpha}\right){\left(\frac{x}{\alpha}\right)}^{\beta -1}{e}^{-{\left(\frac{x}{\alpha}\right)}^{\beta }};x>0,\alpha >0,\beta >0$$

and $F(x)=1-{e}^{-{\left(\frac{x}{\alpha}\right)}^{\beta }}$; x > 0, α > 0, β > 0

Pdf and cdf of odds Weibull distribution (OWD) are respectively (suggested by [48])

$$f(x)=\left(\frac{\alpha \beta}{x}\right){\left(\frac{x}{\theta}\right)}^{\alpha }{e}^{{\left(\frac{x}{\theta}\right)}^{\alpha }}{\left({e}^{{\left(\frac{x}{\theta}\right)}^{\alpha }}-1\right)}^{\beta -1}{\left[1+{\left({e}^{{\left(\frac{x}{\theta}\right)}^{\alpha }}-1\right)}^{\beta}\right]}^{-2};x>0,\alpha <0,0<\left(\beta, \theta \right)$$

and $F(x)=1-{\left[1+{\left({e}^{{\left(\frac{x}{\theta}\right)}^{\alpha }}-1\right)}^{\beta}\right]}^{-1}$; x > 0, α < 0, 0 < (β, θ).

Pdf and cdf of Nadarajah-Haghighi distribution (NHD) are respectively (see [49])

$$f(x)=\left(\alpha \lambda \right){\left(1+\lambda x\right)}^{\alpha -1}{e}^{1-{\left(1+\lambda x\right)}^{\alpha }};x>0,\alpha >0,\lambda >0$$

and $F(x)=1-{e}^{1-{\left(1+\lambda x\right)}^{\alpha }}$; x > 0, α > 0, λ > 0.

Pdf and cdf of Exponentiated Nadarajah-Haghighi distribution (ENHD) are respectively (see [49])

$$f(x)=\left(\alpha \lambda \theta \right){\left(1+\lambda x\right)}^{\alpha -1}{e}^{1-{\left(1+\lambda x\right)}^{\alpha }}{\left(1-{e}^{1-{\left(1+\lambda x\right)}^{\alpha }}\right)}^{\theta -1};x>0,\alpha >0,\lambda >0,\theta >0$$

and $F(x)={\left(1-{e}^{1-{\left(1+\lambda x\right)}^{\alpha }}\right)}^{\theta }$; x > 0, α > 0, λ > 0, θ > 0.

We have estimated the parameters and good fit for the COVID-19 data for WD, OWD, NHD and ENHD, and are reported in Table 10 and depicted in Fig. 2. From Table 10 and Fig. 2 it is noticed that WD shows less AIC, BIC and -2logLL, moreover OWD and NHD are not fitted for COVID-19 data. Hence, Weibull distribution shows a good fit for the COVID-19 data belongs to Italy. The plan parameters for this shape parameter are shown in Tables 11 and 12. For the proposed plan, the shape parameter is ${\hat{\beta}}_N=\left(1+0.04\right)\times 2.2222\approx 2.31$ when I_U =0.04.

Table 10 Estimation and Goodness of fit measures of fitted distribution for daily new deaths

Fig. 2

Comparison of various models for the daily new deaths

Table 11 The plan parameter when $\tilde{\alpha }=0.10;\beta =2.2222$ and a = 0.5

Table 12 The plan parameter when $\tilde{\alpha }=0.10;\beta =2.2222$ and a = 1.0

Suppose that a quality medical practitioner would like to use the proposed repetitive sampling plan for Weibull distribution under indeterminacy to ensure the mean ratio of daily new deaths at least 60 days using the truncated life test for 60 days. Let the producer’s risk be 10% at μ_N/μ_0N =1.1 and the consumer’s risk is 10%. From Table 11, with a = 1.0, $\tilde{\beta }=0.10$ and $\tilde{\alpha }=0.10$ for the repetitive sampling plan, it could be found that the plan parameters are c₁ = 58 c₂ = 66 and ASN = 191.67. Therefore, the plan could be implemented as follows: selecting a random sample of 114 patients from the arrived lot of patients, and doing the truncated life test for 60 days. The proposed sampling plan will be implemented as: accept the null hypothesis H₀ : μ_N = 0.1665 if the average ratio of daily new deaths in 60 days is less than 58, the ratio of daily deaths, but the lot should be rejected as soon as the ratio of daily new deaths exceeds 66. Otherwise, the experiment would be repeated. Table 9 shows the 56 ratios of daily new deaths before the average ratio of daily new deaths of 0.1665. Therefore, the quality medical practitioners would have accepted the arrived lot of patients.

Conclusions

An elaborated investigation of COVID-19 data for Weibull distribution under indeterminacy using time truncated repetitive sampling plan is studied. The proposed design parameters are obtained for known values of the indeterminacy parameters. The plan parameters and corresponding tables are developed for the industrial purposes at specified indeterminacy parametric values. The proposed sampling plan is compared with the existing sampling plans. The result shows that the proposed repetitive sampling plan is more economical than the existing sampling plan. The proposed sampling plan saves time; labor and amount for experimentation, the proposed plan is recommended to apply for testing the average number of deaths due to COVID-19. Also, noticed that if the indeterminacy values increase then the average sample number is decreased. The developed repetitive sampling plan procedure is illustrated with COVID-19 data belong to Italy as an application. The proposed sampling plan can be implemented in various industries covering the packing industry, medical sciences, food industries and electronic industries. Further research can be established to extend our study to group sampling plans, multiple dependent state sampling plans, and multiple dependent state repetitive sampling plans.