**Genetic Algorithms and Group Method of Data Handling-Type Neural Networks Applications in Poultry Science**

Majid Mottaghitalb *Dept. of Animal Science, Faculty of Agri. Uni. of Guilan, Rasht, Iran* 

### **1. Introduction**

20 Will-be-set-by-IN-TECH

218 Real-World Applications of Genetic Algorithms

Zitzler, E. & Thiele, L. (1999). Multiobjective evolutionary algorithms: A comparative

*Computation* 3(4): 257–271.

case study and the strength Pareto approach, *IEEE Transactions on Evolutionary*

The necessity of modeling is well established since the structural identification of a process is essential in analysis, control and prediction. Computer modeling is becoming an important tool in different fields in science including Biology. In Artificial Intelligence research, 'intelligence' is increasingly looked upon not as deliberative reasoning processes alone, but as the ability to exhibit adaptive behavior in a complex world. There have been extensive efforts in recent years to deploy population-based stochastic search algorithms such as evolutionary methods to design artificial neural networks since such evolutionary algorithms are

particularly useful for dealing with complex problems having large search spaces with many local optima(Iba, etal,1996). In recent years, the use of artificial neural networks leads to successful application of different type of algorithm in a broad range of areas in engineering, biology, and economics in which GMDH-type is one.

### **2. Genetic algorithms**

Nature employs the best cybernetic systems that can be conceived. In the neurological domain of living beings, the ecological balance involving environmental feedback, or the regulation of the temperature of the human body, are the examples of cybernetic systems of nature that are fascinating in their accuracy and efficiency (Madala and Ivakhnenko, 1994).

In the 1950s and 1960s several computer scientists independently studied evolutionary systems with the idea that evolution could be used as an optimization tool for engineering problems in different systems (as a collection of interacting, diverse elements that function/ communicate within a specified environment to process information to achieve one or more desired objectives) (Mitchell and Forrest 1994).

Evolution can be considered as the first and highest level of adaptation. It involves the adaptation of a species to global ecological and environmental conditions. This adaptation is a relatively slow process that operates over millennia, although the speed of genetic adaptation may differ widely for individual species.

Genetic Algorithms and Group Method

these nodes are also called synapses (See Fig. 1).

Fig. 1. Schematic structure of a Neuron

Fig. 2. A Schematic representation of neural net

neural networks (NN) (Fig. 2).

**3. Neural networks** 

of Data Handling-Type Neural Networks Applications in Poultry Science 221

Our brain contains about **1011** neurons, each of which is connected to an average of **104** other neurons. This amounts to a total of **1015** connections. If these billions of connections were fully random, it can be shown that the brain would be many times larger than it actually is (Happel and Murre, 1994). Massive regressive events of neuronal connectivity in the vertebrate nervous system can be seen as part of the development and maturation of neural functions. The neuron has set of nodes that, connects it to inputs, output, or other neurons,

A single neuron by itself is not a very useful pattern recognition tool. The real power of neural networks comes when we combine neurons into the multilayer structures, called

Genetic algorithms (GAs) are currently the most prominent and widely used models of evolution in artificial-life systems. GAs have been used both as tools for solving practical problems and as scientific models of evolutionary processes. The intersection between GAs and artificial life includes both, although in this article we focus primarily on GAs as models of natural phenomena. Indeed GAs are optimization algorithms that work according to a scheme analogous to that of natural evolution. Literature review reveled that John Holland (Holland 1975) was the first who offered these principles of natural evolution to artificial systems, more precisely to optimization problems, and came up with the notion of GA. A general definition of these algorithms is (Koza 1980):

*"The genetic algorithm is a highly parallel mathematical algorithm that transforms a set (population) of individual mathematical objects (typically fixed length character strings patterned after chromosome strings), each with an associated fitness value, into a new population (i.e. the next generation) using operations patterned after the Darwinian principle of reproduction and survival of the fittest and after naturally occurring genetic operations (notably sexual recombination)."* 

Genetic algorithms as defined by Goldberg (Goldberg, 1989) is:

...search algorithms based on the mechanics of natural selection and natural genetics."

Goldberg offers four differences between genetic algorithms and other search methods.


The success of nature in solving many problems nowadays recognized as very difficult for the traditional approaches, have led researchers into studying the biological example. In various abstractions and formalizations, biological systems have been theoretically proven to provide robust solutions to these hard problems. However the models used in the area of biological problems are complex, because of their characteristics and processes. This concept leads to the conclusion that the biological activity generates information with special features, most notable being the following (Fernández and Lozano, 2010):


What structure we need depends, of course, on our aims. We may distinguish roughly between operational and physiological models. An operational model aims to describe behavior realistically, but its structure is not intended to resemble the internal structure of particular biological system. Such models are often referred to as black box models to indicate lack of concern about underlying mechanisms. A physiological model, on the other hand, attempts to take into account more of the physiology that produces behavior, e.g., body and nervous system physiology (Dellaert, 1995 ).

### **3. Neural networks**

220 Real-World Applications of Genetic Algorithms

Genetic algorithms (GAs) are currently the most prominent and widely used models of evolution in artificial-life systems. GAs have been used both as tools for solving practical problems and as scientific models of evolutionary processes. The intersection between GAs and artificial life includes both, although in this article we focus primarily on GAs as models of natural phenomena. Indeed GAs are optimization algorithms that work according to a scheme analogous to that of natural evolution. Literature review reveled that John Holland (Holland 1975) was the first who offered these principles of natural evolution to artificial systems, more precisely to optimization problems, and came up with the notion of GA. A

*"The genetic algorithm is a highly parallel mathematical algorithm that transforms a set (population) of individual mathematical objects (typically fixed length character strings patterned after chromosome strings), each with an associated fitness value, into a new population (i.e. the next generation) using operations patterned after the Darwinian principle of reproduction and survival of the fittest and after naturally occurring genetic operations (notably* 

...search algorithms based on the mechanics of natural selection and natural genetics." Goldberg offers four differences between genetic algorithms and other search methods.

3. They only use directly available information provided through a *fitness function*.

4. They rely on *probabilistic transition rules* instead of deterministic rules.

features, most notable being the following (Fernández and Lozano, 2010):

2. They search from a *population* of points in a solution space, rather than from a single

The success of nature in solving many problems nowadays recognized as very difficult for the traditional approaches, have led researchers into studying the biological example. In various abstractions and formalizations, biological systems have been theoretically proven to provide robust solutions to these hard problems. However the models used in the area of biological problems are complex, because of their characteristics and processes. This concept leads to the conclusion that the biological activity generates information with special

1. The obtained information from process presents a non-homogeneous structure since of

2. The information is emerging from the dynamics of change associated with the

What structure we need depends, of course, on our aims. We may distinguish roughly between operational and physiological models. An operational model aims to describe behavior realistically, but its structure is not intended to resemble the internal structure of particular biological system. Such models are often referred to as black box models to indicate lack of concern about underlying mechanisms. A physiological model, on the other hand, attempts to take into account more of the physiology that produces behavior, e.g.,

general definition of these algorithms is (Koza 1980):

Genetic algorithms as defined by Goldberg (Goldberg, 1989) is:

1. Genetic algorithms work with a *coded* parameter set.

the complexity of the objects alive.

functional properties of the studied phenomena.

body and nervous system physiology (Dellaert, 1995 ).

*sexual recombination)."* 

point.

Our brain contains about **1011** neurons, each of which is connected to an average of **104** other neurons. This amounts to a total of **1015** connections. If these billions of connections were fully random, it can be shown that the brain would be many times larger than it actually is (Happel and Murre, 1994). Massive regressive events of neuronal connectivity in the vertebrate nervous system can be seen as part of the development and maturation of neural functions. The neuron has set of nodes that, connects it to inputs, output, or other neurons, these nodes are also called synapses (See Fig. 1).

Fig. 1. Schematic structure of a Neuron

