key: cord-1050196-ck9zlts7 authors: El‐Hadidy, Mohamed Abd Allah; Alfreedi, Ajab A. title: Detection of an appropriate pharmaceutical company to get a suitable vaccine against COVID‐19 with minimum cost under the quality control process date: 2021-04-07 journal: Qual Reliab Eng Int DOI: 10.1002/qre.2881 sha: 82154bce0481a9ebbdc37f72db0b3f93e644508c doc_id: 1050196 cord_uid: ck9zlts7 Some international pharmaceutical companies have succeeded in producing vaccines against COVID‐19. Countries all over the world have aimed to obtain these vaccines with minimum cost. We consider a set of K‐independent Markovian waiting lists. Each list contains a set of countries, where each one of them has an exponential service time and a Poisson arrival process. These companies differ in some characteristics such as the vaccine production cost and the speed of the required quantity delivery. We present a new detection model that helps in providing an appropriate decision to choose a suitable company. Moreover, the concept of balking and the retention of reneged countries is taken into consideration under the quality control process of each waiting list. Under steady state, we face an interesting and difficult discrete stochastic optimization problem. Its solution gives an optimal distribution of the searching effort, which is bounded by a known probability distribution. A simulation study has been derived to get the minimum value of the paid cost random values. The highest service rate, the total expected profit of each queuing system, and the optimum performance measures, which depend on this cost, have been obtained to show the effectiveness of this model. but more than half of them complained of mild or moderate effects such as fatigue, headache, chills, muscle aches, or pain at the injection site. Also, AstraZeneca plc is considered one of the most prominent companies in this field. The vaccine, which was developed by the Oxford University with AstraZeneca plc, helped more than a thousand people to develop antibodies and white blood cells that can fight against COVID-19. Many other companies have also entered in this competition, such as Johnson & Johnson, the American Gilead company, the British Glasgow Smith company, and other international companies and institutions. This virus's danger on the human life has led to a great competition among the countries of the world to quickly obtain a vaccine and with minimum cost. This prompted us to present an efficient model that helps the country in taking an appropriate decision to choose a suitable company among these international companies and institutions. Therefore, the optimal search theory for the lost targets is more effective here to take this decision. This theory has a remarkable importance in our real life where it is used to detect the lost targets (randomly located or moved) with the maximum probability and minimum cost. In an earlier work, many types of search plans have been derived on these line by Beck et al., [1] [2] [3] El-Hadidy et al., [4] [5] [6] [7] [8] [9] [10] [11] [12] Balkhi, 13 and Reyniers. 14, 15 Many studies have presented different search techniques on the plan and space, such as El-Hadidy et al.. [16] [17] [18] [19] [20] [21] [22] [23] [40] [41] For more recent and interesting search plans, see El-Hadidy et al.. [24] [25] [26] All these studies have aimed to detect the lost targets with minimum cost and maximum probability. Another aim is to study the finiteness of the first meeting time between the randomly moving target and the searchers as in El-Hadidy et al.. [27] [28] [29] In this work, we present an effective probabilistic model to overcome the waiting problem to get a vaccine where time is an important factor in containing the spread of the virus, especially in the winter season. This situation is similar to the queuing problem, where clients (countries) seek to obtain a good service at a suitable time with minimum cost. Due to the large number of countries which wish to obtain this vaccine, the country may be balking and retention of reneging the waiting list (queue). Waiting for a long time in a list (queue) may lead to a situation when a country may became impatient. Consequently, it leaves the waiting list before receiving the required service (this state is called a reneged case; see Shortle et al. 39 ). In this model, we can assist countries in making the appropriate decision to choose the appropriate pharmaceutical company in terms of the required cost and time to obtain the vaccine. This choice also depends on the financial capabilities of each country. There exists another benefit of this model helps countries to take an appropriate economic decision in any economic field. Here, we use the same model which was recently used in El-Hadidy, 30 but with quality control for each queue through the presence of an inspector. If any defect happened, then this inspector will return the country to the waiting list again. This puts us in the face of a difficult discrete stochastic optimization problem that was previously studied in El-Hadidy. 25, 26 The main goal here is to obtain the optimal distribution of the searching effort to get the suitable pharmaceutical company with the maximum probability. In the case of multiple M/M/1 cooperative queues where each one of them contains a single server and one access flow, Anily and Haviv 31 provided a collaboration model between them. Also, Timmer and Scheinhardt 32.33 studied another example of improving the service levels with a good number of clients. In addition, the study, which is based on the controlling process of the client access, is derived from Anily and Haviv 34 where the total access flow was divided among all the servers in the system. The ability to manage the available resources in an ideal way is one of the most important reasons for the institution's success. Therefore, the right decision making of managing the waiting list problem to obtain the vaccine is one of the factors to overcome the spread of COVID-19. It also helps the country to succeed in managing its institution where this decision takes into account the minimum waiting time and cost to obtain this vaccine. Thus, the main objective here is to determine the appropriate company that fulfills the country's desires. Our probabilistic model is more different and flexible, where it helps the countries to choose the suitable company. It allows us to select a suitable queue from a set of K-independent M/M/1/N queues (waiting lists). The quality control is applied for each queue by using an inspector. This inspector returns the country (client) back to the service when something is wrong. In addition, this model takes into account the probability of the reneged country and the balking country to enter the waiting list without obtaining the service. This model provided us the ability to find a suitable company with the maximum probability and minimum cost. We consider this cost as a random variable effort that has a known probability distribution. This model contributes to establish a strong link between search theory, queuing theory, operations research, management, and decision-making science. Also, it presents the solution of one of the most difficult stochastic and discrete optimization problems to make an appropriate decision and then increasing the control over the management of the system. This paper is organized as follows: Section 2 provides the assumptions and the notations which are used in our model. Section 3 describes our model and presents a system of basic probabilistic differential-difference equations corresponding to each company. More than solving this system of the equation under a steady-state situation, we provide the detection probability function which depends on this solution. In Section 4, we present a difficult discrete stochastic optimization problem and its solution gives the minimum searching effort and the maximum service rate for the suitable company under the quality control process. Besides that, an application provides the effectiveness of our model in Section 5, the optimization of the economic model which gives the total expected profit of each company. Finally, the results, discussion, and the future work are presented. This system is formulated from a set of K-independent queues , = 1, 2, … , , where the quality is monitored on their performance. We consider the following notations: The queue number , = 1, 2, … , The queue number service rate The queue number arrival rate Steady-state probability when the queue number has units The following assumptions are used to describe the behavior of each , = 1, 2, … , as follows: 1. Countries (units or clients) enter , = 1, 2, … , one by one according to a Poisson process with the rate (> 0) such that 1∕ is the mean interarrival time. 2. The service times of the countries are exponentially distributed random variables (independent and identical distributed random variables) with the rate (> 0) and mean service time 1∕ , where 0 < < . 3. Obtaining the vaccine from each company (queue) is subject to the principle of first come first served. 4. In , = 1, 2, … , , after the country gets the vaccine, it either joins at the end of the original , = 1, 2, … , as a feedback country with probability (1 − ) (the probability that a defect done) or departure the queue with probability , 0 ≤ < 1. 5. The inspection event (responsible for controlling the quality of the service) in the queue number when there are jobs: . This event is defined from a given policy where the inspector examines all serviced units to ensure the quality of service. 1. During this race to obtain the vaccine, the country waits a period until the service will start with a probability 1 − . In this situation, if the country waits for a long period of time and does not start this service, then the country will loss and leave , = 1, 2, … , without obtaining the service with a probability ( − 1) , 2 ≤ ≤ where the time is a very important factor. 2. The country may join , = 1, 2, … , with a probability or balk with a probability 1 − when countries already exist in , 0 ≤ < 1 for 1 ≤ ≤ − 1 and = 1 otherwise. Naturally, there is an increase in the demand for a safe and inexpensive vaccine against COVID-19. Thus, its price will be varied from one company to another. Each country determines the specifications and the vaccine type based on its capabilities. Due to some problems such as delays in performing the service or the high costs, some countries refrain from purchasing the vaccine that they contracted for after a waiting period. The country is striving to choose an appropriate vaccine in terms of the costs and the effectiveness against COVID-19. All companies and institutions are considered independent from each other in this regard. Therefore, the country chooses an appropriate company that guarantees the quality of the vaccine. Each company (i.e., queue or the waiting list) contains a server and an inspector to ensure the quality of the provided services as shown in Figure 1 . Among these independent companies, an appropriate company is selected. According to the previous assumptions, we have the basic probabilistic differential-difference equations corresponding to each queue , = 1, 2, … , (see Kotb and El-Ashkar 35 ): Under the steady-state situation, the above system will take the following form: F I G U R E 1 Choosing an appropriate ∕ ∕1∕ , = 1, 2, … , queue from K-independent queues under the quality control process Kotb and El-Ashkar 35 solved this system by using the iterative method and obtained the probability of units or countries , = 1, 2, … , as follows: where Since and give the average rates of entering and serving the country in , = 1, 2, … , , respectively, then one can consider a measure of traffic congestion (traffic intensity) for -server as where 0 < < 1. Consequently, (9) and (10) will become functions of and as follows: It is clear that the events of choosing one from , = 1, 2, … , are mutually exclusive events. Thus, when there are no countries in the system (i.e., there is no country in 1 , 2 , … , ). If = 1, 2, … , is full then it is normal for all courtiers to leave it upon their arrival. This is one of the important reasons that the country aims to search for a suitable company to get the vaccine. Therefore, we assume that the country moves according to the Markov process under favorable conditions to choose the suitable company. If the service is done, then the conditional probability of detecting the suitable company , = 1, 2, … , will be 1 − ( , , ( ) ) where ( ) shows the amount of the searching effort (see El-Hadidy 30 ). Thus, one can consider the function ( ) = ∑ = 1 ∑ =1 ( ) , as a cost function of detection where 0 ≤ ( ) ≤ ( ) ≤ such that ∑ =1 ( ) ≤ and be a random variable having a known probability distribution. As in El-Hadidy, 30 we can get the probability of detection for the suitable company under the quality control process by From the independence principle, each country aims to minimize the cost of detection for a suitable company. This will participate in maximizing the quality control of the system. As in Hong et al., 36, 37 we consider the exponential detection function, ( , , ( ) ) = − ( ) . Thus, (14) will become From (11) and (12) in (15), we can obtain the probability of no detection for the suitable company under the quality control process bỹ The decision-makers need a set of factors that the country tries to strike a balance between them in order to be able to find a suitable company. One of the important factors is the detection cost of the suitable company under the quality control process. This factor participates in maximizing the service quality. Therefore, we need to maximize ( ) which is equivalent to the minimization of̃( ). This will minimize the values of the searching effort ( ) , = 1, 2, … , . On the other hand, the minimization of the traffic intensity is equivalent to the maximization of the service rate , where = . This is another objective that we strive to find the appropriate (not crowded) queue. This will provide the service as quickly as possible. In our problem, we have the condition ( ) ≤ , then one can let ( ( ) ≤ ) ≤ , ∈ [0, 1]. From the central limit theorem, we can obtain ( we present a difficult convex multiobjective nonlinear discrete stochastic optimization problem (wherẽ( ) is an exponential function): To get the optimal (minimum) values of ( ) , = 1, 2, … , , we apply the necessary Kuhn-Tucker conditions to solve the above problem. Therefore, we have where , = 1, 2, … , is the Lagrange multiplies. Since ( ) > 0 for any , = 1, 2, … , , then from (21) one can get 3 = 0. It is known that < (i.e., < 1) for each , = 1, 2, … , at any time > 0, then from (22), one can have 4 = 0. For optimality, the probability 1 − ( ) should be greater than zero at any time > 0 then from (20), one can (19), one can get 1 ≠ 0. Thus, the case 1 > 0, 2 = 3 = 4 = 0 is the only case that considers to get the minimum cost * ( ) > 0 and the minimum traffic intensity * . Consequently, (17) and (18) and respectively. From (24), we have From (26) in (19), we get Since In the beginning, some companies will announce that they have developed a vaccine, and most countries will turn to them. But during testing the vaccine, a defect was found and it needs some more careful studies. At this time, countries leave these companies and will head towards other companies that have passed all the tests on the vaccine. Therefore, we can assume here that some lists are completely empty (i.e., = 0). Then, from the formula: there exist at least one of , = 1, 2, … , , equals to zero. Hence, (25) and (27) will become and Also, since − ( ) ≠ 0 and = then in (29) and (30), we have and where = ( ) + Ω √ ( ). Consequently, from (31) and (32) we get If the values of and , , , , , , and ( ) , = 1, 2, … , are generated randomly, then one can use (33) to get the minimum values of * and * numerically. The minimum value of the traffic intensity * means the value of the service rate attains its maximum * which eases the crowding in the suitable company under the quality control process. Since COVID-19 swept all the world, competition has raged between international research and pharmaceutical companies to produce the vaccine earlier. These giant companies have intensified their efforts to finish the clinical trials of the potential vaccines. Therefore, we assume that six international companies with different absorptive capacities , = 1, 2, … , 6 reach the production stage of the vaccine. Naturally, the countries choose vaccines that are safer and less expensive. Some vaccines are safer, but they are more expensive because there are some difficulties such as the way of storage (e.g., Pfizer-Pew&Tech COVID-19 vaccine). Each country chooses the appropriate vaccine, according to its financial capabilities. During the clinical trial periods of the vaccine, some countries retreated their requests from some companies, which had already booked and moved to reserve the same quantity from another company. Our model contributes the country to choose the appropriate company to provide it with the appropriate vaccine at minimum time and cost. For all these hypotheses, let the values of , , , , , , and ( ) , = 1, 2, … , be generated randomly by using Maple 13 (see Table 1 ). We use (33) to obtain the maximum values of * and the minimum values of * numerically as given in Table 1 . Figure 2 presents the plotting of three dimensions of (33) . It presents the relationship between and for , = 1, 2, … , 6 as two intersected planes, where all points in the interested curve are satisfied (33) . Thus, we have an infinite number of solutions that give the optimal values * and * . By using Maple 13, we deduce an equation of and (equation of the intersected curve) which presents the intersection between the two planes. Also, the relationship between and appears in Figure 3 . Since 0 < < then one can determine the minimum value * which corresponds to the maximum value * for all , = 1, 2, … , 6. These optimal values depend on the maximum probabilities * ( ) and * 0 for each queue , = 1, 2, … , 6. We use (9) and (10) to get these probabilities (see Table 1 ). In addition, to obtain the maximum probability of detection * ( ), we substitute in (16) to get̃( ) = 0.00330037 and then obtain the maximum value of * ( ) = 1 −̃( ) =0.9966996308. From Table 1 , some countries consider 6 as a suitable company where its cost is minimum ( * 6 ≅ 4.967348633) and it has a service rate * 3 ≅ 0.825 (see Figure 3F ). There are many factors that made this company the best despite its large The intersection between plans which presents the solution of (33) capacity 6 = 40. Although 6 = 27 (it is a large number of countries corresponds to 6 capacity) the probabilities * 06 ≅ 1 and * ( 6 )6 ≅ 8 × 10 −15 (i.e., there is no unit in 6 or all units serve until its arrival). Also, this queue has a high rate of departure where 3 = 0.4822218558. On the other hand, the inspection process was carried out many times with this service rate. This means that this company accomplishes the required quantities with minimum time, cost, and with the highest possible accuracy. But its balking rate 6 = 0.46455488may be high compared with the other companies such as 2 , 3 , and 4 . This may be due to the spread of bad reports about vaccine safety after some clinical trials. Anyway, choosing the suitable company is the decision for the country. The optimal values of * and * for different values of TA B L E 2 The optimal values of the performance measures for each queue , = 1, 2, … , 6 * ≅ * ≅ * ≅ * * * Of course, one can use the probabilities * 0 and * ( ) for all = 1, 2, … , 6 to obtain the optimal values of the performance measures for each queue , = 1, 2, … , 6 (see Table 2 ). If the number of countries in the steady state in is ′ , then we can get the optimal expected number of units in the system number = 1, 2, … , 6, where * = ( ′ ) = ∑ ∞ = 1 * ( ) . Also, the optimal expected number of countries in is * = * − (1 − * 0 ). For , the optimal average number of the occupied service is * = * − * = 1 − * 0 . By using the Littel's formulae, 38 we get the optimal expected waiting time in the system and the queue number = 1, 2, … , 6 by * = * ∕ and * = * ∕ , respectively. Thus, the optimal expected service time is * = * − * = ( * − * )∕ = (1 − * 0 )∕ . Also, the average reneging and retention rates of each queue are given by = ∑ = 1 ( − 1) * ( ) and = ∑ = 1 ( − 1) * ( ) ; respectively. All these values are presented in Table 2 . Also, for the system, one can use (13) to obtain 0 = ∏ 6 =0 0 = 7.906698800 × 10 −60 . Figure 4 presents a nonlinear relationship between * and , = 1, 2, … , 6. This shows the effect on the optimal expected number of units in the system number = 1, 2, … , 6 with the random values . This relationship is insufficient to judge on the suitable queue. In the above calculations, we take into consideration the profit of the clients only, and we neglect it for the company. It is assumed that there will be a profit of the company producing the vaccine in order to continue providing its services. Thus, we need to strike a balance between all these hypotheses. To do this balance, we deal with an optimization model which was used by Kotb and El-Ashkar 35 to select a suitable company. Adequate decision making is a measurement method that presents the governments' ability to accomplish difficult tasks. The previous part concerns only with finding the suitable company with the minimum cost * and the maximum service rate * and do not take into its consideration the profit which returns to the company. This profit is an important factor that contributes to the best decision. With this profit, the company will have the ability to continue with the vaccine production. Thus, we take this condition in our consideration when we choose the suitable company. To identify the suitable company, we calculate (i) the total expected cost per unit time for each queue , = 1, 2, … , 6 from the following: Table 1 and the values , from Table 2 , we consider = 2, 1 = 12, ℎ = 3, = 1, = 2, and = 3. Thus, the values of , , and are given in Table 3 . As in Figure 5 , the companies with the waiting list 4 and 5 have the maximum total expected profit per unit time. With these companies, the country ensures that the production continues with high safety ratios of the produced vaccine. Now, the decision is to choose the best among them. From Table 1 , we found the company with the waiting list 4 is more suitable, where * 4 is smaller than * 5 . The departure rates 4 and 5 are some convergent, respectively. The values of 4 , 4 are smaller than 5 , 5 . This model presents more quality control in finding the suitable company with minimum effort. At the same time, it is keeping the company from losing. A model that presented here increases the speed of finding the suitable company with minimum effort to get a vaccine of COVID-19. It also provides more quality control on the production process of this vaccine. The country's financial and economic conditions are taken into consideration. The waiting list of each company contains a set of countries. Each country has an exponential service time and a Poisson arrival process. Under a steady-state situation, we solved an interesting and difficult discrete stochastic optimization problem to get the minimum detection effort. This effort has a known bounded probability distribution. An application has been presented to show how this model presented more quality control in finding the suitable company with minimum effort. At the same time, it is keeping this company from losing. Also, the optimal performance measures, the highest service rate, and the total expected profit of each queuing system are obtained. In the future work, we can consider the solution of this model under transient behavior to support the decision making about choosing the suitable company with maximum quality control. If we consider the stochastic volatility parameter in this model then the companies will be avoiding the sudden loss. Moreover, this model can be extendable to cover a multiserver setting in the case of multiple production lines for each company. The data used in the application example are available from the MAPLE 13 programs used here can be requested from the authors. Mohamed Abd Allah El-Hadidy https://orcid.org/0000-0002-9407-9586 On the linear search problem More on the linear search problem Son of the linear search problem Coordinated search for a random walk target motion Existence of cooperative search technique to find a Brownian target Cooperative search model for finding a Brownian target on the real line Detection of the diffusion nanoparticle in the turbulent flows using the random walk model Searching for a d-dimensional Brownian target with multiple sensors Study on the three players linear rendezvous search problem Existence of a finite multiplicative search plan with random distances and velocities to find a d-dimensional Brownian target Generalised linear search plan for a D-dimensional random walk target On optimal coordinated search technique to find a randomly located target The generalized linear search problem, existence of optimal search search paths Coordinated search for an object on the line Coordinated two searchers for an object hidden on an interval Coordinated search for a randomly located target on the plane On the coordinated search problem on the plane 3-Dimensional coordinated search technique for a randomly located target On probabilistic modeling and feasibility of collision between a randomly moving meteor and satellite Existence of finite parabolic spiral search plan for a Brownian target Minimizing the expected search time of finding the hidden object by maximizing the discount effort reward search Existence of optimal N-slinky-turn-spiral search Paths for finding N-dimensional Brownian target Fuzzy optimal search plan for N-dimensional randomly moving target Optimal multiplicative Bayesian search for a lost target On maximum discounted effort reward search problem The searching algorithm for detecting a Markovian target based on maximizing the discounted effort reward search Study on the finiteness of the first meeting time between N-dimensional Gaussian jump and Brownian diffusion particles in the fluid Studying the finiteness of the first meeting time between Lévy flight jump and Brownian particles in the fluid reactive anomalous transport A mathematical model for preventing HIV virus from proliferating inside CD4 T brownian cell using Gaussian jump nanorobot Cooperation in service systems How to share the cost of cooperating queues in a tandem network Cost sharing of cooperating queues in a Jackson network Line balancing in parallel m/m/1 lines and loss systems as cooperative games /N queue with balking and retention of reneged customers Optimal search-relocation trade-off in Markovian-target searching A pseudo-polynomial heuristic for path-constrained discrete-time Markovian-target search A proof for the queueing formulae Fundamentals of Queueing Theory Optimal 3-dimensional search model to find the underwater randomly hidden target Spiral with line segment directory for a helix search path to find a randomly located target in the space Mohamed Abd Allah El-Hadidy is an associate professor of mathematical statistics in the Department of Mathematics, Faculty of Science, Tanta University, Tanta, Egypt and an associate professor of mathematical statistics in the Department of Mathematics and Statistics, Faculty of Science, Taibah University, Yanbu, Saudi Arabia. His research interests include search theory, probability theory, applications of stochastic processes, and optimization. He has published many articles related to these areas in international journals. Ajab A. Alfreedi is an assistant professor of statistics in the Department of Mathematics and Statistics, Faculty of Science, Taibah University, Yanbu, Saudi Arabia. His research interests include forecasting and probability theory. He has published many articles related to these areas in international journals.How to cite this article: El-Hadidy MAA, AA Alfreedi Detection of an appropriate pharmaceutical company to get a suitable vaccine against COVID-19 with minimum cost under the quality control process. Qual Reliab Eng Int. 2021;1-19.