One Note for Fractionation and Increase for Mixed-Level Designs When the Levels Are Not Multiple

: Mixed-level designs have a wide application in the ﬁelds of medicine, science, and agricul-ture, being very useful for experiments where there are both, quantitative, and qualitative factors. Traditional construction methods often make use of complex programing specialized software and powerful computer equipment. This article is focused on a subgroup of these designs in which none of the factor levels are multiples of each other, which we have called pure asymmetrical arrays. For this subgroup we present two algorithms of zero computational cost: the ﬁrst with capacity to build fractions of a desired size; and the second, a strategy to increase these fractions with M additional new runs determined by the experimenter; this is an advantage over the folding methods presented in the literature in which at least half of the initial runs are required. In both algorithms, the constructed fractions are comparable to those showed in the literature as the best in terms of balance and orthogonality.


Introduction
One of the goals of experimentation is to establish the form of the relationships that allow accurate data to be obtained for design purposes. Based on this objective, Wilkie in 1962 addressed the need to use mixed-level designs with 6 or 8 levels for one or more factors, and presented a case study, as well as a statistical analysis [1]. Mixed-level designs are commonly used in different applications, especially when factors are qualitative. Mixedlevel designs are defined as those in which the factors have different numbers of levels [2,3]. Included in this definition there are two cases regarding the factors' levels: when the levels are equal for all factors, for example D(4 5 ), these designs are called pure or symmetrical; and when the levels are different for some other factor, for example, D(2 4 3 1 4 1 ), these are called mixed or asymmetrical [4]. Within the group of asymmetrical designs, there are two subgroups with different characteristics: the first of these is a design in which some of its levels are multiples of each other, for example, D(3 1 6 1 7 1 ). The second is a design in which none if its levels are multiples of each other, for example, D(3 1 5 1 7 1 ). For this research, we have focused on the second subgroup, which we have called pure asymmetrical arrays.
Practical success when using mixed-level designs is due to efficient use of experimental runs to study many factors simultaneously [5]. Fractional factorial designs are the most popular designs in experimental investigation [6]. Traditional construction methods for mixed-level fractional factorial designs often make use of complex programming, specialized software, and powerful computer equipment; see [2,3,[7][8][9][10][11][12][13][14]. For this, an important number of criteria has been developed to measure the balance and orthogonality properties as quality attributes [3,7,8,[15][16][17]. Even giving way to some comparisons between them [5,18] with different applications such as those described in [1,10,[19][20][21] as well as techniques to perform augmentations, in which the minimal requirement is to add a number of runs equal to 50% of the initial design size see [22,23]. According to the literature, there is an area of opportunity that must be attended in favor of the development of an algorithm with zero computational cost that allows the construction of fractions with the best levels of balance and orthogonality. Two situations are of particular interest: (1) when these designs form a design themselves; and (2) when these designs are joined to other designs to form a new design for example: an orthogonal fraction of the (n, 2 4 ) and a semi-orthogonal fraction (n, 3 1 5 1 7 1 ) can form a semi-orthogonal fraction for the design (n, 2 4 3 1 5 1 7 1 ); this is a common practice to form a mixed-level fraction [9].
At present, the use mixed-level fractional factorial designs in early stages of experimentation opens an important possibility within the oriented use of resources (i.e., human resources, raw materials, machinery, among others). Allowing the experimenter to scrutinize the influential effects in an economic scenario, with advantages such as allocating resources, obtaining results in a shorter time, reducing the impact of machine deterioration and equipment, among many others. Although the advantages in the use of mixed-level fractional factorial designs are widely known, the use of these designs has been limited because the exiting techniques for generating these fractions require the use of tools that require extensive domain and investment (i.e., complex methods, specialized computer equipment, specialist labor, specialized software, among others.).
There is an interest in the development of an instrument that breaks with the need for these additional resources. This research offers a zero computational cost tool that expands the tools currently offered by the state of the art, providing the experimenter with an easy to understand and apply method that does not require complex programming and can be used by anyone with basic knowledge of statistics, and therefore facilitating the implementation of mixed-level fractional factorial design in different fields of study. Pantoja et al. (2019) developed the NOBA (near-orthogonal balanced array) method to generate mixed-level fractional factorial designs balanced-orthogonal and semi-orthogonal, the study showed that a percentage of the designs analyzed proved to be "infractionable" due to nature of its factors [24]. Several examples of these designs, including 2 to 6 factors, are shown in the Tables 1 and 2. In these designs, several of the levels are not multiples of each other. Therefore, the least common multiple of the levels is equal to the number of runs of the design matrix. When choosing a design of pure asymmetrical arrays to be fractioned, the size of this array stops being a multiple of at least one of the factors levels. Thus, this method is only able to generate near-orthogonal, near-balanced arrays. For this reason, the fractions generated are called near-orthogonal, near-balanced pure asymmetrical arrays (NONBPAs). This group is clearly the least studied since fractions belonging to this group have been only published in [3]. In this work it is possible to see the concept of efficient array (EA), the design with the best possible balance and orthogonality properties. EAs have been obtained from the application of genetic algorithms and the optimization of an objective function resulting from the sum of the standardized J 2 -optimality and the standardized balanced coefficient (Form II). It is in this context, and when considering the possibility that a NONBPA could be required in any field of application just as much as any other design, that the importance of studying NONBPAs became evident.
Consider a shoe manufacturing company in which the implementation of a NONBPA is required. The objective is to evaluate different materials for a new shoe concept, focused on users with foot pathologies. Two response variables are of interest: pathological benefits and production costs. The required design is (2 1 3 1 5 1 7 1 ) and the factors to consider are: buttress material, lining, type sole, and slipper material. Table 3 shows the design levels, in this case, the alternative of running a full factorial (210 runs) was ruled out due to projected costs and required times. The decision was to run a NONBPA consisting in only 20 runs (9.5% of the full factorial).     A notable contribution from this research is the development of two algorithms of zero computational cost. The first algorithm allows the construction of a NONBPA fraction and the second algorithm provides a strategy to increase these fractions with M additional runs. Both designs, the original NONBPA and its augmented version were compared to the EAs presented in [3]. The results showed that the NONBPAs are just as good as the EAs in terms of GBM (general balanced metric), J 2 (orthogonal parameter), and V IFs (Average variance inflation factors).
The paper has been organized as follows: Section 1 presents the introduction and motivation. Section 2 presents several new concepts and two algorithms (NONBPA structure, method to build a NONBPA and an example, as well as a strategy to increase NONBPAs with M additional runs). In Section 3, a comparison of NONBPAs vs. EAs is provided. Section 4 presents a practical application, and finally, the conclusions are presented in Section 5.

Mixed-Level Fractional Factorial Designs
The study of orthogonal arrays has been the focus of many investigations; two desirable properties for these arrays are balance and orthogonality. Orthogonal arrays contain pairs of linearly independent columns and are useful to evaluate the importance of several factors. Orthogonality ensures that the effects can be estimate independently [7]. For a matrix to be balanced, in each column, each possible factor level must appear the same number of times. Columns whose levels do not appear with the same frequency are called unbalanced. The concept of near-balanced denotes that, although not all levels appear equally due to design size limitations, all levels appear with the most similar frequency. The importance of preserving the balance lies in the fact that executing the same number of times each level of a factor in an experiment, results in a uniform distribution of information for each level. Thus, there is consistency in the variances of the difference of observations in pairs of treatment combinations [3].
Mixed-level fractional factorial designs have led to the continued generation of parameters to measure the quality of these arrays. Xu and Wu (2001) developed the generalized minimum aberration (GMA) for comparing asymmetrical fractional factorial designs. This criterion is independent of the choice of treatment contrasts, and thus model free and it is applicable to symmetrical and asymmetrical designs [15]. Xu (2003) proposed the minimum moment aberration (MMA) to assess the goodness of nonregular designs and supersaturated designs [16]. Xu and Deng (2005) proposed the moment aberration projection (MAP) to rank and classify nonregular designs, it measures the goodness of a design through moments of the number of coincidences between the rows of its projection designs [17]. Xu (2002) presented the J 2 parameter (see Section 1.3) [7]. Dean and Lewis (2006) offered an important revision of this criteria from the minimum aberration criterial approach [21]. Liu et al. (2006) generalized χ2 (D) criterial and investigated connections between GMA, MMA, and MAP criteria [5]. Guo et al. (2007) defined the balanced coefficient criterion for main effects and used it as an objective function to measure the degree of balance and orthogonality of a near orthogonal array generated by using genetic algorithms; in this research he presents a catalog of 20 arrays also called EAs; one characteristic of these designs is they require a reduced number of runs while preserving high levels for balance and orthogonality [3].  extended the balance coefficient beyond main effects giving rise to the GBM a minimum aberration criterion that can be used to evaluate and compare mixed-level fractional factorial designs [8] (see, Section 2.3).
Methods for construction of mixed-level fractional factorial designs include Wang and Wu (1992), they proposed an approach for construction of orthogonal designs based upon difference matrices [10]. Wang (1996) presented a method for construction of orthogonal asymmetrical arrays through the generalized Kronecker sum mixed-level matrix and mixed difference matrices [11]. Nguyen (1996) presented a method to augment orthogonal arrays with additional columns in such a way that the resulting design possesses good level for E and other criteria [19]. DeCock and Stufken (2000) designed and algorithm for construction of orthogonal mixed-level design through searching some existing two-level orthogonal designs [25]. Xu (2002) developed an algorithm to add columns sequentially to a design by using the generalized minimum aberration and minimum moment aberration criteria [7]. Salawu (2012) used the balanced coefficient and J 2 optimality criteria to compare the two forms of balanced coefficient methods using the generalized minimum aberration and minimum moment aberration criteria [26]. Fontana (2017) presented a methodology based on the joint use of polynomial counting function, complex counting of levels and algorithms for quadratic optimization [13]. Grömping and Fontana (2018) proposed an algorithm for generation of mixed-level arrays with generalized minimum aberration using mixed integer optimization with conic quadratic constraints [14]. Pantoja et al. (2019) developed the NOBA method, an algorithm based on divisor factors and permuted vectors that can generate mixed-level fractional factorial designs [24].
One consequence of using a fractional factorial design is the aliasing of factorial effects. A standard follows up strategy involves adding a second fraction called foldover. A foldover can be constructed for various reasons. If the analysis of the initial design reveals that a particular set of main effects and interaction are significant, the foldover design can be chosen to resolve confounding problems; if one factor is very important, it should not be confused with other factors. On the other hand, if the goal is to dealias all, or as many as possible main effects from 2FIs, or 2FIs from each other [27,28]. A full foldover consists of adding a second fraction of the same size as the initial fraction, obtained by inverting the signs of one or more columns two-level designs or by rotating one or more columns (for three-level and mixed level designs) [29].
The foldover is only one of several augmentation techniques developed for two-level designs, other techniques include semifold, D-optimal semifold, quarterfold, and R3 algorithm. Sequential experimentation techniques for mixed-level designs include foldover [22] and semifold [23]. The foldover is constructed by rotating columns and the semifold by performing exhaustive research. The foldover technique is computationally more efficient when compared to searching for additional runs in the full factorial, which could not be practical. The main disadvantage of this method is that it requires the same number of runs as the initial array and the size of the augmented design may be large in some situations. In order to reduce the number of runs required by a foldover, the concept of semifold was introduced making it possible to reduce the foldover plans to half the number of runs. [23].

General Balanced Metric and Balanced Columns
Balanced columns contain all levels equally often. Therefore, a balanced matrix for main effects has a value of GBM = 0 (Equation (5)). Columns whose levels do not appear equally often are called unbalanced. The concept of near-balanced denotes that while not all levels appear equally often, due to the size limitations, all levels appear as equally often as possible. Therefore, both balanced and near-balanced designs are considered to have optimal balanced status given the constraint on the number of runs. An unbalanced column is considered not near-balanced when it is neither balanced nor near-balanced [8]. Ghosh and Chowdhury mentioned the importance of balance for achieving some or all treatment contrasts estimated with the same variance, they also mentioned the importance of common variance (CV) designs when the objective is to discriminate between two models having common as well as uncommon parameter. This paper emphasizes the major role played by the uncommon parameters and generalizes the concept of CV designs when there are at most k (≥1) uncommon parameters. They also introduce a new concept of "Robust CV designs for replications" having the possibility of replicated observations and demonstrate the robustness for equally replicated observations. In addition, two general designs for three level symmetric factorial experiments are presented [30].  defines the GBM as a measure of the degree of balance for both, main effects and interactions in a mixed-level design [8]. It is defined as an n × k design matrix d, n is the number of rows and k is the number of factors. Let d t (t = 1, . . . , k) denote matrices including all t-factor interaction columns, and d 1 is the one-factor-interaction matrix for the main effects. Note that d 1 is equivalent to d. Therefore, the whole interaction matrix involves all t-factor interaction matrices d t . That is (see Equation (1)), Let l t j be the number of levels of the j jh column in d t (1 ≤ t ≤ k). Let c t rj be the number of times the rth levels appears in the jth column of d t .
T be the counts for each level for the jth column of d t . The notation H t is used for the balance coefficient of d t . We can employ a distance function to reflect the degree of balance and define the jth columns balance coefficient as shown in Equation (2), for the k-factor interaction matrix, where T t j = n l t j is fixed. Substituting T t j = n l t j , then H t j becomes in the Equation (3), The balance coefficients H t for d t just sum the H t j and are defined as shown in Equation (4), Then, the GBM can be defined as in Equation (5), For two designs d 1 and d 2 , suppose r is the smallest value such that H r (d 1 ) = H r (d 2 ). Say that d 1 is more general balanced than d 2 if H r (d 1 ) < H r (d 2 ). If no design is more general balanced than d 1 , then d 1 is said to be the most general balanced design. To calculate the value of the GBM parameter, consider that H j t (Equation (2)) represents the error between the frequencies with which each level appears with respect to the frequency with which it should appear. Therefore, it is notable that for a semi-balanced column H j t > 0 and said value will tend to increase when the frequency of one or more the levels in that column moves away from the mean, which in this context corresponds to the frequency with each level should appear.

J 2 and VIFs for Orthogonal Arrays
The J 2 optimality parameter was proposed by Xu [7]. For an N x n matrix d = [x ik ], weight w k > 0 is assigned for column k, which has s k levels. For 1 ≤ i, j ≤ N, let (see Equation (6), where δ(x, y) = 1 if x = y and 0 otherwise. The δ i,j (d) value measures the similarity between the ith and jth rows of d. In particular, if w k = 1 is chosen for all k, then δ i,j (d) is the number of coincidences between the ith and jth rows. Defined in the Equation (7), A design is J 2 -optimal if it minimizes J 2 . Obviously, by minimizing J 2 (d), it is desired that the rows of d be as dissimilar as possible.
For an N × n matrix d whose k th columns has s k levels and weight w k , and the equality holds if and only of d is OA (see Equation (8)).
L(n); is the minimum value that is reached by J 2 when a matrix is orthogonal. Therefore, since the NONBPAs are semi-orthogonal arrays, the value of L(n) cannot be considered as a reference point to minimize J 2 . A more direct comparison is achieved by calculating the V IFs.
VIF (variance inflation factor), of the predictor x j is calculated based on the linear relationship between the predictor x j and the other independent variables [x 1 , x 2 , . . . , x j-1 , x j+1 , . . . , x m ]. As shown in Equation (9).
where, R j 2 is the coefficient of determination of the regression of x j on all other independent variables in the data set [x 1 , x 2 , . . . , x j-1 , x j+1 , . . . , x m ] (see Equation (10)).
As it is known if the value of VIF = 1; then el coefficient of determination R j 2 = 0 and the predictors are not correlated, if 1 ≤ VIF ≤ 5; the predictors are moderately correlated and if VIF > 10 indicates that the correlation between predictors is excessively influencing the regression results. VIFs are easy to interpret since the higher the VIFs value, the greater the correlation between the predictors [31,32].

NONBPA Structure
A NONBPA is a fraction of the model matrix formed by k columns and n rows, in which column A has symbols of (1,2, . . . ,la) column B has symbols of (1,2, . . . ,lb) and so on. A balanced N-column is a column that contains NT N number of times each element of the vector VL N , formed by the levels present in the N-column, see Equation (11). In addition, a near-balanced N-column is a column formed by two segments: the first segment (balanced segment) is formed by NT N ; number of times the VL N (vector of levels for the N-factor) and the second segment (non-balanced segment) by the vector of complementary levels for the N-factor (VCL N ); which is formed by elements from 1 to SVCL N (where SVCL N is the size of the vector of complementary levels for the N-factor), see Equation (12). Figure 1 shows the structure for a balanced and a near-balanced N-column. For a given N-column, where NT N value is the ratio of n with respect to the number of levels in the column (ln) (see Equation (11)). If NT N is integer the vector complementary levels does not exist and the column is balanced, in the other case; the column is nearbalanced and the vector of complementary levels exists.
The size of the vector of complementary levels for a column near-balanced is SVCL N , defined in the Equation (12).

Method to Build a NONBPA
This section shows the NONBPA method see Figure 2. The method consists in 3 steps described below. Step 1. Select a pure asymmetrical array. Select a mixed-level design in which factor levels are not multiples of each other.
Step 2. Define the size of the array, n. Determine n ensuring that this is equal or greater to the necessary degrees of freedom needed. For example, for a design with 4 qualitative factors with 5, 6, 7, and 9 levels, the minimum degrees of freedom required to estimate all effects are: 4 + 5 + 6 + 8 (main effects) + 1 (intercept) + 1 (error) = 25, the smallest fraction that can be constructed is size 25.
Step 3. Construction of NONBPA. For the ith factor, replicate the vector 1, . . . , li (with li, the number of levels) until n runs have been assembled. The last vector can be completed (balanced column) or cut before being completed (near-balanced column).
Example, Construction of NONBPA (24,5 1 6 1 7 1 ) Consider a situation that involves three qualitative factors: A, B, and C with 5, 6, and 7 levels, respectively. The experimenter is interested in running a fraction, given that factor levels are not multiple of each other, he decides to use the NONBPA method. Figure 3 shows step by step the construction of a NONBPA. In step 1, the design selected is the D(5 1 6 1 7 1 ). In Step 2, the size of the array is determined by the experimenter, in this case, a size of 24 was chosen, n = 24.
Step 3 consists in the construction of array. First, we will mention column B, note that last vector can be completed, therefore the column is balanced. For columns A and C, the last vector is cut before being completed, for this reason columns A and C are near-balanced.
If the experimenter is interested in keeping a specific factor balanced, this can be done by changing the size of the array. Table 4 shows several possibilities of n for this design. Note that to preserve the balance property, it is advisable to select an array size that is a multiple of la × lb . . . . × ln, where la, lb . . . . ln are the levels of the factors we want to be balanced in the fraction. This table can also be useful to select n, for example, n = 30, keeps factors A and B balanced while maintaining a reasonable fraction size.  Figure 4 shows the GBM calculation for the design (3 1 5 1 7 1 ) note that the column of factor B is a balanced column. Therefore, H j t = 0. That is, all the levels for this column appear with the same frequency C t j = [4,4,4,4,4,4] T . On the other hand, since columns A and C are semi-balanced, they have values of H j t > 0. That is, they have a contribution of 0.8 and 1.71 respectively and GBM = 0.80 + 0 + 1.71 = 2.51. Regarding orthogonality, J 2 =112 while V IFs = 1.01. Therefore, NONBPA (5 1 6 1 7 1 ) is a semi-orthogonal array in which the predictors are minimally correlated.

Augmentation Strategies vs. NONBPAs
Ghosh and Rao (1996) presented a comprehensive study on sequential assembly of fractions (see Figure 5) [33]. They used the design presented by Box [34] (p. 394), this design is presented in T 1 ; it contains 7 factors and 8 runs. T 2 is obtained from T 1 by switching the signs of column 4. T 3 is obtained from T 1 by switching the signs of all columns and T 4 is presented as T 2 with the columns (1,2,3) as (4, 5, 6), (4, 5, 6) as (1,2,3) and the runs (rows) are also in different order. T 4 is known as a Search Design [35]. Then a series of augmentation options are presented and evaluated in terms of balance and orthogonality.
The coincidence between these matrices and the NONBPAs is notable. Because in the NONBPAs in a similar way; when renaming the levels of one or more factors it is possible to construct a new fraction. Original NONBPAs and new NONBPAs whit renamed levels have similar properties for balance and orthogonality. Figure 6 shows the level rotation of factor B of NONBPA (2 1 3 1 5 1 ), note that two additional designs have been generated.
NONBPAs, like many other designs can be augmented. In this section we show a simple strategy to augment the NONBPAs with M additional new runs determined by the experimenter. To increase a NONBPA, it is enough to decide M and then to add the M additional rows by rotating factor levels.   Example, Augmenting a NONBPA (15,2 1 3 1 5 1 7 1 11 1 ) with 7 Additional Runs Figure 7 shows a NONBPA (15,2 1 3 1 5 1 7 1 11 1 ) augmented with M = 7 additional runs. Note that the first 15 runs belong to the original NONBPA and the new M = 7 runs are added by rotating factor levels. Note that creating a NONBPA of size n + M from scratch would produce the same fraction. Table 5 shows a comparison between the NONBPA (15,5 1 6 1 7 1 ) and its augmented versions. In both designs the number of balanced columns is 2. Increasing the number of runs directly benefits columns A and E, minimally effects B and C and has no effect on D. This provoked a significant reduction in the GBM value from 3.90 to 2.72 improving the balance. Regarding the orthogonality property, the values of J 2 and V IFs were calculated; as it is already known, J 2 increases as the number of runs increases. Therefore, a direct comparison is not possible. A direct comparison is achieved with the V IFs; from 1.03 to 1.01, note that the design becomes more orthogonal as the number of runs increases.    Table 5 shows a comparison between the NONBPA (15,5 1 6 1 7 1 ) and its augmented versions. In both designs the number of balanced columns is 2. Increasing the number of runs directly benefits columns A and E, minimally effects B and C and has no effect on D. This provoked a significant reduction in the GBM value from 3.90 to 2.72 improving the balance. Regarding the orthogonality property, the values of J2 and were calculated; as it is already known, J2 increases as the number of runs increases. Therefore, a direct comparison is not possible. A direct comparison is achieved with the ; from 1.03 to 1.01, note that the design becomes more orthogonal as the number of runs increases.

Results Comparison of NONBPAs vs. EAs
NONBPAs were compared with its competitors, the EAs, presented in [3]. To perform the comparison, four EAs were selected, including EA(21,3 1 4 1 7 1 ), EA(20,3 1 4 1 5 1 ), EA(24,5 1 6 1 7 1 ), and EA(15,2 1 3 1 5 1 7 1 ). The equivalent NONBPAs were constructed and augmented so that the number of runs of the NONBPAs were equal to the number of runs of the EAs. In this way, a more direct comparison was possible. To compare the designs, the

Results Comparison of NONBPAs vs. EAs
NONBPAs were compared with its competitors, the EAs, presented in [3]. To perform the comparison, four EAs were selected, including EA(21,3 1 4 1 7 1 ), EA(20,3 1 4 1 5 1 ), EA(24,5 1 6 1 7 1 ), and EA(15,2 1 3 1 5 1 7 1 ). The equivalent NONBPAs were constructed and augmented so that the number of runs of the NONBPAs were equal to the number of runs of the EAs. In this way, a more direct comparison was possible. To compare the designs, the balance (GBM) and orthogonality (J 2 and V IFs) are measured. Figures 8 and 9 show that the EAs and the NONBPAs have similar levels for balance and orthogonality, the difference is only minimal only with respect to V IFs values. It was also observed that GBM and J 2 remain equal for EAs and NONBPAs even if the size of the array changes. Figure 10 shows an efficient array for a D(3 1 5 1 7 1 ) compared to its corresponding NONBPA. The comparison was made by using 15, 21, and 30 runs. Note that in all cases, the GBM and the J 2 were identical for EAs and NONBPAs, and V IFs values are very similar.

Practical Application
To demonstrate de capacity of the NONBPAs to estimate factorial effects, the method was compared to full factorial and EA using simulated data ( Figure 11). The design selected to perform the comparisons was the D(210,2 1 3 1 5 1 7 1 ). To generate the simulated data, a simple model with the form Y = 6 [A] + 10 [C] + ε(0, σ 2 ), was used and an experimental error was introduced, which is a random variable with zero mean (u = 0) and variance (σ 2 ) [27]. Based on research performed by Ríos et al. (2011), the size of the variance must be one third of the regression coefficient for the regressor to be reported as significant [23].  Table 6 shows the ANOVA tables and the optimization for the three designs. Results are very consistent for the three methods; ANOVA tables look similar, and all the designs were able to detect A and C as significant effects. Regarding the optimization, a desirability function for maximization was used and the three designs produced the same recommended levels for factors A and C, which are two and five, respectively. In all cases, the optimal value for the response is very similar.

Conclusions
Industrial experiments often involve situations in which categorical and numerical factors with different numbers of levels are present, these experiments are commonly known as mixed-level designs. Mixed-level designs require a high number of runs and are difficult to carry out because of the cost and time required. One alternative to avoid running a full factorial is to run a mixed-level fractional factorial design. Unfortunately, these fractions are not easy to construct because they often require complex programming techniques, specialized software, and expensive computer equipment.
The new method presented here, called NONBPA, is an algorithm capable of generating mixed-level fractional factorial designs when the factor levels are not multiple of each other. The near-orthogonal near-balanced pure symmetrical arrays generated are extremely flexible in run size and possess high levels of balance and orthogonality. The arrays generated with this method were compared to the EAs presented in [3] and the results showed that the balance and orthogonality property were identical for both methods. In addition to the construction method, a method to perform augmentations was also provided, this method allows augmenting any NONBPA with M additional runs while preserving the balance and orthogonality properties.
The main advantages of the NONBPA method are that it is easy to understand and apply, it does not require complex programming, the computational cost is low and it can be used by any person with basic knowledge in statistics.
GBM and V IFs are parameters that allow to compare respectively, balance and orthogonality between arrays with the same or different number of runs. On the other hand, J 2 only allows the comparison of the level of orthogonality between arrays that have same number of runs. A disadvantage in the use of J 2 for NONBPAs is that it is not possible to know the minimum value of L(n) for semi-orthogonal arrays. Therefore, for the NONBPAs the use of V IFs is recommended.
Future research for the NONBPA will focus on evaluating balance and orthogonality beyond main effects, opening a greater number of possibilities for experimenters in the various fields of application.  Acknowledgments: We thank the editor and reviewers for their helpful comments and suggestion that greatly improved the content and quality of paper.

Conflicts of Interest:
The authors declare no conflict of interest.

Balance coefficient L(n)
The minimum value that is reached by J 2 when a matrix is orthogonal CV Common variance R j 2 Coefficient of determination of the regression NT N Number of times the VL N VL N Vector of levels for the N-factor VCL N The vector of complementary levels for the N-factor SVCL N The size of the vector of complementary levels for a column near-balanced