A single neuron by itself is not a very useful pattern recognition tool. The real power of neural networks comes when we combine neurons into the multilayer structures, called neural networks (NN) (Fig. 2).

Fig. 2. A Schematic representation of neural net

Genetic Algorithms and Group Method

(Madala and Ivakhnenko, 1994).

identification.

perceptron type".

of Data Handling-Type Neural Networks Applications in Poultry Science 223

variants. The framework of these methods differs slightly in some important respects

A major difficulty in modeling complex systems in such unstructured areas as economics, ecology, sociology, and others is the problem of the researchers introducing their own prejudices model. In the mid 1960's the Russian mathematician and cyberneticist, A.G. Ivakhnenko, introduced a method (Ivakhnenko, 1966), based in part on the Rosenblatt perceptron (Rosenblatt, 1958), that allows the researchers to build model of complex systems without making assumptions about internal working. The idea is to have the computer construct a model of optimal complexity based only on data and not on any preconceived ideas of the researchers; that is,by knowing only simple input-output relationship of the system. Ivakhnenko's GMDH algorithm will construct a self-organizing model (an extremely high-order polynomial in the input variable) that can be used to solve prediction,

The algorithm was developed for identifying nonlinear relationships between inputs and outputs. The algorithm provides an optimal structure, obtained in an iterative procedure of partial descriptions of the data by adding new layers. The number of neurons in each layer, the number of layers and the input variables are automatically determined to minimize a criterion of prediction error and thus organizes an optimal NN architecture using a selfheuristics, which is the basis of the GMDH algorithm. (Ivakhnenko,1971). This method is particularly successful in solving problems of modeling multiple entries for a single output (Mutasem, 2004). The main idea of GMDH is to build an analytical function in a feedforward network based on a quadratic node transfer function whose coefficients obtained by using a regression technique (Farlow, 1984). By means of the GMDH-type NN algorithm, a model can be represented as a set of neurons in which different pairs in each layer are connected through a quadratic polynomial and thus produce new neurons in the next layer, and therefore can be used to map inputs to outputs. Such an NN identification process needs some optimization method to find the best network architecture. This sub-model of ANN is considered as a Self - organizing approach by which gradually more complex models are generated from their performance evaluation (Lemke and Mueller, 2003). The unique feature of GMDH-type NN is that it facilitates, systematically and autonomously, developing optimal complex models by performing both variable and structure

Incorporating Genetic Algorithm to GMDH-type NNs, each neuron is represented as a string, which can be mutated or crossed with each other to form new generations. Thus GA has been used in feed-forward GMDH-type NN for each neuron searching its optimal set of connections with the preceding layer (Vasechkina & Yarin 2001; Nariman-Zadeh et al. 2003). In the early stage of the development of GMDH theory the similarity between NNs and multilayer GMDH algorithms had been highlighted. (Ivakhnenko 1970) in one of the introductory articles claims that since the differences between perceptron and GMDH are neither significant nor fundamental it is appropriate to call GMDH systems as "systems of

During the modeling procedure, GMDH algorithm involves four heuristics that represent

the main features of GMDH theory (Anastasakis and Mort, 2001):

identification, control synthesis, and other system problems (Farrow, 1981).

In a nutshell, a NN can be considered as a black box that is able to predict an output pattern when it recognizes a given input pattern. Once trained, the NN is able to recognize similarities when presented with a new input pattern, resulting in a predicted output pattern. The process of evolution is used as a 'real-world model' that serves as a source of ideas for solving practical and theoretical problems in modeling and optimization.

Researchers from a wide range of fields have discovered the benefits of applying NNs to pattern recognition problems in various systems (including biological system). Artificial NN (ANNs) is a system loosely modeled based on the human brain and are considered as a branch of the field known as "Artificial Intelligence" (AI). The techniques of AI is being applied in this field significantly in recent decades, and among them those known as ANNs, are characterized by their properties of learning and generalization. It's often necessary to take into account their potential for induction, which can be implemented by software (Miroslav Šnorek, 2006).

To understand behavior of model system, we need ways of describing behavior maps and state transition equations. Ideally, behavioral models should fulfill the following requirements (Dellaert, 1995):


NNs are a powerful technique to solve many real world problems. They have the ability to learn from experience in order to improve their performance and to adapt themselves to changes in the environment. In addition to that they are able to deal with incomplete information or noisy data and can be very effective especially in situations where it is not possible to define the rules or steps that lead to the solution of a problem. Once trained, the NN is able to recognize similarities when presented with a new input pattern, resulting in a predicted output pattern. NNs are applied in many fields to model and predict the behavior of unknown systems or systems with complexity (or both) based on given input–output data. Using NNs does not require a priori equation or model. This characteristic is potentially advantageous in modeling biological processes (Dayhof & DeLeo 2001).

There are several methods to obtain inductive models. The Group Method of Data Handling methods, GMDH (Ivakhnenko AG, 1971) is well known and has recently gained popularity as a self-organizing and powerful tool to express complex input–output dependencies.

### **4. Group Method of Data Handling method (GMDH)**

Generally, the connection between input-output variables can be approximated by Volterra functional series, the discrete analogue of which is Kolmogorov-Gabor polynomial. Ivakhnenko (Ivakhnenko, 1966), inspired by the form of Kolmogorov-Gabor polynomial, developed a new algorithm, known as Group Method of Data Handling (GMDH), which are also called inductive learning methods, self-organization, sorting out, and heuristic methods. This approach is substantially different from deductive methods used commonly for modeling. It has inductive nature, i.e., it finds the best solution by sorting-out of possible

In a nutshell, a NN can be considered as a black box that is able to predict an output pattern when it recognizes a given input pattern. Once trained, the NN is able to recognize similarities when presented with a new input pattern, resulting in a predicted output pattern. The process of evolution is used as a 'real-world model' that serves as a source of

Researchers from a wide range of fields have discovered the benefits of applying NNs to pattern recognition problems in various systems (including biological system). Artificial NN (ANNs) is a system loosely modeled based on the human brain and are considered as a branch of the field known as "Artificial Intelligence" (AI). The techniques of AI is being applied in this field significantly in recent decades, and among them those known as ANNs, are characterized by their properties of learning and generalization. It's often necessary to take into account their potential for induction, which can be implemented by software

To understand behavior of model system, we need ways of describing behavior maps and state transition equations. Ideally, behavioral models should fulfill the following

NNs are a powerful technique to solve many real world problems. They have the ability to learn from experience in order to improve their performance and to adapt themselves to changes in the environment. In addition to that they are able to deal with incomplete information or noisy data and can be very effective especially in situations where it is not possible to define the rules or steps that lead to the solution of a problem. Once trained, the NN is able to recognize similarities when presented with a new input pattern, resulting in a predicted output pattern. NNs are applied in many fields to model and predict the behavior of unknown systems or systems with complexity (or both) based on given input–output data. Using NNs does not require a priori equation or model. This characteristic is

potentially advantageous in modeling biological processes (Dayhof & DeLeo 2001).

**4. Group Method of Data Handling method (GMDH)**

There are several methods to obtain inductive models. The Group Method of Data Handling methods, GMDH (Ivakhnenko AG, 1971) is well known and has recently gained popularity as a self-organizing and powerful tool to express complex input–output dependencies.

Generally, the connection between input-output variables can be approximated by Volterra functional series, the discrete analogue of which is Kolmogorov-Gabor polynomial. Ivakhnenko (Ivakhnenko, 1966), inspired by the form of Kolmogorov-Gabor polynomial, developed a new algorithm, known as Group Method of Data Handling (GMDH), which are also called inductive learning methods, self-organization, sorting out, and heuristic methods. This approach is substantially different from deductive methods used commonly for modeling. It has inductive nature, i.e., it finds the best solution by sorting-out of possible

ideas for solving practical and theoretical problems in modeling and optimization.

