Application of Central Composite Design to Optimize Spawns Propagation

Despite the increased recognition of the nutritional value of the Oyster mushroom, coupled with its ability to tolerate a wide range of climatic conditions, its production is still at infancy stage with low adoption rate in Kenya. The low uptake could be attributed to the cost of spawns or lack of skills for spawns preparations coupled with poor knowledge on oyster mushroom consumption benefits. The objective of this study was to optimize Pleurotus ostreatus (Oyster mushroom) spawns production. To achieve the objective, the spawns propagation was optimized by varying the temperature level, sterilization time and culture media concentration in order to establish the feasible levels which minimized the days of mycelium full development using central composite designs. Based on the study findings, 26.29 ̊C, 17.36 minutes and 60.95 g/L of temperature level, sterilization time and culture media concentration levels respectively minimized the days to full coverage of mycelium in a petri dish. Central composite designs for controlling temperature, sterilization time and culture media concentration were recommended for spawns maximum production. A further research on multiple response optimizations aimed at achieving resistance to bacterial diseases and yield by varying the strain in the culture were recommended.

k is the number of independent variables in which the design can be replicated f n times. The star points consist of 2k units on the axis of each factor at a distance α, from the centre of the design [(±α, 0, …, 0), (0, ±α, 0, …, 0)], whose selection is based on the orthogonality and rotatability criterion, and takes one observation at each of the vector can be replicated s n times, where i e is the i-th Euclidean unit vector and 0 α > . The centre points designated as "0" points and can be replicated c n times. Thus by letting n denote the total number of experimental runs in the CCD, based on k design factors, the runs size can be determined by using Equation (1) where: n is the desired sample size. c n is the number of times centre points that were replicated. There are various types of central composite designs depending on where the star points are placed, namely: circumscribed, inscribed and face centered central composite designs. Circumscribed Composite Designs (CCC) has a circular, spherical or hyperspherical symmetry. In the circumstances where the limits specified for settings are truly the limits, then it is referred to as the Inscribed Composite Design (CCI) which uses the factor settings as the star points and creates a factorial or fractional factorial design within those limits. Then lastly Face Centered Composite Design (CCF) has an alpha of 1, such that the axial points are at the center of each face of the factorial space, hence, 1 α ± . Figure 1 shows the graphical representation of the three variants of central composite designs;

Spawns and Oyster Mushroom Production
Mushrooms are macrofungi with a distinctive fruiting body which can be either epigeous or hypogeous. The macrofungus has fruiting body large enough to be seen with the naked eye and can be picked up by hand [3]. Spawn, which is the reproducing part of the mushrooms, should be multiplied in a pure culture media in a sterilized environment. The pure culture can either be raised by tissue culture or spore culture. In tissue culture a well grown mushroom with membrane covering the gills is selected, from which a small bit of mushroom from gill portion is taken using forceps and inoculated on Potato Dextrose Agar (PDA) or Malt Extract Agar (MEA). Under spore culture method, the spores are collected from well-developed fruiting body by "spore mapping technique" and then the spores are inoculated to the PDA or MEA slants as in tissue culture under aseptic condition. The mycelia (spawns) are supposed to cover the entire surface after which they are placed in a sterilized jars or bottles with a boiled grain seeds like wheat or millet for further multiplication. It is at this stage when the seeds can be used as mushroom spawn and can produce mushrooms if mixed with a well compost substrate.

