 Research
 Open access
 Published:
Application of the multiple dependent state sampling strategy to late adolescent suicide rates
BMC Medical Research Methodology volume 23, Article number: 192 (2023)
Abstract
A multiple dependent state sampling plan (MDSSP) is designed when the lifetimes of the variables follow New Lomax Rayleigh Distribution (NLRD). The decision to accept or reject a lot in the proposed methodology is based on the quality of the given present or previous lots. A binomial modelbased operating characteristic curve (OC curve) for continuous lots of variables under similar settings in healthcare is used in finding the probability of acceptance, acceptance number, rejection number, and the number of preceding (succeeding) lots to consider. Time truncated life test based on the specified median of the NLRD is used in designing the current acceptance sampling plan. For specified values of the parameters of NLRD, quantile ratios, consumer’s risk and producer’s risk, average sampling number (ASN) and probability of acceptance of a lot are reported in tables. Real data on worldwide suicide rates of 15–19 years in the year 2019 from the World Health Organization (WHO) website is considered to illustrate this methodology. The minimum sample size required from the selected data to comment on worldwide suicide rates in late adolescents is explained with MDSSP. The results of the proposed acceptance sampling method are compared with the singlestage sampling plan.
Introduction
Quality in the manufacturing or the service sector is a determinant for attracting demand. Continuous inspection and improving quality in the healthcare sector gained significance for the benefit of the patients and service providers. In this sector, quality is in two ways healthrelated and hospitality related. The former is objective, and the latter is subjective, it is equally important to measure both in the healthcare setting for healthy clinical outcomes and the satisfaction of the customer. Measuring, monitoring and controlling are parts of health policy design. Many health indicators are the targets of sustainable development goals (SDGs). Reducing deaths due to suicides is one of the priorities of the WHO mental health action plan 2013–2030 and is also goal 3.4 of SDGs. This paper aims to introduce the use of MDSSP in estimating the required sample size in monitoring and commenting on worldwide health indicators.
Sampling plans specify parameters of the sampling process like the number of samples to inspect from a selected lot, average Sampling Number (ASN), acceptance criteria: acceptable quality level (AQL), limiting quality level (LQL), producer’s (service providers) risk and consumer’s risk. Acceptance sampling plan (ASP), a tool of statistical quality control (SQC) spread its application wings to the manufacturing and service industry after its first use in World War II. In the service sector, inspect a lot with minimum cost, communicate with management to necessitate action. Acceptance sampling plans operating procedure is designed differently in lotbylot inspections and continuous flow processes (services). In the first case, single sampling and double sampling plans are useful, for the second type twostage sampling, chain sampling, and multiple state dependent sampling plans are suitable.
The most common phenomenon in life testing is truncating a life test after a predefined time \(t_{0}\) and noting the number of non confirming members/items in the inspecting lot before the pre specified time.Acceptance sampling plans literature for different symmetric and skewed distributions under truncated life testing is structured and coherent. The first single acceptance sampling plan based on exponential lifetimes was introduced by [1]. Truncated life test sampling plans based on normal and lognormal distributions are designed by [2]. Average life based sampling plans for loglogistic distribution are given by [3]. Rayleigh distribution based time truncated acceptance sampling plans were discussed by [4]. ASP for BirnbaumSaunders distribution based on median life is illustrated by [5]. Percentile based acceptance sampling plans are given by [6, 7] for BirnbaumSaunders and Burr typeXII distribution. AlNasser et al. [8] discussed Power Lomax acceptance sampling based on truncated life testing. Group acceptance sampling plans for resubmitted lots, sizebiased Lomax distribution and MDSSP for exponentiated loglogistic distribution are given by [9, 10] respectively. Acceptance sampling plans for new WeibullPareto distribution based on percentiles are given by [11].
Application of statistical quality control tools in health care
Unlike the manufacturing industry service sector suffers from more variability in outcome. Especially in health care service, input variability is more due to inevitable biological variability among human beings. So it is a great challenge to the healthcare sector to deliver quality clinical outcomes. Health care sector started employing SQC tools way back in the early nineteenth century mentioned in [12] or additional information on the SQC techniques usefulness in managing quality indicators in hospitals refer to [12,13,14]. Callahan et al. [15] used control charts to analyze patient waiting times and time from registration to the physician’s first orders. Rachmania et al. [16] discussed the application of a lot acceptance sampling plan in screening women and concluded that this procedure is useful in deciding the quality level of the indicator. Clemente et al. [17] mentioned the use of statistical process control in the emergency medicine department studying doortoreperfusion time of myocardial infarction patients. Ray et al. [18] presented a case study in Indonesia on the use of control charts in hospitals. Rao et al. [19] discussed the use of control charts usefulness in monitoring surgical site infections and comparing mortality rates in CABG surgery and other surgeries. Quality indicators and benchmarks for different ICUs are given by [20]. The use of a multiplestate dependent sampling plan in assessing COVID19 mortality rate lots is discussed by [21]. Statistical process control use in economic analysis and safety of chemotherapy batches in hospital pharmacies is mentioned by [22]. After this literature review, it is observed that distributionbased control charts and acceptance sampling plans were not of much use in healthcare quality analysis. With this motivation, we developed the MDSSP based on NLRD and its application in analyzing quality gaps and acceptability of worldwide suicide rates of late adolescents in 2019.
Multiple dependant state sampling plans (MDSSP)
Conditional sampling plans are suitable for continuous process monitoring and improvement; it is first introduced by [23]. Chain sampling and multiple dependent state sampling come under continuous sampling plans. Different variations in chain sampling plans were studied by [24]. In 1976, [25] introduce another sampling plan, named MDSSP belonging to the conditional sampling family. MDSSP for measurement data was discussed by [26]. Aslam et al. [27] designed novel multiple dependent state sampling plan based on process loss consideration. The Multiple dependent state repetitive group sampling plan for normally distributed quality characteristics is given by [28]. Rao et al. [29] gave novel control charts based on multiple dependent state sampling plans. Some recent works to mention on MDSSP are [9, 10, 30] developed MDSSP based on repetitive group sampling for exponentiated halflogistic distribution and exponentiated Weibull distribution respectively. Modified multiple dependent state sampling plan is introduced and its economic efficiency is studied by [31]. The MDSSP efficiency over other sampling plans in reducing ASN is given by [32].
In the following sections, medianbased probability of failure formulas for NLRD distribution and MDSSP when quality (health) indicator variable follows NLRD is derived.
New Lomax Rayleigh distribution
NLRD is a new TX family distribution proposed by [33] considering Lomax distribution, a special case of Pareto TypeII and Rayleigh distribution, its suitability in estimating life time data was verified. Different generalizations of Lomax distribution are available in the literature to make it suitable for studying lifetime data. To estimate drilling machine lifetime, [34] used the LomaxG generator to derive LomaxRayleigh distribution and applied it to drilling machine data. Different types of estimators of ParetoRayleigh distribution are given by [35]. Different traditional and heuristic methods of estimation for the LomaxRayleigh distribution generated by the LomaxX generator are given by [36]. Rady et al. [37] used Power Lomax distribution and tested its applicability in estimating the mean remission times of bladder cancer patients. Bladder cancer remission times data is studied with Pareto Weibull distribution by [38].
Let T be a quality indicator following NLRD with parameters \(\theta ,\lambda\) and \(\sigma\) then the cumulative distribution function (CDF) and probability density function (PDF) of the variable T are given as follows:
where \(\lambda\) and \(\sigma\) are the scale parameters and \(\theta\) is the shape parameter.
The \(q^{th}\) quantile of NLRD is
Let \(\eta_{q} = \left[ {2\lambda \left[ {\left[ {1  q} \right]^{{\frac{  1}{\theta }}}  1} \right]} \right]^{\frac{1}{2}}\), then \(t_{q} = \sigma \eta_{q} \Rightarrow \sigma = \frac{{t_{q} }}{{\eta_{q} }}\).
For specified values of \(\lambda = \lambda_{0} ,\theta = \theta_{0}\), \(t_{q}\) is a function of the scale parameter \(\sigma = \sigma_{0}\).
Therefore \(\sigma_{0} = \frac{{t_{q} }}{{\left[ {2\lambda \left[ {\left[ {1  q} \right]^{{\frac{  1}{\theta }}}  1} \right]} \right]^{\frac{1}{2}} }}\).
Median of the NLRD distribution, 50^{th} percentile is
Let p be the probability of failure of the variable following NLRD with terminating time \(t_{0}\) and truncating time \(t_{q}\), then \(p = F(t_{0} )\). Experiment termination time can be expressed as \(t_{0} = kt_{{q_{0} }}\).
By substituting \(t_{0} = kt_{{q_{0} }}\) and \(\sigma_{0}\) values in Eq. (5)
Here it is important to consider both consumer’s risk and the producer’s risk. Quantile ratio plays a significant role in decision making about lot acceptance. As per Producer’s angle of lot acceptance quantile ratio must be at least one. If \({{(t_{q} } \mathord{\left/ {\vphantom {{(t_{q} } {t_{{q_{0} }} }}} \right. \kern0pt} {t_{{q_{0} }} }}) > 1\), then the probability (\(p_{1}\)) obtained is considered as AQL and the consumer wants this ratio must be almost the pre defined level of the consumer’s risk \(\beta\). LQL is the probability of this ratio \({{(t_{q} } \mathord{\left/ {\vphantom {{(t_{q} } {t_{{q_{0} }} }}} \right. \kern0pt} {t_{{q_{0} }} }}) = 1\) being one(\(p_{2}\)). NLRD parameters can be obtained from previous data on mortality rates. MDSSP is designed by substituting estimated parametric values of NLRD at different quantile ratios, terminating ratios, consumer’s risks and producer’s risk; the average sample number required to sentence a lot of world suicide rates is estimated with the WHO fact sheet data on suicides in 2019.
Operating procedure for NLRD based MDSSP
The operating procedure and designing methodology of MDSSP based on NLRD are described in detail in the subsequent subsections.
Operating procedure
MDSSP consists of the following symbols and notations:
N Lot size, n sample size,
\(P_{a} (p)\)Probability of acceptance for a given p,
\(p_{1}\)Probability of rejecting lot that has acceptable quality level (AQL), related to producers risk.
\(p_{2}\)Probability of acceptance of lot with limiting quality level (LQL), associated with consumer’s risk.
\(\alpha\) Producers risk, \(\beta\)Consumers risk.
ASN Average Sample Number
\(c_{1} ,c_{2}\) Unconditional and conditional acceptance numbers respectively.
m Required number of past or future lots to accept the present lot.
d Number of units failed before the terminating time \(t_{0}\) in the sample selected from the lot under test.
Step1 From the lot, select a sample of units with size n, and put all these n items under test till predefined time \(t_{0}\).
Step2 Observe d, the number of units that failed before the test terminating time \(t_{0}\).
Step3Compare d with unconditional and conditional acceptance number \(c_{1} ,c_{2}\) and make a decision.

a.
If \(d \le c_{1}\) accept the lot otherwise reject the lot.

b.
If \(c_{1} < d \le c_{2}\), then consider the number of failures in the preceding lots before terminating \(t_{0}\) and this number must be less than \(c_{1}\) in all the m previous lots.
The operating characteristic (OC) function of MDSSP is given by
The probabilities in Eq. (7) are obtained by using the binominal theorem and the decision rule designed about lot acceptance at p is as follows:
Designing methodology
The objective of designing any sampling plan is to avoid laborious 100% inspection and achieve ASN that minimizes the resources of the organizations. MDSSP also reduces the inspection time and cost by deciding the minimum number of units required to verify whether the sampled lot of quality indicator is acceptable or not. If the lot is rejected communication should be given to the management about the process (service). In this section, we tried to achieve minimum ASN through an optimization problem. The proposed optimization procedure for MDSSP based on NLRD is given below:
Minimize \(ASN(p) = n\),
Subject to
\(P_{a} \left( {p_{2} } \right) \le \beta\),
The probability of AQL and LQL is obtained from the following equations:
In this paper median quantile ratio \(t_{q} /t_{{q_{0} }}\) at the consumer’s risk must be at least 1, \(t_{q} /t_{{q_{0} }} = 2,4,6,8,10\) are considered at the producer’s risk. The optimal parameters of the proposed plan for NLRD \(\left( {\theta ,\lambda } \right) = \left( {1.5,1.5} \right),\left( {0.5,0.5} \right),\left( {1.5,1.0} \right),\left( {1.5,2.0} \right)\) are presented in the Tables 1, 2, 3 and 4. Values assumed for consumer’s risk are \(\beta = 0.25,0.10,0.05,0.01\) and producer’s risk is considered as \(\alpha = 0.05\) at 50^{th} percentile. The values considered for the termination ratio are \(k = 0.5,0.7,1.0\). The computational work is done by using R software and R code is provided in Supplementary file.
Results observed from Tables 1, 2, 3 and 4 when parametric combinations are fixed are as follows.

1.
An inverse relationship is observed between sample size and consumer’s risk. In all the Tables 1, 2, 3 and 4 sample size increases when the consumer’s risk decreases.

2.
Sample size decreases as the termination ratio \(k\) increases from 0.5 to 1.0

3.
Probability of lot acceptance increases along with the quantile ratio. As \(t_{q} /t_{{q_{0} }}\) approaches ‘10’ probability of lot acceptance also increases and approximates to almost ‘1’.
Applications of MDSSP for real data
In this subsection, we demonstrate the proposed methodology with two real data sets one related to public health issue and engineering application related to failure time of the components.
Worldwide suicide rates data
Suicide is a preventable death. As it is a global public health issue impacting all age groups, suicide prevention is a high priority condition in the WHO Mental Health Action Plan 2013–2030. Losing an adolescent in suicide is a great loss to the family as well as the country. Reasons for suicide in middle and late adolescence may be due to not addressing mental health issues. The adolescent state is more vulnerable to mental health problems due to psychophysiological changes, environment and circumstances. Late adolescence in India has become stressful due to their senior secondary school education under heavy pressure towards productivity and achievement without any physical activity. According to WHO Facts on Suicide 2019, suicide is the fourth leading cause of death among late adolescents. The major contribution (77%) of global suicidal deaths is from low and middle income countries (data available in Supplementary Material file).
As per the WHO fact sheet on 15–19 years suicides in the year 2019, overall suicide rates are high in Western specific, Europe and American countries (Guyana (40.37 per 100,000)) in India it is reported as 10.4 per 100,000. In females highest percentage of suicides is reported in the American country Guyana (45.71 per 100,000), Indian female late adolescents suicide rate is 14.55 per 100,000 which is double that of males (6.71 per 100,000). The highest suicide rate in males was reported in the Western Pacific country, Kiribati (51.93 per 100,000). Worldwide, it is the need of the hour to design effective suicide prevention strategies to reduce its contribution to premature deaths due to non communicable diseases. From the above discussion it is observed that, world suicide rates of late adolescents’ data are highly skewed, the applicability of NLRD distribution proposed by [33] to heavily rightskewed data is verified and compared with some other existing distributions.
In this section worldwide suicide mortality rates downloaded from WHO factsheet data on suicides in 2019 in the age group 15–19 are used to estimate the parameters of the NLRD. The goodness of fit of the NLRD distribution to this data is tested by using the Kolmogorov–Smirnov test, the value of D is 0.0291 and the p value is 0.7908. Estimated values of the parameters of the NLRD are obtained by using the maximum likelihood estimation method and they are \(\widehat{\theta } = 1.0559\), \(\widehat{\lambda } = 0.6867\) and \(\hat{\sigma } = 3.7817\). Figure 1 displayed the visual presentation of fitted model, these graphs explain to us about the goodness of fit of the NLRD distribution to worldwide suicide rates of all the countries.
For instance, the industrialist would like to use the developed multiple dependent state sampling plans to implement the median life percentile of the product where the product lifetime follows NLRD with the shape parameters are \(\widehat{\theta } = 1.0559\) and \(\widehat{\lambda } = 0.6867\). By 2030, suicide prevention and reducing the suicide mortality rate by one third is Goal 3 of Sustainable Development Goals. It may be necessary to check whether the suicide death rates are reduced year by year or not. The medical practitioners suggest that the given median suicide rate is 0.95 whereas the medical practitioners expected that the median suicide rate is 1.9. The consumer’s risk is 0.05 if the actual median suicide rate is 0.95 and the producer’s risk is 0.10 if the actual median suicide rate is 1.9. With these constraints, the optimal parameters selected from Table 5 are n = 37, \(c_{1}\) = 3, \(c_{2}\) = 9, and m = 2 with values of \(\widehat{\theta } = 1.0559\), \(\widehat{\lambda } = 0.6867\),\(t_{{q_{0} }}\) = 0.75, \(\alpha\) = 0.05, \(\beta\) = 0.05, \(t_{q} /t_{{q_{0} }}\) = 2 at k = 0.5. The MDSSP are illustrated as follows:
In this context, a lot is a group of suicide rates of some of the countries selected at random from worldwide suicide rates. For example, a sample of 37 countries suicide rates of late adolescents in the age group of 15–19 years will be selected at random for the group of young people and check their suicide rate is 0.95. If the suicide rate before 0.95 is in 3 countries in a group of the population will be accepted and the group of the population will be rejected if it is greater than 9 countries in a group. There will be indecision of the group of the population is deferred until the 2 preceding the group of the population will be tested in case of the suicide rate of countries in younger people in a group of the population is between 3 and 9. For this real example, there are 16 countries’ suicide rates in younger people in a group of the population before the suicide rate before 0.75. Hence, reject the suicide rates of the countries in late adolescent people in a group of the population. Thus medical practitioners could suggest to the government or public that the median suicide rate of countries in late adolescent people in a group of the population is at an unacceptable level.
Performance of NLRD is compared with five other related distributions; Rayleigh distribution (RD), Lomax distribution (LD), Pareto Type–II(PD2), Power Lomax distribution (PLD) and Lomax Rayleigh distribution (LRD). Estimates of 2log likelihood, AIC, BIC, CAIC and HQIC of NLRD compared with other above mentioned distributions is given below. A better fit of the distributions to real data sets is estimated with CramerVon Mises statistic (W^{*}), Anderson–Darling statistic (A^{*}) and KolmogrovSmirnov statistic (KS).The distribution with lower A,W, and KS values is a better fit. Maximum likelihood estimates of parameters of the NLRD and KS test pvalues are presented in Table 6 and AIC, BIC, CAIC, HQIC, W^{*} and A^{*} values are reported in Table 7. From Tables 6 and 7 it is evident that, the intended model is a satisfactory model among several models. The Performance of NLRD is compared with five other related distributions RD, LD, PD2, PLD and LRD. In both real data sets, our anticipated model fits perfectly when compared to all other distributions. As a result, the anticipated model is an alternative model for real data modeling in lifetime data, specifically in the medical field.
Failure time data
We consider the real time application to illustrate the suggested MDS sampling plan for engineering industrial use, if lifetime of an item comes from NLRD with unknown shape parameters. The following data is acquired from [39] and this data set corresponds to time between failures of 30 repairable components. Failure time (days): 0.11, 0.30, 0.40, 0.45, 0.59, 0.63, 0.70, 0.71, 0.74, 0.77, 0.94, 1.06, 1.17, 1.23, 1.23, 1.24, 1.43, 1.46, 1.49, 1.74, 1.82, 1.86, 1.97, 2.23, 2.37, 2.46, 2.63, 3.46, 4.36, 4.73.
The demonstration of the goodness of fit for the given model is shown in Fig. 2, the empirical and theoretical cdfs and QQ plots for the NLRD for the time between failures of repairable items and also in Fig. 2, Estimates of the density functions for the time between failures of repairable items. The maximum likelihood estimation of the parameters of NLRD for the time between failures of 30 repairable items is are \(\hat{\theta }\) = 1.9273 and \(\hat{\lambda }\) = 1.8958 and the maximum distance between the real time data and the fitted of NLRD was found from the Kolmogorov–Smirnov test is 0.0862 and also the pvalue is 0.979.
Assume that the investigator wish for implement the median life of the product for the proposed MDS sampling scheme when the product lifetime comes from NLRD with estimated shape parameter is \(\hat{\theta }\) = 1.9273 and \(\hat{\lambda }\) = 1.8958. Let the specified quantile time between failures of repairable items is 0.5 i.e. \(t_{q}^{0} = 0.65\) and the experiment termination time is 0.65, i. e. \(t_{0}\) = 0.65. Hence, the termination constant, \(k = {{t_{0} } \mathord{\left/ {\vphantom {{t_{0} } {t_{q}^{0} }}} \right. \kern0pt} {t_{q}^{0} }} = 1.0\). Assuming that,\(\alpha = 0.05\), \(\beta = 0.10\),\(k = 1.0\) and \({{t_{q} } \mathord{\left/ {\vphantom {{t_{q} } {t_{q}^{0} }}} \right. \kern0pt} {t_{q}^{0} }}\) = 2, so from Table 8 the design parameters are n = 30, \(c_{1}\) = 2, \(c_{2}\) = 4 and m = 1. Thus, the design could be carried out as persists: obtain 30 samples at random from the current lot. Accept the lot if two failures occur before 0.65 time between failures of repairable items. If more than 4 failures the lot should be rejected. Whereas, if the number of failures are between 2 and 4, the outlook of the lot will be dependent until the preceding one lot is tested. If the preceding one lot is accepted, then product could be accepted. Otherwise, reject present lot. For the present example six failures are noted before the time between failures of repairable items 0.65. Therefore, the lot is rejected.
Comparison of proposed MDSSP with the singlestage sampling plan
The comparison is made for single and MDSSP when quality control follows NLRD, the OC curve is used here to show the efficiency of the plan. The curve has displayed the difference in probabilities of accepting a good lot as well as rejecting a bad lot. Table 9, has revealed the efficiency of the proposed MDSSP over the single sampling plan (SSP) while assuming the underlying distribution of data to follow NLRD. Considering the quantile ratio \({{t_{q} } \mathord{\left/ {\vphantom {{t_{q} } {t_{{q_{0} }} }}} \right. \kern0pt} {t_{{q_{0} }} }} = 2,4,6,8,10\) for each consumer's risk \(\beta = 0.25,0.10,0.05,0.01\) while keeping the producer’s risk at \(\alpha = 0.05\). The comparison is basically on the sample size \(n\) and probability of acceptance \(P_{a} \left( {p_{1} } \right)\). The acceptance sample size for the proposed MDSSP is smaller than the existing single sampling plan for several set parameters see Table 9. For quantile ratio 2, the plan parameters for the MDSSP are n = 21, \(c_{1}\) = 2, \(c_{2}\) = 4, m = 2 whereas for SSP the design parameters are n = 33, c_{1}= 4 with a corresponding probability of acceptances are 0.9792 and 0.9720 respectively when \(\beta = 0.25\), r = 2, \(\theta = 1.5, \, \lambda { = 1}{\text{.5}}\). The sample size is smaller for the MDSSP as compared with SSP. As the quantile ratio increased the acceptance sample size decreased for both sampling plans. Figure 3 depicts the OC curve for comparison of MDSSP with plan parameters n = 40, \(c_{1}\) = 3, \(c_{2}\) = 13, m = 2 and SSP with n = 61, \(c_{1}\) = 6 when \(\beta = 0.05\). It is noticed that the MDSSP is convincingly greater efficient than SSP in terms of sample size.
Conclusions
Multiple dependent state sampling plans are developed when the lifetimes follow the New Lomax Rayleigh distribution. Lifetimes are truncated at a specified time point. At predefined producer’s and consumer’s risk optimal parameters of the proposed sampling plans are estimated. As suicide rates are continuously observed over years, and reducing suicide deaths is one of the SDGs, application of the proposed MDSSP to check acceptability of the countries’ suicide rates worldwide in late adolescents is verified when the data follows NLRD. We conclude that the proposed MDSSP is more efficient than the existing plans and it can be applied to verify the acceptability of different variables in health statistics. In future applications, multiple dependent state sampling plans may be applied in understanding health indicators quality, hospital equipment manufacturing industries can use this plan to sentence on manufactured lot quality (a lot of health quality indicators) and the average sampling number needs to decide on the quality of the lot. Some examples where it can be applied to evaluate the quality in the utilization of operation theatre timings, waiting times in outpatient departments, waiting times for lab reports and quality care in lifesaving treatments like the door to balloon times in cardiology, etc. The proposed MDSSP helps to understand the quality of services provided today not only based on today’s data but also from previous days’ lots to ensure consistency in quality health care. The future scope of the work is to develop NLRD based attribute control charts, control charts based on multiple dependent state sampling, multiple dependent state repetitive sampling and modified multiple dependent state sampling plan for NLRD distributed quality indicators.
Availability of data and materials
Data is available in Supplementary Material file. Source of the data link is also provided. The datasets generated and/or analysed during the current study are available in the [WHO website] repository, [who.int/data/gho/data/themes/mentalhealth/suiciderates].
Abbreviations
 ASN:

Average sampling number
 AQL:

Acceptable quality level
 CDF:

Cumulative distribution function
 PDF:

Probability density function
 LQL:

Limiting quality level
 ASP:

Acceptance sampling plan
 MDSSP:

Multiple dependent state sampling
 NLRD:

New Lomax Rayleigh Distribution
 OC curve:

Operating characteristic curve
 SDGs:

Sustainable development goals
 SQC:

Statistical quality control
 WHO:

World health organization
References
Epstein B. Truncated Life tests in the exponential case. Ann Math Stat. 1954;19(3):1403–33.
Gupta SS. Life test sampling plans for normal and lognormal distributions. Technometrics. 1962;4(2):151–75. https://doi.org/10.1080/00401706.1962.10490002.
Kantam RRL, Rosaiah K, Rao GS. Acceptance sampling based on life tests: loglogistic model. J Appl Stat. 2001;28(1):121–8.
Baklizi A, ElMasri AQ, ALNasser A. Acceptance sampling plans in the rayleigh model. Commun Stat Appl Methods. 2005;12(1):11–8.
Balakrishnan N, Leiva V, López J. Acceptance sampling plans from truncated life tests based on the generalized BirnbaumSaunders distribution. Commun Stat Simul Comput. 2007;36(3):643–56.
Lio YL, Tsai TR, Wu SJ. Acceptance sampling plans from truncated life tests based on the birnbaumsaunders distribution for percentiles. Commun Stat Simul Comput. 2009;39(1):119–36.
Lio YL, Tsai TR, Wu SJ. Acceptance sampling plans from truncated life tests based on the Burr type XII percentiles. J Chinese Inst Ind Eng. 2010;27(4):270–80.
AlNasser AD, ulHaq MA. Acceptance sampling plans from a truncated life test based on the power Lomax distribution with application to manufacturing. Stat Transit. 2021;22(3):1–13.
Rao GS, Rosaiah K, Rameshnaidu C. Design of multipledeferred state sampling plans for exponentiated half logistic distribution. Cogent Math Stat. 2020;7:1857915. https://doi.org/10.1080/25742558.2020.1857915.
Rao GS, Jilani SD, Rao AV. Designing of multiple dependent state repetitive sampling plan for typeII generalized half logistic distribution. Int J Syst Assur Eng Manag. 2021;12(6):1397–409.
Shrahili M, AlOmari AI, Alotaibi N. Acceptance sampling plans from life tests based on percentiles of new weibull–pareto distribution with application to breaking stress of carbon fibers data. Processes. 2021;9(11):1–11.
Benneyan JC. Statistical quality control methods in infection control and hospital epidemiology, part II: chart use, statistical properties, and research issues. Infect Control Hosp Epidemiol. 1998;19(4):265–83.
Finison LJ, Finison KS. Applying control charts to quality improvement. J Healthc Qual. 1996;18(6):32–41.
Jané AC, Cintas PG. Lot sampling plans in the measure of quality of care indicators. Int J Qual Heal Care. 1999;11(2):139–45.
Callahan CD, Griffen DL. Advanced statistics: Applying statistical process control techniques to emergency medicine: a primer for providers. Acad Emerg Med. 2003;10(8):883–90.
Rachmania I, Setyaningsih S, Rakhmaniar M, Basri M. Application of quality tools and techniques in hospital: case study in Bandung Indonesia. 2012. p. 1–12.
Clemente F, Papi M, Pontecorvi L, Menichetti A. Evaluation of indices for the measurement of quality in health systems. Int J Metrol Qual Eng. 2016;7(4):1–14.
Ray B, Samaddar DP, Todi SK, Ramakrishnan N, John G, Ramasubban S. Quality indicators for ICU: ISCCM guidelines for ICUs in India. Indian journal of critical care medicine : peerreviewed, official publication of Indian society of critical care medicine. Indian J Crit Care Med. 2009;13(4):173–206.
Rao GS, Aslam M. Inspection plan for COVID19 patients for Weibull distribution using repetitive sampling under indeterminacy. BMC Med Res Methodol. 2021;21(1):1–15.
Paci A, Borget I, Mercier L, Azar Y, Desmaris RP, Bourget P. Safety and quality assurance of chemotherapeutic preparations in a hospital production unit: acceptance sampling plan and economic impact. J Oncol Pharm Pract. 2012;18(2):163–70.
Wortham AW, Mogg JM. Dependent stage sampling inspection. Int J Prod Res. 1970;8(4):385–95. https://doi.org/10.1080/00207547008929855.
Stephens KS, Dodge HF. Comparison of chain sampling plans with single and double sampling plans. J Qual Technol. 1976;8(1):24–33.
Wortham AW, Baker RC. Multiple deferred state sampling inspection. Int J Prod Res. 1976;14(6):719–31.
Balamurali S, Jun CH. Multiple dependent state sampling plans for lot acceptance based on measurement data. Eur J Oper Res. 2007;180(3):1221–30.
Aslam M, Yen CH, Chang CH, Jun CH. Multiple dependent state variable sampling plans with process loss consideration. Int J Adv Manuf Technol. 2014;71(5–8):1337–43.
Balamurali S, Aslam M. Determination of multiple dependent state repetitive group sampling plan based on the process capability index. Seq Anal. 2019;38(3):385–99. https://doi.org/10.1080/07474946.2019.1648930.
Aslam M, Nazir A, Jun CH. A new attribute control chart using multiple dependent state sampling. Trans Inst Meas Control. 2015;37(4):569–76.
Jeyadurga P, Balamurali S. Optimal designing of multiple deferred (dependent) state repetitive group sampling plan for variables inspection. Commun Stat  Theory Methods. 2020;51(3):1–17. https://doi.org/10.1080/03610926.2020.1814815.
Rao GS, Rosaiah K, Naidu CR. Design of multipledeferred state sampling plans for exponentiated half logistic distribution. Cogent Math Stat. 2020;7(1):1857915. https://doi.org/10.1080/25742558.2020.1857915.
Rao GS, Fulment AK, Peter JK. Design of multiple dependent state sampling plan application for COVID19 data using exponentiated weibull distribution. Complexity. 2021;2021:1–10.
Aslam M, Jeyadurga P, Balamurali S, Azam M, Almarshadi A. Economic determination of modified multiple dependent state sampling plan under some lifetime distributions. J Math. 2021;2021:1–13. https://doi.org/10.1155/2021/7470196.
Yen CH, Chang CH, Lee CC. A new multiple dependent state sampling plan based on onesided process capability indices. Int J Adv Manuf Technol. 2023;2023:1–13.
Nagasaritha K, Rao GS, Rosaiah K. Survival analysis of cancer patients using a new Lomax Rayleigh distribution. J Appl Math Stat Informatics. 2023;19(1):19–46.
Venegas O, Iriarte YA, Astorga JM, Gómez HW. LomaxRayleigh distribution with an application. Appl Math Inf Sci. 2019;13(5):741–8.
Jebeli M, Deiri E. Estimation methods for the probability density function and the cumulative distribution function of the ParetoRayleigh distribution. Statistics (Ber). 2020;54(1):135–51.
AlAnber NJ. LomaxRayleigh distribution: traditional and heuristic methods of estimation. J Phys Conf Ser. 2020;1591(1):1–12.
Rady EHA, Hassanein WA, Elhaddad TA. The power Lomax distribution with an application to bladder cancer data. Springerplus. 2016;5(1):1–22.
Rana MS, Shahbaz SH, Shahbaz MQ, Rahman MM. Paretoweibull distribution with properties and applications: a member of ParetoX family. Pakistan J Stat Oper Res. 2022;18(1):121–32.
Hassan AS, Elgarhy M, Ahmad Z. Type II generalized ToppLeone family of distributions: properties and applications. J Data Sci. 2019;17(4):638–59.
Acknowledgements
The authors are deeply thankful to the editor and reviewers for their valuable suggestions to improve the quality of the paper.
Funding
The author(s) declared that no grants were involved in supporting this work.
Author information
Authors and Affiliations
Contributions
Kolli write a manuscript, Rosaiah supervised the work and methodology, Rao developed computational work and methodology, Kirigiti constructed figures and tables.
Author’s information
None.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This study was based on published data, so ethical approval was not required for published data.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Kolli, N., Rosaiah, K., Rao, G.S. et al. Application of the multiple dependent state sampling strategy to late adolescent suicide rates. BMC Med Res Methodol 23, 192 (2023). https://doi.org/10.1186/s12874023020072
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12874023020072