(Miroslav Šnorek, 2006).

1. Versatility 2. Robustness 3. Learning 4. Ontogeny 5. Evolution

requirements (Dellaert, 1995):

variants. The framework of these methods differs slightly in some important respects (Madala and Ivakhnenko, 1994).

A major difficulty in modeling complex systems in such unstructured areas as economics, ecology, sociology, and others is the problem of the researchers introducing their own prejudices model. In the mid 1960's the Russian mathematician and cyberneticist, A.G. Ivakhnenko, introduced a method (Ivakhnenko, 1966), based in part on the Rosenblatt perceptron (Rosenblatt, 1958), that allows the researchers to build model of complex systems without making assumptions about internal working. The idea is to have the computer construct a model of optimal complexity based only on data and not on any preconceived ideas of the researchers; that is,by knowing only simple input-output relationship of the system. Ivakhnenko's GMDH algorithm will construct a self-organizing model (an extremely high-order polynomial in the input variable) that can be used to solve prediction, identification, control synthesis, and other system problems (Farrow, 1981).

The algorithm was developed for identifying nonlinear relationships between inputs and outputs. The algorithm provides an optimal structure, obtained in an iterative procedure of partial descriptions of the data by adding new layers. The number of neurons in each layer, the number of layers and the input variables are automatically determined to minimize a criterion of prediction error and thus organizes an optimal NN architecture using a selfheuristics, which is the basis of the GMDH algorithm. (Ivakhnenko,1971). This method is particularly successful in solving problems of modeling multiple entries for a single output (Mutasem, 2004). The main idea of GMDH is to build an analytical function in a feedforward network based on a quadratic node transfer function whose coefficients obtained by using a regression technique (Farlow, 1984). By means of the GMDH-type NN algorithm, a model can be represented as a set of neurons in which different pairs in each layer are connected through a quadratic polynomial and thus produce new neurons in the next layer, and therefore can be used to map inputs to outputs. Such an NN identification process needs some optimization method to find the best network architecture. This sub-model of ANN is considered as a Self - organizing approach by which gradually more complex models are generated from their performance evaluation (Lemke and Mueller, 2003). The unique feature of GMDH-type NN is that it facilitates, systematically and autonomously, developing optimal complex models by performing both variable and structure identification.

Incorporating Genetic Algorithm to GMDH-type NNs, each neuron is represented as a string, which can be mutated or crossed with each other to form new generations. Thus GA has been used in feed-forward GMDH-type NN for each neuron searching its optimal set of connections with the preceding layer (Vasechkina & Yarin 2001; Nariman-Zadeh et al. 2003).

 In the early stage of the development of GMDH theory the similarity between NNs and multilayer GMDH algorithms had been highlighted. (Ivakhnenko 1970) in one of the introductory articles claims that since the differences between perceptron and GMDH are neither significant nor fundamental it is appropriate to call GMDH systems as "systems of perceptron type".

During the modeling procedure, GMDH algorithm involves four heuristics that represent the main features of GMDH theory (Anastasakis and Mort, 2001):

Genetic Algorithms and Group Method

identification and control problems.

conclusion remarks of this study were reported as:

nutritional factors by optimizing methods such as the GA.

appropriate for the TMEn prediction of a feedstuff.

meal and poultry offal meal (POM) based on their chemical analysis.

production costs.

inputs.

of Data Handling-Type Neural Networks Applications in Poultry Science 225

methods is ANNs, which have shown great ability in solving complex nonlinear system

The optimal structures of the evolved 2- hidden-layer GMDH-type NN that were suggested by GA for performance index (PI) as the system output modeling, were found with 2, 4, and 4 hidden neurons for growth periods 1, 2, and 3, respectively. In the first period, the structure obtained appeared with the GA, which was less complex than in the second and third periods, in which the GA suggested 2 hidden neurons to fit the network. All models constructed from this data set were characterized by a superb response for all input variables from the learning set. The partial descriptions of the GMDH-type NN were found with 2 hidden layers and 2 hidden neurons for growth period, whereas it appeared with 2 hidden layers and 4 hidden neurons for growth periods 2 and 3. In fact, these results revealed the quantitative relation between input (ME, Met, and Lys) and output (PI) variables under investigation, which meant GMDH-type NN may be considered as a promising method for modeling the relationship between dietary concentrations of nutrients and poultry performance, and therefore can be used in choosing and developing special feeding programs to decrease production costs. Also, it can enhance our ability to predict other economic traits, make precise predictions of the nutrition requirements, and achieve optimal performance in poultry production systems. The

4-1- Knowledge of an adequate description of broiler ME and AA requirements can help in establishing specific feeding programs, defining optimal performance, and reducing

4-2- Calculated statistics indicate that GMDH-type NN provide an effective means of efficiently recognizing the patterns in data and predicting a PI based on investigating

4-3- The genetic approach could be used to provide optimal networks in terms of hidden layers, the number of neurons and their configuration of connectivity, or both so that a polynomial expression for dependent variables of the process can consequently be achieved. 4-4-The polynomials obtained could be used to optimize broiler performance based on

In animal and poultry production, feed composition is very important item for diet formulation. Since conventional laboratory techniques for feed analysis is expensive and time consuming, it would be advantageous if a simple means of estimating feed composition could be developed. One year after the first work another study (Ahmadi *etal.,* 2008) was done. The purpose of this study was to examine the validity of GMDH-type NN with a genetic algorithm method to predict the True Metabolizable Energy corrected for nitrogen (TMEn) of feather

All the previously TMEn prediction models reported for poultry by-product meals were based on the regression analysis methods using their CP, ether extract (EE), and ash content. In this study, a soft-computing method of artificial NN (ANN) seemed to be more

