The Construction of a Model-Robust IV-Optimal Mixture Designs Using a Genetic Algorithm

: Among the numerous alphabetical optimality criteria is the IV-criterion that is focused on prediction variance. We propose a new criterion, called the weighted IV-optimality. It is similar to IV-optimality, because the researcher must ﬁrst specify a model. However, unlike IV-optimality, a suite of “reduced” models is also proposed if the original model is misspeciﬁed via over-parameterization. In this research, weighted IV-optimality is applied to mixture experiments with a set of prior weights assigned to the potential mixture models of interest. To address the issue of implementation, a genetic algorithm was developed to generate weighted IV-optimal mixture designs that are robust across multiple models. In our examples, we assign models with p parameters to have equal weights, but weights will vary based on varying p . Fraction-of-design-space (FDS) plots are used to compare the performance of an experimental design in terms of the prediction variance properties. An illustrating example is presented. The result shows that the GA-generated designs studied are robust across a set of potential mixture models.


Introduction
Industrial product formulations (e.g., food processing, chemical formulations, textile fibers, and pharmaceutical drugs) frequently involve the blending of multiple mixture components.Mixture experiments form a special class of response surface experiments in which the product under investigation is comprised of several components.In this research, we assume that the response of interest is a function only of the proportions of the components that are present in the mixture.The levels of the q experimental factors (x i ; i = 1, 2, . . ., q) in a mixture experiment are component proportions.Thus, each x i is between zero and one, and the sum of the q component proportions is one.Under these conditions, the experimental region involving q proportions is a regular q − 1-dimensional simplex.Typically, there are also single-component constraints (SCCs) defined by lower L i and upper U i bounds on the component proportions: If SCCs are imposed on the component proportions, the experimental region will now be an irregularly shaped polyhedron within the simplex.For a general review of the design and analysis of mixture experiments, see [1,2].When the mixture region is a simplex, the standard mixture designs such as simplex lattice and simplex centroid designs are often suitable for implementation.However, they are not applicable in situations involving SCCs.When SCCs exist, an extreme vertices design (McLean and Anderson [3]) is one possibility.However, these designs may be very inefficient with respect to prediction variance properties.To improve the prediction variance properties, an algorithmically generated design focused on optimizing an optimality design criterion is commonly used.Optimal designs, however, are, in general, optimal only for a specified model.Numerous approaches have been developed for constructing optimal designs, such as ACED [4], XVERT [5], CONSIM [6], Fedorov [7], DETMAX [8], exchange algorithms [9,10], and genetic algorithms (GAs) [11][12][13][14][15]. GAs have the advantage of being adaptive search algorithms [14].Recent applications of GAs provide alternative approaches to classic exchange-point algorithms to generate designs.Examples of using GAs to generate designs can be found in Borkowski [11], Heredia-Langner et al. [12,13], Park et al. [15], and Limmun et al. [16].
In this research, we develop and employ a GA that extends the work in Limmun et al. [16] that focused on a single model when constructing an IV-optimal design where the experimental region was an irregularly shaped polyhedral region that is a subspace of a simplex.Their procedure offered a degree of flexibility in its way of constructing designs that allowed it to overcome restrictions that may limit the applicability of other algorithms.The first aim of our research is to consider model misspecification and introduce a weighted IV-optimality criterion determined over a set of potential mixture models, and it will serve to address our second aim: develop a design-generating GA using weighted IV-optimality as its objective function.
The rest of this paper is organized as follows.In Section 2, the Scheffé mixture models and relevant theory relating to mixture experiments of interest are reviewed, and weighted IV-efficiency is defined.Section 3 includes a brief introduction to GAs for constructing designs and an illustration of the steps used in our GA.In Section 4, the proposed scheme is demonstrated with examples, and with concluding statements in Section 5.