Problem Statement
Despite the increased recognition of the nutritional value of the Oyster mushroom, coupled with its ability to tolerate a wide range of climatic conditions, its production in Kenya is still at infancy stage with low adoption rate [4]. The low uptake could be attributed to the cost of buying spawns which could be substituted by the acquisition of the skills to prepare them, coupled with poor knowledge on its production and consumption benefits.
Past studies indicates that, beside the correctly compost substrate, the oyster mushroom does well under certain levels of temperature, humidity, light and ventilation, while the quality of spawn is affected by the mother culture [5], media culture concentrations and spawn substrate preparations [6]. Other several studies indicate that oyster mushrooms yield depends on the type of substrate used and the quality of the spawns [7] [8].
Hence this study focused on determining the optimal level of temperature, sterilization time for media culture and the media culture concentration settings, that minimized the time (days) to form mycelium culture establishment during the oyster mushroom spawn preparation, by using central composite design.
This would reduce the process, production costs and narrow the gap of opportunistic competitors' invasion, hence improve on the oyster mushroom's spawn quality, production and consumption levels.

Literature Review
Central Composite Design has widely been used in various fields of research; it Open Journal of Optimization was used in optimizing machine parameters in drilling hybrid metal matrix composites, and it was found that the predicted values and measured values were fairly close, indicating that the developed models could be effectively be used to predict the responses in the drilling of hybrid metal matrix composites, [9].
CCD's are very popular 2 nd -order designs because; they are extremely simple to use, and allows estimation of all the parameters in a full second-order model coupled with the fact that among the exact (integer) designs, they often have high efficiencies under the commonly used A-, D-and E-optimality criteria [10].
Composite Design and Response Surface Methodology were applied to optimize the conditions of submerged cultivation of the Fomesfomentarius mushroom. Concentration of dextrose, yeast extract and time of cultivation in days were selected and their correlative effect on mushroom multiplication established. The maximum yield of dry weight biomass was 23.74 g/L with 0.8 -7.5 g/L concentration after 11 days of submerged cultivation [11].
CCD was used to show that, the evaporative cooling of the air surrounding the pileus created convective airflows capable of carrying spores at speeds of centimetres per second. That work revealed how mushrooms tolerate and even benefit from crowding and explained their high water needs. It was evident that spores continuously flow out from thin gaps, even in the absence of external winds [12]. CCD was used to optimize the Degradation of Tylosin from Aqueous solution by Photo-Fenton Reaction. The interaction effects and optimal parameters were obtained and the significance of the independent variables and their interactions was tested by means of analysis of variance (ANOVA) at 5% significance level. The results showed that the concentration of the ferrous ion and pH were the main parameters affecting Total Organic Carbon (TOC) removal, while peroxide concentration had a slight effect on the reaction. The optimum operating conditions to achieve maximum TOC removal were determined. There was a good agreement between the model prediction and experimental results confirming the soundness of the developed model [13].
Response surface methodology was applied to optimize the main factors which significantly affected exopolysaccharide (EPS) production. The concentrations of glucose and peptone were found to be the main effective factors for EPS production by the fractional factorial design (FFD) and central composite design experimental analysis. Verification experiment confirmed the validity with the actual EPS yield as 1.97 g/L, which was 6.29-fold in comparison with the (2.22 g/L) in the original basal medium [14].
Central composite design and response surface methodology were used in the advanced treatment of olive oil processing wastewater using Fenton's peroxida- Under this study, the central composite design was used for theoretical examination of optimal temperature, sterilization time and the media culture concentration experimental space that minimized the time (days) to full spawn growth.

Methodology
The growth process of the spawns was a function of temperature level, sterilization time and the culture media concentration level as shown in Equation (2) ( ) , , The objective function was to determine the level of the regressors' experimental setting that minimized the mycelia full coverage in the petri dish area through the response surface methodologies.
A second order response surface design should be rotatable, to make it possible to extract the most information regarding the dependent variable, and leave the least amount of uncertainty for the prediction of future values [16].
Rotatable design provides the preferred property of constant prediction variance at all points that are equidistant from the design center, hence improving the quality of the prediction. It is a desirable property, especially when there is a need to optimize ( )y x over the region of interest. A design is said to be orthogonal if it can provide independent information about the effects of the various terms in the model.
In order to ensure orthogonality and rotatability of the design, circumscribed composite design was preferred. This ensured that any non-allowable operating conditions at two or more of the extremes of the design region were encompassed.