The parameters of interest in this multi-input, single-output system that influenced the TMEn were CP, EE, and ash content of the samples. The raw data were divided into 2 parts of training and validation sets. Thirty input-output data lines (12 from FM and 18 from POM


A detailed description of a GMDH-type NN terminology, development, application, and examples of using this approach were reported by several researchers (Farrow, 1984; Mueller and Lemke, 2000; Lemke and Mueller, 2003; Nariman-Zadeh et al., 2005). Recently the GEvoM software for GMDH-type NN training (GEvoM 2009) was developed in University of Guilan, Iran.

### **5. Applications of GMDH-type algorithms in animal and poultry production systems**

Contributions to GMDH type of NN, have come from many research areas of different disciplines, and recently, the use of such self-organizing networks has led to a successful application of the GMDH-type algorithm in a broad range of areas in engineering, science, and economics (Amanifard et al., 2008). However, very little research has been conducted on modeling animal and /or poultry growth and production using ANNs.

A series of studies have been conducted to examine the potential use of ANNs in various poultry subjects, such, prediction of ascites in broilers (Roush *et al.,* 1996; Roush and Wideman,2000), the estimation of production variables in the production phase of broiler breeders (Salle *et al.,* 2003), and the comparison of Gompertz and NN models of broiler growth (Roush *et al.,* 2006).

However no attempt was made to use GMDH-type NN in animal agriculture, until 2007, when the results of study was published based on the first work of my group in University of Guilan, Iran (Ahmadi, *etal.,* 2007). The idea behind this work was that, when considering the effects of nutrition on broiler performance, several nutrients may influence the breast meat yield, feed : gain ratio, and number of days required to produce the market body weight; among them, Metabolizable Energy (ME) and Amino Acid(AA) , such as Lysine(Lys) and Methionine (Met) (Hruby and Hamre, 1996 ; Gous, 1998) . In terms of AA, whatever system is used to describe the essential AA requirements for broiler chickens, predicting the performance to be used in deciding the most advantageous dietary AA patterns in practical and useful terms is still difficult, even when the digestibility or availability of AA is specified (NRC, 1994; Sibbald, 1987 ]. This difficulty is partly due to the nonlinearity of growth responses related to changes in dietary AA concentrations [Hruby and Hamre, 1996, Phillips, 1981]. A more useful method is to model the system, which in turn requires an explicit mathematical input-output relationship. Such explicit mathematical modeling is, however, very difficult and is not readily tractable in poorly understood systems. Alternatively, soft-computing methods, which concern computation in an imprecise environment, have gained significant attention. One of the soft-computing

2. Divide the observations into two groups. The first will be used to estimate the coefficients of model while the second will separate the information embedded in the data into either useful or harmful. Strictly speaking: "no partition of the data, no

3. Create a set of elementary functions where complexity will increase through an iterative

4. Acording to Gödel's incompleteness theorem, apply an external criterion to choose the

A detailed description of a GMDH-type NN terminology, development, application, and examples of using this approach were reported by several researchers (Farrow, 1984; Mueller and Lemke, 2000; Lemke and Mueller, 2003; Nariman-Zadeh et al., 2005). Recently the GEvoM software for GMDH-type NN training (GEvoM 2009) was developed in

**5. Applications of GMDH-type algorithms in animal and poultry production** 

modeling animal and /or poultry growth and production using ANNs.

Contributions to GMDH type of NN, have come from many research areas of different disciplines, and recently, the use of such self-organizing networks has led to a successful application of the GMDH-type algorithm in a broad range of areas in engineering, science, and economics (Amanifard et al., 2008). However, very little research has been conducted on

A series of studies have been conducted to examine the potential use of ANNs in various poultry subjects, such, prediction of ascites in broilers (Roush *et al.,* 1996; Roush and Wideman,2000), the estimation of production variables in the production phase of broiler breeders (Salle *et al.,* 2003), and the comparison of Gompertz and NN models of broiler

However no attempt was made to use GMDH-type NN in animal agriculture, until 2007, when the results of study was published based on the first work of my group in University of Guilan, Iran (Ahmadi, *etal.,* 2007). The idea behind this work was that, when considering the effects of nutrition on broiler performance, several nutrients may influence the breast meat yield, feed : gain ratio, and number of days required to produce the market body weight; among them, Metabolizable Energy (ME) and Amino Acid(AA) , such as Lysine(Lys) and Methionine (Met) (Hruby and Hamre, 1996 ; Gous, 1998) . In terms of AA, whatever system is used to describe the essential AA requirements for broiler chickens, predicting the performance to be used in deciding the most advantageous dietary AA patterns in practical and useful terms is still difficult, even when the digestibility or availability of AA is specified (NRC, 1994; Sibbald, 1987 ]. This difficulty is partly due to the nonlinearity of growth responses related to changes in dietary AA concentrations [Hruby and Hamre, 1996, Phillips, 1981]. A more useful method is to model the system, which in turn requires an explicit mathematical input-output relationship. Such explicit mathematical modeling is, however, very difficult and is not readily tractable in poorly understood systems. Alternatively, soft-computing methods, which concern computation in an imprecise environment, have gained significant attention. One of the soft-computing

1. Collect a set of observations that seems to be relevant to the object

procedure producing different models.

GMDH".

optimum model.

University of Guilan, Iran.

growth (Roush *et al.,* 2006).

**systems** 

methods is ANNs, which have shown great ability in solving complex nonlinear system identification and control problems.

The optimal structures of the evolved 2- hidden-layer GMDH-type NN that were suggested by GA for performance index (PI) as the system output modeling, were found with 2, 4, and 4 hidden neurons for growth periods 1, 2, and 3, respectively. In the first period, the structure obtained appeared with the GA, which was less complex than in the second and third periods, in which the GA suggested 2 hidden neurons to fit the network. All models constructed from this data set were characterized by a superb response for all input variables from the learning set. The partial descriptions of the GMDH-type NN were found with 2 hidden layers and 2 hidden neurons for growth period, whereas it appeared with 2 hidden layers and 4 hidden neurons for growth periods 2 and 3. In fact, these results revealed the quantitative relation between input (ME, Met, and Lys) and output (PI) variables under investigation, which meant GMDH-type NN may be considered as a promising method for modeling the relationship between dietary concentrations of nutrients and poultry performance, and therefore can be used in choosing and developing special feeding programs to decrease production costs. Also, it can enhance our ability to predict other economic traits, make precise predictions of the nutrition requirements, and achieve optimal performance in poultry production systems. The conclusion remarks of this study were reported as:

4-1- Knowledge of an adequate description of broiler ME and AA requirements can help in establishing specific feeding programs, defining optimal performance, and reducing production costs.

4-2- Calculated statistics indicate that GMDH-type NN provide an effective means of efficiently recognizing the patterns in data and predicting a PI based on investigating inputs.

4-3- The genetic approach could be used to provide optimal networks in terms of hidden layers, the number of neurons and their configuration of connectivity, or both so that a polynomial expression for dependent variables of the process can consequently be achieved.

4-4-The polynomials obtained could be used to optimize broiler performance based on nutritional factors by optimizing methods such as the GA.

In animal and poultry production, feed composition is very important item for diet formulation. Since conventional laboratory techniques for feed analysis is expensive and time consuming, it would be advantageous if a simple means of estimating feed composition could be developed. One year after the first work another study (Ahmadi *etal.,* 2008) was done. The purpose of this study was to examine the validity of GMDH-type NN with a genetic algorithm method to predict the True Metabolizable Energy corrected for nitrogen (TMEn) of feather meal and poultry offal meal (POM) based on their chemical analysis.

All the previously TMEn prediction models reported for poultry by-product meals were based on the regression analysis methods using their CP, ether extract (EE), and ash content. In this study, a soft-computing method of artificial NN (ANN) seemed to be more appropriate for the TMEn prediction of a feedstuff.

The parameters of interest in this multi-input, single-output system that influenced the TMEn were CP, EE, and ash content of the samples. The raw data were divided into 2 parts of training and validation sets. Thirty input-output data lines (12 from FM and 18 from POM

Genetic Algorithms and Group Method

of appropriate coefficients of quadratic expressions.

in nutrition and management (Havenstein et al., 2007).

protein versus fat accretion, parts yield, and feed intake.

a strong effect of all selected input variables on output prediction.

of Data Handling-Type Neural Networks Applications in Poultry Science 227

The success of poultry meat production has been strongly related to improvements in growth and carcass yield, mainly by increasing breast proportion and reducing carcass fat. In addition to its measurement in the laboratory using wet chemistry, carcass composition of broiler chickens has been predicted by means of allometric equations, real-time ultrasonography, radioactive isotopes and specific gravity studies (Pesti & Bakalli 1997; Toledo et al.2004; Rosa et al.2007; Makkar 2008). Conventional laboratory techniques for determining carcass composition are expensive, cumbersome and time consuming. Therefore, it would be useful if a simple means of estimating carcass composition could be developed. In this respect, the potential advantages from modeling growth are considerable. Results obtained from two above mentioned studies urged our group to think about third study in 2010 (Faridi, *etal.*, 2012), aimed at applying the GMDH-type NNs to data from two studies with broilers in order to predict carcass energy (CEn, MJ/g) content and relative growth (g/g of body weight) of carcass components (carcass protein, breast muscle, leg and thigh muscles, carcass fat, abdominal fat, skin fat and visceral fat). The effective input variables involved in the prediction of CEn and carcass fat content using data from the first study were dietary metabolizable energy (ME, kJ/kg), crude protein (CP, g/kg of diet), fat (g/kg of diet) and crude fibre (CF, g/kg of diet). For this purpose, in the current study, GA were deployed to design the whole architecture of the GMDH-type NN, i.e. the number of neurons in each hidden layer and their configuration of connectivity to find the optimal set

Quantified values of bias in this study showed very little under- and over-estimation by the models proposed by the GMDH-type NN, which revealed close agreement between observed and predicted values of CEn and carcass components. The value of R2, a measure of the relation between the actual and predicted values, was high for both studies indicating

In conclusion, the results of the current study showed that a GMDH type NN modeling approach can be a simple but very effective method for predicting carcass composition in broiler chickens based on dietary input variables. This is in agreement with previous studies aimed at investigating the effects of different dietary nutrients on body composition in

Selection pressure applied by industry geneticists has greatly reduced feed conversion ratio and age to slaughter as well as increased growth rate and yield of edible meat for commercial turkeys. These genetic improvements have occurred along with improvements

There has been extensive research conducted to clarify protein, essential amino acids, and energy requirements in poultry. To avoid conventional laboratory and field based techniques problems for determining nutrient requirements alternative methods was offered using GMDH – type NN (Mottaghitalab,*etal.,* 2010). In determining nutrient requirements, the potential benefits from modeling growth in poultry are considerable. This approach has the potential to provide information in several areas for poultry production, including prediction of growth rate and market weights, determination of factors that are truly of economic importance to the operation, general knowledge about the systems involved in production, and determination of more precise nutrient requirements based on sex, strain,

broilers (Fraps 1943; Donaldson et al. 1956; Kubena et al. 1972; Edwards et al. 1973).

Samples ) were randomly selected and used to train the GMDH-type NN model as a training set. The validation set consisted of the 7 remaining data lines (3 from FM and 4 from POM samples), which were used to validate the prediction of the evolved NN during the training processes. The data set was imported into a GEvoM for GMDH-type NN training (GEvoM, 2008). Two hidden layers were considered for prediction of the TMEn model. A population of 15 individual values with a crossover probability of 0.7, mutation probability of

0.07, and 300 generations was used to genetically design the NN (Yao, 1999). It appeared that no further improvement could be achieved for this population size. A quantitative verifying fit for the predictive model was made using error measurement indices commonly used to evaluate forecasting models. The goodness of fit or accuracy of the model was determined by R2 value, adjusted R2, mean square error, residual standard deviation, mean absolute percentage error, and bias (Oberstone,1990).

The results of this study revealed that the novel modeling of GMDH-type NN with an evolutionary method of GA can be used to predict the TMEn of FM and POM samples based on their CP, EE, and ash content( See Fig. 3). The advantage of using the GMDH-type NN to predict an output from the input variables is that there is no need to preselect a model or base the model entirely on the fit of the data. It is concluded that the GMDH-type NN may be used to accurately estimate the nutritive value of poultry meals from their corresponding chemical composition.

Fig. 3. The comparison of observed and model predicted TMEn values obtained from training (1 to 12 and 13 to 30 are feather and offal samples, respectively) and validation (31 to 33 and 34 to 37 are feather and offal samples, respectively) sets.

Samples ) were randomly selected and used to train the GMDH-type NN model as a training set. The validation set consisted of the 7 remaining data lines (3 from FM and 4 from POM samples), which were used to validate the prediction of the evolved NN during the training processes. The data set was imported into a GEvoM for GMDH-type NN training (GEvoM, 2008). Two hidden layers were considered for prediction of the TMEn model. A population of 15 individual values with a crossover probability of 0.7, mutation