Notation and Models
The most common forms of mixture models are the Scheffé (canonical) polynomials.For example, the Scheffé linear model is given by and the Scheffé quadratic model is given by In these models, y is the response variable; each β i coefficient represents the expected response when the proportion of the i-th component equals unity; each β ij coefficient represents the nonlinear blending properties of the i-th and the j-th component proportions, and ε is a random error.In this paper, we focus on the Scheffé quadratic model mixture model and possible model misspecification.
The Scheffé mixture models are linear models and can be expressed in the matrix form as where y = (y 1 , y 2 , . . . ,y n ) is the response vector, X is the N × p model matrix with columns associated with the p model terms (such as linear and cross-products terms), β is the p × 1 vector of model parameters, and ε = (ε 1 , ε 2 , . . . ,ε n ) is the vector of random errors associated with natural variation of y around the underlying surface assumed to be independent and identically normally distributed with zero mean and variance σ 2 I n .The prediction properties of a design, specifically the scaled prediction variance, is dependent on the chosen model through model matrix X.The scaled prediction variance v(x 0 ) is defined as where x 0 is the expansion of a mixture x = x 1 , x 2 , . . ., x q to vector form corresponding to the p columns of X; N is the design size, which penalizes prediction variances for larger designs, and Var( ŷ(x)) is the variance of the estimated response at x.For example, if there are three components and the model is the Scheffé quadratic model, then x = (x 1 , x 2 , x 3 ), and X is an N × 6 model matrix.

Optimality Criteria
Design optimality criteria are single valued criteria that represent different variance or parameter estimation properties of a design.Several criteria have been advanced with the purpose of constructing and comparing and designs.Four commonly used optimality criteria are the D, A, G, and IV criteria.These four criteria are functions of the information matrix (X X) for a given design.The D and A criteria are focused on parameter estimation, while G and IV criteria are focused on the prediction variance.When a design is being considered for implementation, several of its properties can be evaluated by computing measures of design efficiency.Common D, A, G, and IV design optimality measures are where χ is the design space and V is the volume of χ.The D and A criteria are the simplest to handle computationally, and they were the criteria considered in the earliest design generating algorithms (e.g., exchange algorithms).Because designs with high IV or G efficiencies also tend to have good D-and A-efficiencies, this paper will focus on designs that minimize the average prediction variance (i.e., minimize the denominator in the IV-efficiency over the entire experimental region).Designs that minimize the IV-optimality criterion include the IV-optimal design, the Q-optimal design [17], the V-optimal design [18], and the I-optimal design [19,20].In a review of software approaches to evaluating average prediction variances (APVs), Borkowski [21] showed that it is common for software packages to estimate the APV = χ v(x 0 )dx 1 dx 2 . . .dx q by taking the sample mean of v(x 0 ) over a fixed set of points (e.g., the candidate set for an exchange algorithm).He demonstrated that estimation of the APV based instead on a random set of evaluation points is unbiased and superior to estimation based on a fixed set of evaluation points.This is one flaw in using software using exchange algorithms to generate IV-optimal designs: the estimated APV is an overestimation of the integral.In our proposed GA, the APV measure is calculated by averaging v(x 0 ) over a random set of points and will provide an unbiased estimate of the IV-optimality criterion.Additionally, the variability of the estimator will decrease as the size of the random evaluation set increases.In this paper, we use 5000 random points in the evaluation set.Several authors have provided results for the IV-optimal designs, including Borkowski [11], Syafitri et al. [20], and Coetzer and Haines [22].Further details on the motivation and uses of the optimality criteria can be found in Box and Draper [23], Atkinson et al. [18], and Fedorov [7].