Rotatability Conditions
The design matrix for a CCD experiment involving k independent variables of L levels, is a matrix derived from the values corresponding to the three types of experimental runs stated in Equation (2) to form a design with N of N L treatment combinations matrix, such that usually referred to as the design matrix. Where, the treatment combinations are the points of the design, such that i x denotes the i th factor at j th treatment in a given level. Such a design will be a rotatable design of order d if a response polynomial surface is of Equation (3) form [16].
which can be so fitted such that the variance, implying the variance of ˆu y is a function of only the distance of x from the origin but not of the direction and therefore constant at all points that are equidistant from the design center.
When the model is of the second degree, that is 2 d = , then such constancy of variance is attained if the design points are chosen to satisfy the following conditions: where 3 λ and 4 λ are constants of any design. Hence, where the summation of the relations in Equation (4)  Therefore by using Equation (4) and the rotatability conditions, Equation (5) was obtained [

The First-Order Model
The first order model aids in getting the direction of the relationship between the response and independent variables. It leads the experimenter sequentially and efficiently along the path of steepest ascent/descent. The first order model is  (6).
where i x and k β 's are the design variables and regression coefficients respectively.
When the curvature appears, the second-order model is used to approximate the response values precisely.

The Second-Order Model
The central composite design was used to explore the region for fitting in the first and second order models in this study. The prevailing region of exploration was 20˚ -30˚ degree Celsius for temperature, 10 -20 minutes of sterilization time and 35 -65 grams of agar per 1000 cubic centimeter of PDA. To simplify the calculations the independent variables were coded (−1, 1) interval such that given 1 ξ , 2 ξ and 3 ξ as the natural variables for temperature, time and culture media concentration levels respectively, then the coded variables should be derived at by using the Equation (7), relations [16] [20].
To enable in the identification of the best CCD by creating a grid of all combinations of the design choices, the information function for the second-order (quadratic) model had to be rotatable, the α values necessary for orthogonality and rotatability were computed. Ideally, in order to make an unbiased estimate of pure error, the CCD should comprise of three to five centre points [21]. The central composite design is commonly used to fit the second order response surface model of the form expressed in Equation (8) where; u y is the response obtained from the u th combination of factors ( x denotes the level of the i th factor ( 1, 2, ,

The Significant Regression Coefficients
To test the significance of the individual regression coefficient j β , the P-Value and the t test were used. The t test statistics is based on the t distribution and has the form of Equation (9) Open Journal of Optimization where ˆj β is the least square estimator of the parameter j β ( 1, 2, , j q =  ) and ( )j s β is the estimated standard error of ˆj β . The standard error of each parameter ˆj β is given by the square root of the diagonal elements of the matrix, and X is the matrix of values of explanatory variables referred to as the design matrix, such that.

Location and Characterization of the Stationary Points
The temperature level ( 1 x ), Sterilization time ( 2 x ) and Culture media concentration level ( 3 x ), that minimized the days of Pleurotus ostreatus spawns development were established by carrying out the partial derivatives from the function ( ) and by letting 0 , , , whereby the second order model in matrix notation is as summarized in Equation (12) such that the fitted values in matrix form are; whereby b is a (k × 1) vector of first order regression coefficients and B is a (k × k) symmetric matrix whose main diagonal elements are pure quadratic coefficients (ˆi i β ) and whose off-diagonal elements are one-half the mixed quadratic Open Journal of Optimization Therefore at the stationary points ˆ2 0 y b Bx x ∂ = + = ∂ and the stationary points would be and the predicted response at the stationary points would be determined by Equation (13) 1 0 The second order model normally takes the maximum, minimum, saddle point or stationary ridge. The coordinates and the nature of each stationary point were determined by carrying out the partial derivatives and represented by using the contour maps.

Canonical Analysis
The where { }

Optimization of Spawns Propagation
The operating optimal levels of temperature, sterilization time and culture media concentration that minimized time in days of the mycelia full coverage in a petri dish area were determined through central composite designs. The colonised media in a petri dish is displayed in Figure 2  The Petri dishes were regularly inspected and the contaminated ones were immediately discarded. The inset (b) of Figure 2 shows a contaminated Petri dish while the inset (d) is a fully colonized Petri dish. The mycelia growth covered the entire substrate within the 10 -20 days of inoculation, by which the spawns were ready for use.

Analysis of a Second-Order Model
A model that incorporated the curvature was appropriate to approximate the response, given that the region of optimization was identified, the second order model was fitted to determine the optimum set of operating levels for the temperature, time of sterilization and the culture media concentration that minimized the days to full mycelium colonization.

Rotatable Second-Order Designs
To make a design both (approximately) orthogonal and rotatable, the axial distance for rotatability should be chosen, and then the center points are added [18]. Rotatability is a desirable and necessary property in response surface methodology for second order models exploration [16].
Under this design criterion, equidistant from the design center ( 1, 2, , i k =  ). Therefore by applying Equation (1) the experimental runs were 19 and it was possible to have a rotatable and orthogonal design by adding the star-points and centre points replicated five times to the simple square factorial design points; Such that given k = 3, p = 0, respectively for temperature, as per Equation (7).
The actual data value for each input variable was sequentially altered and noted, as the process and the response surface was observed. The length of the colony diameter of mycelia was noted daily and the time to full colonization was observed as summarized in Table 1. The rows in Table 1 represented the runs while the columns represented the levels of the factors. The first row was run at the "low" level of all of the factors, the second row was run at the "high" level of temperature as a factor but the "low" levels of sterilization time and agar media concentration.
Although the variable names remained 1 x , 2 x and 3 x and design order 1 -19, the actual experiment implementation was randomized. All the necessary properties for CCD second-order (quadratic) model such as rotatability and orthogonality were observed. To make an unbiased estimate of pure error, the CCD comprised of five centre points.
Alteration of the temperature level, sterilization time and the culture media concentration that shortened the mycelia colonization time (days), were as summarized in Table 2.
The shortest time observed was six (6) days with nineteen (19) days being the longest time. By using the coded second model data, it was possible to get the design matrix (X), the model equation and carry out the computation of the criteria and their efficiencies as discussed in the following section.

Optimality Criteria
The general three factors rotatable CCD was modified by varying the number of Table 2. Coded variable design data.
By using A-criterion, whereby a design is said to be A-optimal if it minimizes the sum of the main diagonal element (trace) of ( ) The sum of the diagonal was 0.7445928 in all the designs. The value changed only due to the number of experimental runs.
With regard to the D-optimal, the value that maximized X X ′ was 7.48 and it was constant among different number of design centre points. Finally, regarding the E-optimal the value that minimized the maximum eigenvalue of ( )  Table 3.
Therefore the addition of center points on CCD second-order designs had no significant effects on the optimality. However, the five center points were chosen to improve in making an unbiased estimate of pure error, since for an unbiased estimate the CCD should comprise of three to five centre points [22].

The Second Order Model Regression Equation
The second degree order model was computed using the R commands and the outputs were as given in Tables 4-7.
Therefore based on the second order model regression output, the model Eq-

Significant Predictors
Based on the t and the p-values it was observed that the entire predictor variables were statistically significant for influencing the time to spawn full development but at a varied levels of significance. Temperature was only significant at 5% but sterilization time and media culture concentration at 1% significance level. The temperature quadratic term was statistically significant at 1% significance     Table 6. Stationary point of response surface. level. The coefficient of determination R 2 indicated that the model could approximate the data on the design points at 80.75% (after adjusting 61.5%). Table 5 outputs display the first-order models (FO), which specifies the first-order response surface (i.e., a linear function), the canonical analysis of the Open Journal of Optimization Since the lack of fit was statistically insignificant ( Pr 0.078 0.05 = > ) then the model could predict the response variable appropriately.
The response surface stationary points output was in coded form, which could be transformed into their natural value through the equations; The information in Tables 4-7 was computer generated. However, the same results could be arrived at through manual calculations as discussed in the following section.

Canonical Analysis
The fitted quadratic model provided a noticeable response variation with clarity and the canonical analysis gave the coordinates of the estimated stationary points and the canonical directions from the points.
The fitted canonical form model is characterized by ( ) Now suppose we let   x w x Pw

Stationary Points
To determine the levels of temperature ( 1 x ), sterilization time ( 2 x ) and culture media concentration level ( 3 x ) that minimized time to full colonization of the agar medium, the partial derivatives defined in Equation (10) were conducted such that;    Table 8.
Therefore the response surface could be predicted by using the following equ-

Diagnostics and Plots of the Estimated Response Surfaces
The basic assumptions on errors which include, independence, normality and the constancy of variance for errors in order to apply the second order statistical model and ANOVA test statistics were tested. The findings are summarized in Figure 3. The fact that there were equally spread residuals around the horizontal line without any distinct patterns in the residual versus fitted values plot, then that was an indication of residuals linear or non-linear relationships nonexistence, thereby satisfying the assumption of residuals independence. The normal quantile-quantile plot shows the residuals were normally distributed given that the observed values were along the reference line hence ascertained normality. The scale-Location plot showed that the residuals were equally spread along the ranges of predictors and hence the assumption of equal variance (homoscedasticity) was not violated. Based on the residuals versus leverage plot, most of the predictor variables were well inside of the Cook's distance lines hence few identified outliers. The following response surface graphs display the predictor variables  on the x-and y-axes and then a continuous surface that represents the fitted response values on the z-axis. For each of the surface map, two variables were displayed at a time holding the third variable constant. Figure 4 shows the response surface map and a contour plot for temperature and time of substrate sterilization for spawn's multiplication.
They represent a minimum response surface thus indicating a point of optimum operating conditions reaching a minimum. By inspection the plots are non-linear, implying a strong temperature and sterilization interaction effect on the time to spawns maturity and the response surface is curved due to the fact that the model contains quadratic terms that were statistically significant.
The contour lines appeared to follow the direction of movement along the path of minimum response from a reference point. The contour lines were within the response values in the data which ranged between 6 to19 nevertheless it is worth noting that the accuracy of the response surface plot depends on how well the model represents the true relationships among the variables.
The lowest values of days to full colonization of the substrate in the petri dish were in the upper side of the contour plot as it appears in Figure 5, which corresponds with high values of both agar media concentration and temperature, time of sterilization was held constant for this plot.
As with the contour plot the lowest values of the response surface were along the side opposite the temperature axis, which corresponded with high values agar media concentration and temperature  The following and the last is the contour plot and then the response surface map for agar media concentration and time of sterilization.
As in the contour plot in Figure 6, the lowest values of days to full substrate colonization are in the lower right corner of the response surface, which corresponded with high values of agar media concentration.
The lowest values of days to full substrate colonization are in the upper right corner of the plot, which corresponds with high values of agar media concentration but average time of sterilization, temperature was held constant in this particular analysis and representation. Open Journal of Optimization

Conclusion and Recommendations
The CCD provided the statistical elements necessary for the evaluation of Temperature, agar concentration and sterilization time during spawn propagation.
The contour plots and the response surface maps generated characterized surface and helped a great deal in locating the optimum points precisely at a glance.
Routinely training on how one would process the spawns using the locally available agricultural by-products or forest wastes, using the scientifically proven formulations was timely and would not only provide a protein-rich food but reduces the environmental pollution. It is a transformative link for inedible wastes into edible biomass of high monetary value.
Multiple response optimizations aim to achieve maximal nutritional value and yield against minimal duration of spawns propagation and the spawns quality is suggested and recommended.