0.07, and 300 generations was used to genetically design the NN (Yao, 1999). It appeared that no further improvement could be achieved for this population size. A quantitative verifying fit for the predictive model was made using error measurement indices commonly used to evaluate forecasting models. The goodness of fit or accuracy of the model was determined by R2 value, adjusted R2, mean square error, residual standard deviation, mean

The results of this study revealed that the novel modeling of GMDH-type NN with an evolutionary method of GA can be used to predict the TMEn of FM and POM samples based on their CP, EE, and ash content( See Fig. 3). The advantage of using the GMDH-type NN to predict an output from the input variables is that there is no need to preselect a model or base the model entirely on the fit of the data. It is concluded that the GMDH-type NN may be used to accurately estimate the nutritive value of poultry meals from their corresponding

Fig. 3. The comparison of observed and model predicted TMEn values obtained from training (1 to 12 and 13 to 30 are feather and offal samples, respectively) and validation

(31 to 33 and 34 to 37 are feather and offal samples, respectively) sets.

absolute percentage error, and bias (Oberstone,1990).

probability of

chemical composition.

The success of poultry meat production has been strongly related to improvements in growth and carcass yield, mainly by increasing breast proportion and reducing carcass fat.

In addition to its measurement in the laboratory using wet chemistry, carcass composition of broiler chickens has been predicted by means of allometric equations, real-time ultrasonography, radioactive isotopes and specific gravity studies (Pesti & Bakalli 1997; Toledo et al.2004; Rosa et al.2007; Makkar 2008). Conventional laboratory techniques for determining carcass composition are expensive, cumbersome and time consuming. Therefore, it would be useful if a simple means of estimating carcass composition could be developed. In this respect, the potential advantages from modeling growth are considerable.

Results obtained from two above mentioned studies urged our group to think about third study in 2010 (Faridi, *etal.*, 2012), aimed at applying the GMDH-type NNs to data from two studies with broilers in order to predict carcass energy (CEn, MJ/g) content and relative growth (g/g of body weight) of carcass components (carcass protein, breast muscle, leg and thigh muscles, carcass fat, abdominal fat, skin fat and visceral fat). The effective input variables involved in the prediction of CEn and carcass fat content using data from the first study were dietary metabolizable energy (ME, kJ/kg), crude protein (CP, g/kg of diet), fat (g/kg of diet) and crude fibre (CF, g/kg of diet). For this purpose, in the current study, GA were deployed to design the whole architecture of the GMDH-type NN, i.e. the number of neurons in each hidden layer and their configuration of connectivity to find the optimal set of appropriate coefficients of quadratic expressions.

Quantified values of bias in this study showed very little under- and over-estimation by the models proposed by the GMDH-type NN, which revealed close agreement between observed and predicted values of CEn and carcass components. The value of R2, a measure of the relation between the actual and predicted values, was high for both studies indicating a strong effect of all selected input variables on output prediction.

In conclusion, the results of the current study showed that a GMDH type NN modeling approach can be a simple but very effective method for predicting carcass composition in broiler chickens based on dietary input variables. This is in agreement with previous studies aimed at investigating the effects of different dietary nutrients on body composition in broilers (Fraps 1943; Donaldson et al. 1956; Kubena et al. 1972; Edwards et al. 1973).

Selection pressure applied by industry geneticists has greatly reduced feed conversion ratio and age to slaughter as well as increased growth rate and yield of edible meat for commercial turkeys. These genetic improvements have occurred along with improvements in nutrition and management (Havenstein et al., 2007).

There has been extensive research conducted to clarify protein, essential amino acids, and energy requirements in poultry. To avoid conventional laboratory and field based techniques problems for determining nutrient requirements alternative methods was offered using GMDH – type NN (Mottaghitalab,*etal.,* 2010). In determining nutrient requirements, the potential benefits from modeling growth in poultry are considerable. This approach has the potential to provide information in several areas for poultry production, including prediction of growth rate and market weights, determination of factors that are truly of economic importance to the operation, general knowledge about the systems involved in production, and determination of more precise nutrient requirements based on sex, strain, protein versus fat accretion, parts yield, and feed intake.

Genetic Algorithms and Group Method

*nutrient intake" (Faridi et al., 2011).* 

training set.

of Data Handling-Type Neural Networks Applications in Poultry Science 229

Zeuthen et al., 2008). The calculated values of CE model error measurement showed that the testing set for toms yielded lower values of MS error, mean absolute deviation, mean absolute percentage error, mean relative error, and higher values of R2 compared with the

Fig. 5. Evolved structure of the generalized group method of data handling-type NNs for feed efficiency in tom turkeys. The letters *a, b, c, d*, and *e* stand for input variables of age (wk), ME (kcal/g), CP (% of diet), Met (% of diet), and Lys (% of diet), respectively. This figure illustrates the generated relationship between input variables to reach output.