Weighted IV-Optimality
In this paper, we develop and propose the weighted IV-efficiency, which is defined as where ..dx q is the IV-efficiency for model i, χ is the design space, V is the volume of χ, N is the number of design points, r is the number of reduced models for a given full model, and w i is the weight for model i.
In terms of design generation, the goal is to use weighted IV-optimality to find a set of points that will minimize the weighted average of the average of the scaled prediction variance over the design region across a set of reduced models.Practically speaking, the goal is to generate a design that protects against having a final model with poor prediction variance properties.Similar to finding an IV-optimal design, the researcher must specify a model.This serves as the "full" model when we consider the weighted IV-optimality.However, unlike IV-optimality, a suite of "reduced" models are also proposed if the full model is misspecified via over-parameterization.Examining Equation (6), note that weighted IV-optimality is not restricted to a particular full model nor is it restricted to mixture models.It easily generalizes to other response surface designs and the associated models that can be fitted by those designs.
The experimenter has the freedom to choose the weighting factors for the full model, the most parsimonious or smallest model, as well as all other intermediate models.To exemplify the use of weighted IV-optimality, a "full" quadratic mixture model and a "smallest" linear mixture model will be considered.
Although there are numerous ways to assign weights in the weighted optimality criterion, in this paper, we assume (i) that not all models should be weighted equally and (ii) that only models having an equal number of model parameters receive equal weight.Here are the reasons for assuming (i) and (ii).Before running the experiment, the experimenter believes that the full model may be the most appropriate, so he/she chooses the maximum weighting factor for the full model.In the analysis phase, however, the full model might be inappropriate because of misspecification due to overparameterization. Thus, it seems reasonable that the weights be reduced as we move further away from the full (or the experimenter's most likely a priori) model via model reduction, and stop with the smallest weight assigned to the model with the fewest number of parameters.We also treat each parameter as equally important under the assumption that the researcher has no prior knowledge or make an educated guess regarding which terms would most likely be removed if a model reduction occurs.Therefore, we assign uniform weights among all models resulting after removal of one term, and then assign a decreasing but uniform weight to all models resulting when two terms are moved, and so on.Therefore, models with more parameters have weights that reflect their greater importance relative to models with fewer parameters.For the proposed method of calculating a weighted IV-efficiency, it is necessary to have a complete enumeration of the set of subset models (reduced models) of interest.If an experimenter has justification for another weighting scheme, then it should certainly be implemented.
One specific weighting scheme consistent with (i) and (ii) above is to use weights for each reduced model based on the numbers of model parameters.Suppose model i has j parameters.The weight we assigned to model i is then defined as where ψ j is the jth weighting factor, and j = q, q + 1, . . ., s − 1, s, where q is the number of parameters of the linear mixture model in Equation (2); s is the number of parameters of the quadratic model in Equation (3), and m j is the number of reduced models with j parameters.The model i weight (w i , i = 1, 2, . . ., r) is nonnegative and the weights satisfy ∑ r i=1 w i = 1.The relative weighting factor for computing the weighting factors (supplied by the experimenter) is defined as Math.Comput.Appl.2018, 23, 25

of 13
In this paper, we assume that the experimenter weights the quadratic mixture model 100 times the weight to be assigned to the linear mixture model; therefore, we use R = 100 as the relative weighting factor.The weighting factors are positive and subject to the restriction There are k equispaced levels of the weighting factor, and k = s − q + 1 is the number of levels for a weighting factor.Therefore, the range of the weighting factor is defined as The increment value of the weighting factor can be expressed as The values of the weighting factor can be represented as follows: Again, the weighting factor must sum to one.The weighting factor can be rewritten in the form: Therefore, the minimum weighting factor can be expressed as After simplifying, the increment value of the weighting factor is defined as For example, suppose the full model is the quadratic mixture model with three component proportions.This model can be expressed as There are 8 reduced models (the full model and the seven reduced models as shown in Table 1).A 1 or a 0 in any "Terms in Model" column represents the presence or absence of that term in the model, respectively.The column p is the number of model parameters.Using 100 as the relative weighting factor, the increment value of the weighting factor is 0.1634.The weighting factor value and the weight for each model are shown in the ψ j and w i columns, respectively.If the 2nd, 3rd, and 4th models have five parameters, then m 5 = 3, and the 5-parameter model weighting factor ψ 5 = 0.3318 and the individual model weights w 2 , w 3 , and w 4 = ψ 5 /3 = 0.1106.

Development of the Genetic Algorithm
Heuristic optimization has been used to solve a variety of experimental design problems.One popular optimization algorithm that is based on the general principle of local improvement is the genetic algorithm (GA).The genetic algorithm was developed by John Holland in 1975 and was popularized through the work of Goldberg in 1989 [24].Since then, GAs have been used to solve optimization problems for many applications because it is very efficient over a variety of search spaces.Bäck et al. [25] mention that GAs often yield excellent results when applied to complex optimization problems where other methods are either not applicable or turn out to be unsatisfactory.GAs have been successfully applied to experimental design problems using various optimality criteria.For an introduction to GAs, see Michaelewicz [26] and Haupt and Haupt [27].
A GA is a search and optimization technique developed by mimicking the evolutionary principles and the chromosomal processing in the natural selection.A GA takes an initial population of potential solutions to a problem (parent chromosomes).Through evolutionary reproduction operators, the current parent population then passes some of its properties (genes) to produce offspring chromosomes.A subset of the best parent and offspring is retained for the next generation, and the reproduction process is repeated for many generations until an acceptable chromosome has evolved.A GA uses an objective function as a measure of a chromosome's fitness as a solution to the problem of interest.Although a GA is not guaranteed to find the global optimum in a finite number of generations, it is less likely to get trapped at a local optimum than traditional gradient-based search methods when the objective function is not smooth and generally well behaved.
GAs can be very useful for construction of a design when the optimality criterion (objective function) is difficult to work with and/or where there are constraints on the experimental region such as in mixture experiments.We now describe a GA used to generate optimal designs for a mixture experiment with single component constraints (SCCs) that is based on weighted IV-optimality, which extends and modifies the GA approach of Limmun et al. [16] who generated optimal designs for a specified mixture model.Throughout this research paper, we have encoded GA chromosomes using real-value encoding instead of another encoding (e.g., binary) because real-value encoding is easy to interpret, can be modified for many applications, and is flexible enough to allow for a unique representation for every variable.
A chromosome C will represent a potential design (solution) and is encoded as an N × q matrix, where N is the number of the design points, and q is the number of mixture components.Each row in C is a gene x i = x i1 x i2 . . .x iq and represents one experimental mixture.The objective function F measures a chromosome's fitness as a solution to the function that we wish to optimize.That is, F takes a chromosome as input and outputs an objective function value.Our GA's objective function F is a weighted IV-efficiency.The goal using a weighted IV-efficiency is to find a set of points that will minimize the weighted average of the scaled prediction variance throughout the entire experimental region for a set of reduced models representing possible model misspecification.Specifically, the goal is to find a design that maximizes objective function F where where ..dx q is the IV-efficiency for model i, χ is the design space, V is the volume of χ, r is the number of reduced models for a given full model, and w i is the weight for model i.The GA that we employed includes several genetic operators.Prior to running a GA, the experimenter must specify the r models of interest and the model weights w i .The steps and operators in our GA are as follows: Step 1 Specify the GA parameters: population size (M), number of iterations (I), selection method (elitism with random parent pairing), blending rate (α b ), crossover rates (α bc and α wc ), and mutation rate (α m ).
Step 2 Generate an initial population of M chromosomes (mixture designs).We use the real-value encoding with four decimal places to encode each chromosome.Assume M is odd.To generate the initial population, a random sample is first taken in a hypercube.Then each sampled point in a hypercube is mapped to the constrained mixture space by applying the function used by Borkowski and Piepel [28].Each experimental mixture is recorded to four-decimal place accuracy.
Step 3 Calculate the IV-efficiency objective function F for each chromosome in the initial population.
Step 4 Find the elite chromosome which is the chromosome that has the largest weighted IV-efficiency F. The remaining M − 1 chromosomes are randomly paired for the reproduction process.
Step 5 Produce offspring of the next generation by using genetic operators: blending, between-parent crossover, within-parent crossover, and mutation.Larger values of genetic parameters are used for the early iterations and the smaller values of genetic parameters are used for the later.
Step 6 Calculate objective function F for each parent/offspring pair.The fitter chromosome in the pair is retained for the next generation.
Step 7 Repeat Steps 5 and 6 for I generations.
Step 8 Apply a local grid search to the best design to further improve objective function F yielding the IV-optimal design.A local grid search searches designs in a small neighborhood of the best design.This is accomplished by perturbing the component proportions by small increments to search for further improvements in F. This continues until no further improvement is found.

Illustration: A Genetic Algorithm for Constructing an Optimal Mixture Design
To illustrate our algorithm, consider a three-component mixture (q = 3) with the goal of generating a weighted IV-optimal design having N = 7 design points.We create an initial population of M = 5 chromosomes (C 1 , C 2 , C 3 , C 4 , C 5 ) each having 7 genes where the genes of each chromosome (mixture design) are drawn randomly from the mixture space.One possible realization of the initial population and the associated objective function values are presented in Table 2. C 1 is the elite chromosome because it has the largest F value.Next, the remaining four chromosomes are randomly paired for the reproduction process.Suppose that (C 5 , C 4 ) and (C 2 , C 3 ) are the randomly formed pairs.Applying the reproduction process to C 5 and C 4 yields offspring C * 5 and C * 4 .Table 3 contains a set of random uniform deviates for C 5 and C 4 .Boldfacing indicates a probability test is passed (PTIP) at rates α b , α bc , α wc and α m = 0.02, and that reproduction operation will be performed where α b , α bc , α wc , and α m represent the success probabilities of blending, between-parent crossover, within-parent crossover, and mutation, respectively.The first generation is summarized in Table 4 where boldfacing indicates that a gene was changed.For the blending operator, Row 3 of C 5 is blended with a random row (e.g., Row 4) of C 4 with random blending value β = 0.4438.For the between-parent crossover, Row 6 of C 5 will crossover with a random row (e.g., Row 3) of C 4 .For the within-parent crossover, Row 3 of C 5 and Row 7 of C 4 will have within-row crossovers based on random permutations of (1, 2, 3) (e.g., (3,2,1), and (1, 3, 2), respectively).For mutation, Row 4 of C 5 and Row 1 of C 4 will be mutated.Then, Components 1, 2, and 3 are randomly selected from Row , which now has the largest F value.Hence, chromosome C 1 will be a part of the reproduction process in the second generation.This process will continue for I generations.

A Three-Component Mixture Example
Consider the three-component mixture experiment in Crosier [29].The three SCCs are The boundary is formed by six vertices.The full model under consideration is the quadratic mixture model.Weighted IV optimal designs with N = 7, 10, and 14 points were generated by the GA using a relative weighting factor R = 100.These designs are labeled GA7, GA10, and GA14, respectively.The set of possible reduced models are those in Table 1.We choose M = 25 designs to comprise the GA population of chromosomes.In a preliminary study, we investigated the choice of GA parameter values and the number of generations to observe convergence when running the GA.Based on these results, the number of generations was set to 6000, and the ranges of genetic parameter values were 0.02 ≤ α b ≤ 0.2; 0.03 ≤ α bc , α wc ≤ 0.15; 0.05 ≤ α m ≤ 0.2, where α b , α bc , α wc , and α m represent the success probabilities of blending, between-parent crossover, within-parent crossover, and mutation, respectively.Initially, the parameter values are set to the largest level, and after every 1500 generations, these parameter values are systematically reset to smaller values.The optimal choice operators and the GA parameters across operators is an open research area.
Because no software package can generate weighted IV-optimal designs, we studied how well our weighted IV-optimal designs would fare with respect to D-, A-, G-, and IV-efficiencies.Therefore, we generated IV-optimal designs using the design-generating statistical software package Design-Expert 11 (Stat-Ease (2017)) to generate their versions of IV-optimal designs.These 7, 10, and 14 point designs are labeled DX7, DX10, and DX14, respectively.
The mixtures from the GA and DX designs in the constrained SCC mixture space are shown in Figure 1.The patterns of points are quite different.The majority of points in the GA designs are distributed around the boundary with replications near the overall centroid, while the DX designs also tend to place design points on the boundary but at different locations and the number of replicates differ from the GA designs.Additionally, the DX designs select multiple mixtures as interior points, while the interior points of the GA designs are concentrated with replicates at only one mixture.The fact that the mixtures selected differ between DX and GA designs is based on the fact the DX designs were constructed only to optimize with respect to one model (the quadratic mixture model), while the GA design considered eight potential models when optimizing.Additionally, DX designs were generated using exchange algorithms that restrict mixture selection from only a finite candidate set of mixtures, while GA designs explore the continuum of points in the entire mixture space.The boundary is formed by six vertices.The full model under consideration is the quadratic mixture model.Weighted IV optimal designs with  = 7, 10, and 14 points were generated by the GA using a relative weighting factor  = 100.These designs are labeled GA7, GA10, and GA14, respectively.The set of possible reduced models are those in Table 1.We choose  = 25 designs to comprise the GA population of chromosomes.In a preliminary study, we investigated the choice of GA parameter values and the number of generations to observe convergence when running the GA.Based on these results, the number of generations was set to 6000, and the ranges of genetic parameter values were 0.02 ≤   ≤ 0.2; 0.03 ≤   ,   ≤ 0.15; 0.05 ≤   ≤ 0.2, where   ,   ,   , and   represent the success probabilities of blending, between-parent crossover, within-parent crossover, and mutation, respectively.Initially, the parameter values are set to the largest level, and after every 1500 generations, these parameter values are systematically reset to smaller values.The optimal choice operators and the GA parameters across operators is an open research area.

Initial
Because no software package can generate weighted IV-optimal designs, we studied how well our weighted IV-optimal designs would fare with respect to D-, A-, G-, and IV-efficiencies.Therefore, we generated IV-optimal designs using the design-generating statistical software package Design-Expert 11 (Stat-Ease (2017)) to generate their versions of IV-optimal designs.These 7, 10, and 14 point designs are labeled DX7, DX10, and DX14, respectively.
The mixtures from the GA and DX designs in the constrained SCC mixture space are shown in Figure 1.The patterns of points are quite different.The majority of points in the GA designs are distributed around the boundary with replications near the overall centroid, while the DX designs also tend to place design points on the boundary but at different locations and the number of replicates differ from the GA designs.Additionally, the DX designs select multiple mixtures as interior points, while the interior points of the GA designs are concentrated with replicates at only one mixture.The fact that the mixtures selected differ between DX and GA designs is based on the fact the DX designs were constructed only to optimize with respect to one model (the quadratic mixture model), while the GA design considered eight potential models when optimizing.Additionally, DX designs were generated using exchange algorithms that restrict mixture selection from only a finite candidate set of mixtures, while GA designs explore the continuum of points in the entire mixture space.The next comparison is summarized in Table 6, which contains D-, A-, G-and IV-efficiencies for the six designs.Based on the D-, A, G-, and IV-efficiency, the GA designs are superior to the DX designs.What is interesting about these results is that, even though Design-Expert generated designs with the goal of optimizing the IV-efficiency for the quadratic mixture model, the GA design was more IV-efficient (despite weighting over eight models) for the quadratic model.The final comparison uses fraction-of-design-space (FDS) plots to examine model robustness.A model-robust design should perform consistently well with respect to the scaled predictions in the design region for all possible reduced models.The vertical axis represents possible scaled prediction variance v(x 0 ) values (Equation ( 5)) and the horizontal axis represents the fraction of the mixture design space that has scaled prediction variance values less than or equal to v(x 0 ).Thus, a FDS plot is equivalent to a cumulative distribution function plot but with the axes reversed.For additional information of the FDS plots for examining model robustness, see Ozol-Godfrey et al. [30].
Although an FDS plot can be made for each model, we restrict the study to the quadratic (full) mixture model and the most parsimonious linear mixture model because these curves give a manageable summary of that design's prediction variance performance at the extremes in terms of model parameters.These FDS plots are presented in Figure 2.
The next comparison is summarized in Table 6, which contains D-, A-, G-and IV-efficiencies for the six designs.Based on the D-, A, G-, and IV-efficiency, the GA designs are superior to the DX designs.What is interesting about these results is that, even though Design-Expert generated designs with the goal of optimizing the IV-efficiency for the quadratic mixture model, the GA design was more IV-efficient (despite weighting over eight models) for the quadratic model.The final comparison uses fraction-of-design-space (FDS) plots to examine model robustness.A model-robust design should perform consistently well with respect to the scaled predictions in the design region for all possible reduced models.The vertical axis represents possible scaled prediction variance (  ) values (Equation ( 5)) and the horizontal axis represents the fraction of the mixture design space that has scaled prediction variance values less than or equal to (  ).Thus, a FDS plot is equivalent to a cumulative distribution function plot but with the axes reversed.For additional information of the FDS plots for examining model robustness, see Ozol-Godfrey et al. [30].
Although an FDS plot can be made for each model, we restrict the study to the quadratic (full) mixture model and the most parsimonious linear mixture model because these curves give a manageable summary of that design's prediction variance performance at the extremes in terms of model parameters.These FDS plots are presented in Figure 2.For the linear mixture model, the GA and the DX designs have similar scaled prediction variance distributions for most of the design space.However, for the quadratic mixture model, the DX designs performs poorly for most of the design space compared to the GA designs.This provides evidence for why the IV-efficiencies are smaller for the GA designs compared to the DX designs.For most of the design space, the GA designs appear to have good robustness properties because it has much flatter and lower FDS curves for both the quadratic and linear mixture models, indicating more desirable prediction variance distributions.

Conclusions
In this paper, we introduced the concept of a weighted IV-efficiency, which was then used as the objective function of a GA when generating IV-efficient designs robust to a set of reduced models.The goal of the weighted IV-optimality criterion is to find a set of points that will minimize the weighted average of the average of the scaled prediction variance over the design region across a set of reduced models.The IV-efficiencies are calculated over a set of points randomly selected in the experimental region rather than a fixed set of points to get unbiased estimates of weighted IV-efficiencies.Additional assumptions-(i) that not all models should be weighted equally and (ii) that only models having an equal number of model parameters receive equal weight-were adopted for calculating the weighted IV-efficiency.
We studied the problem of generating weighted IV-optimal designs for mixture experiments where the experimental region in an irregularly shaped polyhedron is formed by single component constraints.The proposed methodology allows movement through a continuous region that includes highly constrained mixture regions and does not require selecting points from a user-defined candidate set of mixtures unlike traditional exchange algorithm approaches.
The GA presented in this research paper is effective for generating model-robust designs.The results from the example show that GA designs performed better than the Design Expert designs when examining the model robustness using FDS plots, as well as the D-, A-, G-and IV-efficiencies.When the experimenter believes that the initial model may not turn out to be the model adopted in the final analysis of the experimental data, GA designs based on weighted IV-efficiency are suggested because they protect against the possibility of having a final model with poor prediction variance properties.
Future research will look for alternatives to using random points to estimate the integrals in the definition of the IV-criterion, additional full and reduced model situations, and alternative schemes for weighting these models when calculating IV-efficiency.

Figure 2 .
Figure 2. The fraction-of-design-space (FDS) plots for all competing designs.Figure 2. The fraction-of-design-space (FDS) plots for all competing designs.

Figure 2 .
Figure 2. The fraction-of-design-space (FDS) plots for all competing designs.Figure 2. The fraction-of-design-space (FDS) plots for all competing designs.

Table 1 .
The set of reduced models the Scheffé quadratic model with three components.

Table 2 .
Initial population of five chromosomes.

Table 3 .
Random deviates for probability tests on C 5 and C 4 .

Table 4 .
First generation reproduction for chromosomes C 5 and C 4 .

Table 5 .
First generation summary.