Conducting a sensitivity analysis (SA) on the obtained polynomial equations reveals the sensitivity of model output to input variables. Hence it is necessary to do sensitivity analysis for any proposed model. In other words, SA increases confidence in the model and its predictions by providing an understanding of how the model responds to changes in its inputs. Moreover, the SA identifies critical regions in the space of the inputs, establishes priorities for research and simplifies the model (Castillo et al.,2008; Saltelli et al., 2008).

*"Sensitivity analysis of an early egg production predictive model inbroiler breeders based on dietary* 

Although the use of NN and SA techniques has led to successful application in a broad range of areas (Seyedan&Ching 2006;Lee &Hsiung 2009), the use of SA along with NN models is appeared uncommon in poultry science. The aim of the present study was to use the GMDH-type NN to model early egg production (EEP) in broiler breeders (BB) based on the dietary intake levels of ME, CP, and the two first limiting amino acids, methionine (Met)

For such reason and in line with previous work, another study was designed, titled:

The structures of the 2 hidden layers GMDH-type NN evolved for CE and FE are shown in Figures 4 and 5, respectively. These figures correspond to the genome representation of (abceadaa) and (eeabacdd) for the CE and FE models, respectively, and illustrate the generated relationship between input variables to reach the output. As Figures 4 and 5 show, the optimal structure of the evolved 2 hidden layer GMDH-type NN suggested by GA was found with 5 and 4 hidden neurons for CE and FE, respectively. In most GMDHtype NN, the neurons in each layer are only connected to neurons in the adjacent layer (Farlow, 1984), but for GMDH-type NN developed here, variable a of the input layer for CE is connected to adaa in the second hidden layer by directly passing through the first hidden layer. The same process happens for d and e input variables in the FE model. Such repetition occurs whenever a neuron passes some adjacent hidden layer and connects to another neuron in the next following hidden layer.

Fig. 4. Evolved structure of the generalized group method of data handling-type NNs for caloric efficiency in tom turkeys. The letters *a, b, c, d*, and *e* stand for input variables of age (wk), ME (kcal/g), CP (% of diet), Met (% of diet), and Lys (% of diet), respectively. This figure illustrates the generated relationship between input variables to reach output.

It appears that all selected input variables in both models had a strong effect on output prediction, which is in agreement with previous studies (Lemme et al., 2006 for amino acid; Noy and Sklan, 2004 for energy and amino acid; Potter et al., 1966 and Waibel et al., 1995 for Met and Lys; and Bowyer and Waldroup, 1986 for protein). Figure 1 shows a very strong effect of age on CE. This result is similar to previous studies aiming to describe the growth pattern of animals with age using growth functions (Darmani-Kuhi et al., 2003; Schulin-

The structures of the 2 hidden layers GMDH-type NN evolved for CE and FE are shown in Figures 4 and 5, respectively. These figures correspond to the genome representation of (abceadaa) and (eeabacdd) for the CE and FE models, respectively, and illustrate the generated relationship between input variables to reach the output. As Figures 4 and 5 show, the optimal structure of the evolved 2 hidden layer GMDH-type NN suggested by GA was found with 5 and 4 hidden neurons for CE and FE, respectively. In most GMDHtype NN, the neurons in each layer are only connected to neurons in the adjacent layer (Farlow, 1984), but for GMDH-type NN developed here, variable a of the input layer for CE is connected to adaa in the second hidden layer by directly passing through the first hidden layer. The same process happens for d and e input variables in the FE model. Such repetition occurs whenever a neuron passes some adjacent hidden layer and connects to another

Fig. 4. Evolved structure of the generalized group method of data handling-type NNs for caloric efficiency in tom turkeys. The letters *a, b, c, d*, and *e* stand for input variables of age (wk), ME (kcal/g), CP (% of diet), Met (% of diet), and Lys (% of diet), respectively. This figure illustrates the generated relationship between input variables to reach output.

It appears that all selected input variables in both models had a strong effect on output prediction, which is in agreement with previous studies (Lemme et al., 2006 for amino acid; Noy and Sklan, 2004 for energy and amino acid; Potter et al., 1966 and Waibel et al., 1995 for Met and Lys; and Bowyer and Waldroup, 1986 for protein). Figure 1 shows a very strong effect of age on CE. This result is similar to previous studies aiming to describe the growth pattern of animals with age using growth functions (Darmani-Kuhi et al., 2003; Schulin-

neuron in the next following hidden layer.

Zeuthen et al., 2008). The calculated values of CE model error measurement showed that the testing set for toms yielded lower values of MS error, mean absolute deviation, mean absolute percentage error, mean relative error, and higher values of R2 compared with the training set.

Fig. 5. Evolved structure of the generalized group method of data handling-type NNs for feed efficiency in tom turkeys. The letters *a, b, c, d*, and *e* stand for input variables of age (wk), ME (kcal/g), CP (% of diet), Met (% of diet), and Lys (% of diet), respectively. This figure illustrates the generated relationship between input variables to reach output.

Conducting a sensitivity analysis (SA) on the obtained polynomial equations reveals the sensitivity of model output to input variables. Hence it is necessary to do sensitivity analysis for any proposed model. In other words, SA increases confidence in the model and its predictions by providing an understanding of how the model responds to changes in its inputs. Moreover, the SA identifies critical regions in the space of the inputs, establishes priorities for research and simplifies the model (Castillo et al.,2008; Saltelli et al., 2008).

For such reason and in line with previous work, another study was designed, titled:

*"Sensitivity analysis of an early egg production predictive model inbroiler breeders based on dietary nutrient intake" (Faridi et al., 2011).* 

Although the use of NN and SA techniques has led to successful application in a broad range of areas (Seyedan&Ching 2006;Lee &Hsiung 2009), the use of SA along with NN models is appeared uncommon in poultry science. The aim of the present study was to use the GMDH-type NN to model early egg production (EEP) in broiler breeders (BB) based on the dietary intake levels of ME, CP, and the two first limiting amino acids, methionine (Met)

Genetic Algorithms and Group Method

65(Suppl. 1):16. Abstr.

91(Suppl.8), 1615–1635.

University. January, 1995

215

287

of Data Handling-Type Neural Networks Applications in Poultry Science 231

Castillo, E., R, Minguez, and C, Castillo.(2008). Sensitivity analysis in optimization and

Darmani- Kuhi, H. Kebreab, H., E. Lopez, S. and France J. (2003). An evaluation of different

Dayhoff, J. E., & de Leo, J.M.(2001). Artificial NNs: opening the black box. Cancer

Dellaert, F.(1995).Toward a biologically defensible model of development.Submitted in

Donaldson,W.E.,Coombs, G.F. & Romoser, G. L.(1956).Studies on energy levels in poultry

utilization and body composition of chicks. Poultry Science 35, 1100–1105. Edwards, H.M.,Denman,F., Abov-Ashour,A. & Nugara, D.(1973).Carcass composition

Faridi, A. ,. Mottaghitalab, M. , Darmani-kuhi, H. France, J.& Ahmadi. H. (2011). Predicting

Farlow,S. J. (1981).The GMDH algorithm of Ivakhnenko. The American Statistician.35, 4:210-

Farrow S.J . (1984). The GMDH algorithm. In self organizing Methods in Modeling: GMDH Type Algorithms ( Ed .S. J. Farrow), pp. 1–24. New York: Marcel Dekker. Fernández,F.H. & F. H. Lozano.(2010). GMDH Algorithms implemented in the intelligent

Fraps, G. S. (1943).Relation of the protein,fat, and energy of the ration to the composition of

GEvoM.(2009).GMDH-type NN Design by Evolutionary Method for Modelling. Rasht, Iran:

Goldberg, D. E.(1989).Genetic Algorithms in Search, Optimization and Machine Learning.

Happel, B. L. M. & J. M.J. Murre. 1994. The Design and Evolution of Modular NN

Gous, R. M. 1998. Making progress in the nutrition of broilers. Poult.Sci. 77:111–117.

carcass fatty acid composition. Poultry Science 52, 934–948.

intake.Journal of Agricultural Science. 150, 87-93

chickens. Poultry Science 22, 421–424.

University of Guilan. Available online at:

Addison-Wesley,Reading, Massachusetts.

Architectures.NNs. 7: 985-1004.

http://research.guilan.ac.ir/gevom/(verified 20 October 2010).

reliability problems. Reliability Eng. and Syst. Safety 93:1788–1800.

and egg strains of chicken. Poult. Sci. 82:1536–1543.

Report No. 813. October 2001. *Department of Automatic Control & Systems Engineering The University of Sheffield Mappin St, Sheffield, S1 3JD, United Kingdom.*  Bowyer, B. L., and P. W. Waldroup. (1986). Evaluation of minimal protein levels for growing

turkeys and development of diets for estimating lysine requirements. Poult. Sci.

growth functions for describing the profile of live weight with time (age) in meat

partial fulfillment of the requirements for the degree of Master of Science. Department of Computer Engineering and Science,Case Western Reserve

rations.1. The effect of calorie-protein ratio of the ration on growth, nutrient

studies. 1. Influences of age, sex, and type of dietary fat supplementation on total

carcass energy content and composition in broilers using the group method of data handling-type NNs. Journal of Agricultural Science, Vol. 149, No2 , pp. 249-254 Faridi A, M.Mottaghitalab & H. Ahmadi .(2012).Sensitivity analysis of an early egg

production predictive model in broiler breeders based on dietary nutrient

identification of a bioprocess.ABCM Symposium Series in Mechatronics-Vol.4:.278-

and lysine (Lys). The SA method was utilized to evaluate the relative importance of input variables on model output and to determine the optimum levels of nutrient intake for obtaining the maximum EEP in BB.

In this study, the GMDH-type NN with GA method was used to develop the EEP in BB. By means of the GMDH algorithm, a model can be represented as a set of quadratic polynomials. In this way, GA are deployed to assign the number of neurons (polynomial equations) in the network and to find the optimal set of appropriate coefficients of the quadratic expressions.

The variables of interest in this model were the dietary intake levels of ME (MJ/bird/day), CP (g/bird/day), Met (g/bird/day), Lys (g/bird/day) and weekly egg production (eggs/bird) during early production (from 24 to 29 weeks of age). Datasets were imported into the software GEvoM for GMDH-type NN training (GEvoM 2009).

Results of the developed GMHD-type NN models revealed close agreement between observed and predicted values of EEP. Results showed that the evolved GMDH-type NNs have been successful in obtaining a model for the prediction of EEP in BB. All input variables were accepted by the model, i.e. the GMDH-type NN provides an automated selection of essential input variables and builds polynomial equations to model EEP.

The advantage of using GMDH-type NN is that which polynomial equations obtained can be used to analyze the sensitivity of output with respect to input variables. SA discusses how and how much changes in the input variables modify the optimal objective function value and the point where the optimum is attained. The simple approach to SA is easy to do, easy to understand, easy to communicate, and applicable with any model.

### **6. Conclusion**

The conclusion was that, genetic algorithm in general and GMDH-type NN in particular may be used as a powerful tool to enhance our ability to predict economic traits, make precise prediction of nutrition requirements, and achieve optimal performance in poultry nutrition and production.

### **7. References**


and lysine (Lys). The SA method was utilized to evaluate the relative importance of input variables on model output and to determine the optimum levels of nutrient intake for

In this study, the GMDH-type NN with GA method was used to develop the EEP in BB. By means of the GMDH algorithm, a model can be represented as a set of quadratic polynomials. In this way, GA are deployed to assign the number of neurons (polynomial equations) in the network and to find the optimal set of appropriate coefficients of the

The variables of interest in this model were the dietary intake levels of ME (MJ/bird/day), CP (g/bird/day), Met (g/bird/day), Lys (g/bird/day) and weekly egg production (eggs/bird) during early production (from 24 to 29 weeks of age). Datasets were imported

Results of the developed GMHD-type NN models revealed close agreement between observed and predicted values of EEP. Results showed that the evolved GMDH-type NNs have been successful in obtaining a model for the prediction of EEP in BB. All input variables were accepted by the model, i.e. the GMDH-type NN provides an automated

The advantage of using GMDH-type NN is that which polynomial equations obtained can be used to analyze the sensitivity of output with respect to input variables. SA discusses how and how much changes in the input variables modify the optimal objective function value and the point where the optimum is attained. The simple approach to SA is easy to do,

The conclusion was that, genetic algorithm in general and GMDH-type NN in particular may be used as a powerful tool to enhance our ability to predict economic traits, make precise prediction of nutrition requirements, and achieve optimal performance in poultry

Ahmadi, H.,M.Mottaghitalab,& N.Nariman-Zadeh.(2007).Group Method of data handling-

AHMADI, H., M. Mottaghitalab,N. Nariman-Zadeh & A. Golian. (2008). Predicting

Amanifard,N., N. Nariman-Zadeh, M. Borji, A. Khalkhali, & A. Habibdoust. (2008).

Energy, Methionine, and Lysine. J. Appl. Poult. Res. 16:494–501

Type NN Prediction of Broiler Performance Based on Dietary Metabolizable

performance of broiler chickens from dietary nutrients using group method of data handling-type neural networks. British Poultry Science Volume 49, Number 3, pp.

Modelling and pareto optimization of heat transfer and flow coefficients in microchannels. Energy Conversion & Management, issue 2, Vol.49, pp 311-325. Anastasakis, L.& N. Mort.(2001).The development of self-organization techniques in

modeling: A review of the group method of data handling (GMDH). Research

selection of essential input variables and builds polynomial equations to model EEP.

into the software GEvoM for GMDH-type NN training (GEvoM 2009).

easy to understand, easy to communicate, and applicable with any model.

obtaining the maximum EEP in BB.

quadratic expressions.

**6. Conclusion** 

**7. References** 

nutrition and production.

315—320

Report No. 813. October 2001. *Department of Automatic Control & Systems Engineering The University of Sheffield Mappin St, Sheffield, S1 3JD, United Kingdom.* 


http://research.guilan.ac.ir/gevom/(verified 20 October 2010).


Genetic Algorithms and Group Method

York,NY.

186.

301.

of Data Handling-Type Neural Networks Applications in Poultry Science 233

Nariman-Zadeh, N.,Darvizeh ,A.& Ahmad-Zadeh,G.R.(2003).Hybrid genetic design of

Mechanical Engineers B: Journal of Engineering Manufacture 217, 779–790. Noy, Y., and D. Sklan.(2004). Effects of metabolizable energy and amino acid levels on turkey performance from hatch to marketing. J. Appl. Poult. Res. 13:241–252. NRC. (1994).Nutrient Requirements of Poultry.9th rev. ed. Natl. Acad. Sci.,Washington,DC. Oberstone,J.(1990).Management Science-Concepts,Insights,and cations. West Publ. Co.,New

Pesti, G.M. & Bakalli, R. I.(1997).Estimation of the composition of broiler carcasses from their

Phillips, R. D.(1981). Linear and nonlinear models for measuring protein nutritional

Potter, L. M., A. T. Leighton Jr., and C. E. Howes. (1966). Methionine and lysine

Rosa, P.S.,Faria Filho,D. E., Dahlke, F., Vieira, B. S., Macari, M. & FURLAN, R. L.(2007).

Rosenblatt, F.(1958). The perceptron: A probabilistic model for information storage and

Roush,W.B.& Wideman, R.F.(2000) Evaluation of broiler growth velocity and acceleration in relation to pulmonary hypertension syndrome. Poultry Science, 79: 180–191. Roush, W.B., Dozier III, W.A. & Branton, S.L. (2006) Comparison of Gompertz and NNs

Roush, W.B., Kirby, Y.K., Cravener, T.L. & Wideman JR, R.F.(1996) Artificial NN prediction

Salle, C.T., Guahyba, A.S., Wald, V.B., Silva, A.B., Salle, F.O. & Aacimento, V.P. (2003) Use of

Saltelli, A., M. Ratto, T. Andres, F. Campolongo,J. Cariboni, D. Gatelli, M. Saisana and S. Tarantola.2008.Global Sensitivity Analysis.The Primer. John Wiley&Sons,Ltd. Schulin-Zeuthen, M. Kebreab, E. Dijkstra, J. Lopez, S. Bannink, A. Darmani Kuhi, H. and

Seyedan,b. & Ching, C. Y. (2006). Sensitivity analysis of freestream turbulence parameters

Sibbald, I. R.(1987). Estimation of bioavailable amino acids in feeding stuffs for poultry and

Toledo, G. S. P., Lopez, J. & Costa,P. T. C.(2004). Yield and carcass composition of broilers

describing growth in pigs. Anim. Feed Sci. Technol. 143:314–327.

artificial NNs to estimate production variables of broilers breeders in the

France, J. (2008). A comparison of the Schumacher with other functions for

on stagnation region heat transfer using a netural network. International journal of

pigs: A review with emphasis on balance experiments. Can. J. Anim. Sci. 67:221–

fed with diets based on the concept of crude protein or ideal protein.Revista

organization in the brain.Psychological Review. 68: 386-408

models of broiler growth. Poultry Science, 85: 794–797.

production phase. British Poultry Science, 44: 211–217.

Heat and Fluid Flow 27,1061-1068.

Brasileira de Ciência Avícola 6,219–224.

of ascites in broilers. Poultry Science, 75: 1479–1487.

supplementation of turkey starting diets containing varying levels of protein. Poult.

Performance and carcass characteristics of broiler chickens with different growth potential and submitted to heat stress. Revista Brasileira de Ciência Avícola 9, 181 –

specific gravity. Poultry Science 76, 948–951.

quality.J.Nutr.111:1058–1066.

Sci. 45:1117–1118.

GMDH-Type NNs using singular value decomposition for modeling and prediction of the explosive cutting process. Proceedings of the Institution of


Havenstein, G. B., P. R. Ferket, J. L. Grimes, M. A. Qureshi, and K. E. Nestor. 2007.

Holland, J. H.(1975).*Adaptation in Natural and Artificial Systems*. Ann Arbor, MI: University of

Hruby, M., M. L. Hamre,& C. N. Coon.(1996).Non-linear and linear functions in body

Iba, H., deGaris, H. and Sato,T. (1996). A numerical Approach to Genetic Programming for

Ivakhnenko A.G.(1970).Heuristic self-organization in problems of engineering

Ivakhnenko A.G.(1971).Polynomial theory of complex systems. IEEE Transactions on

Ivakhnenko, A.G. (1966).Group Method of Data Handling-A rival of the method of

Kubena, L .F.,Lott,B.D., Deaton, J.W., Reece, F.N. & May, J.D. (1972). Body composition of

Lee,C.J. & Hsiung, T. K. (2009).Sensitivity analysis on a multilayer perceptron model for recognizing liquefaction cases. Computers and Geotechnics 36, 1157-1163. Lemme, A.,U. Frackenpohl, A. Petri, & H. Meyer.(2006).Response of male BUT Big 6 turkeys

Lemke, F.,and J.A.Mueller. 2003. Medical data analysis using self-organizing data mining

Madala,H.R.& A.G.Ivakhnenko.(1994).Inductive learning algorithms for complex Systems Modeling. CRC Press:, ISBN:0-8493-4438-7.Boca Raton, Ann Arbor,London Tokyo. Makkar,H.P.S.(2008).A review of the use of isotopic and nuclear techniques in animal

Miroslav Šnorek,P.K.(2006).Inductive modeling World Wide the State of the Art.Report of investigation. Dept.of Computer Science and Engineering, Karlovo nam. Mitchell, M. & S.Forrest (1994).Genetic Algorithms and Artificial Life. Artificial Life, Vol. 1,

Mueller,J.A.,&F.Lemke.(2000).Self-Organizing Data Mining:An Intelligent Approach to

Mottaghitalab, M., Faredi, A., Darmani-kuhi, H., France , J. & Ahmadi, H. (2010). Predicting

Mutasem H. N. M.(2004).An evolutionary method for term selection in the Group Method of

Nariman-Zadeh,N.,A.Darvizeh,A.Jamali, & A.Moieni.(2005). Evolutionary design of generalized

caloric and feed efficiency in turkeys using the group method of data handling-

Data Handling. Automatic Control & Systems Engineering, University of Sheffield,

polynomial NNs for modeling and prediction of explosive forming process. J.

Extract Knowledge from Data.Libri Publ. Co., Hamburg, Germany.

chicks as influenced by environmental temperature and selected dietary factors.

System Identification. Evolutionary Computation 3(4):417- 452,.

stochastic approximation.*Soviet Automatic Control*.13:43-71. Koza,J.R.(1980).*Genetic Programming*. Cambridge, MA: The MIT Press.

to varying amino acid feeding programs. Poult.Sci. 85:652–660.

production. Animal Feed Science and Technology 140, 418–443.

technologies. Syst. Anal. Model.Simul.10: 1399–1408.

type neural networks. Poultry Science 89, 1325–1331.

Mater.Process. Technol. 164–165:1561–1571.

conversion. Poult. Sci. 86:232–240.

cybernetics.Automatica.l.6:.207-219,.

Poultry Science 5, 517–522.

No. 3, 267-289

11-14.

protein growth.J. Appl. Poult. Res. 5:109–115.

Systems, Man,and Cybernetics.SMC-1, 364–378

Michigan Press.

Comparison of the performance of 1966- versus 2003-type turkeys when fed representative 1966 and 2003 turkey diets: Growth rate, livability, and feed


**12** 

*1USA 2Russia 3Canada* 

**New Approaches to Designing Genes by** 

*1The I. M. Sechenov Institute of Evolutionary Physiology and Biochemistry,* 

The field of Evolutionary Computation (EC) has been inspired by ideas from the classical theory of biological evolution, with, in particular, the components of a population from which reproductive parents are chosen, a reproductive protocol, a method for altering the genetic information of offspring, and a means for testing the fitness of offspring in order to include them in the population. In turn, impressive progress in EC – understanding the reasons for efficiencies in evolutionary searches - has begun to influence scientific work in the field of molecular evolution and in the modeling of biological evolution (Stemmer, 1994a,b; van Nimwegen et al. 1997; 1999; Crutchfield & van Nimwegen, 2001). In this chapter, we will discuss how developments in EC, particularly in the area of crossover operators for Genetic Algorithms (GA), provide new understanding of evolutionary search efficiencies, and the impacts this can have for biological molecular evolution, including

GA approaches have five particular elements: encoding (the 'chromosome'); a population; a method for selecting parents and making a child chromosome from the parents' chromosomes; a method for altering the child's chromosomes (mutation and crossover/recombination); criteria for fitness; and rules, based on fitness, by which offspring are included into the population (and parents retained). We will discuss our work and others' on each of these aspects, but our focus is on the substantial efficiencies that can be found in the alteration of the child chromosome step. For this, we take inspiration from real

Traditional GA, using random point mutations, indicates that such a mechanism would be too slow to account for the observed speed of biological evolution (e.g. Shapiro, 2010). This suggests that other more complicated mutational mechanisms are acting (Shapiro, 1999,

**1. Introduction** 

directed evolution in the test tube.

biological reproduction mechanisms.

**1.1 Biological evolution by random point mutations?** 

**Evolution in the Computer** 

*3British Columbia Institute of Technology,* 

*2Stony Brook University,* 

Alexander V. Spirov and David M. Holloway

