# Multiattribute Response of Maize Genotypes Tested in Different Coastal Regions of Brazil

Authors: Lúcio Borges de Araújo; Mario Varela Nualles; Mirian Fernandes Carvalho Araújo; Carlos Tadeu dos Santos Dias
Journal: International Journal of Agronomy
 (2011)

---

## Abstract

This work applies the three mode principal components analysis to analyze simultaneously the multiple attributes; to fit of models with additive main effects and multiplicative interaction effects (AMMI models) and the regressions models on sites (SREG models); to evaluate, respectively, the multivariate response of the genotype × environment interaction and the mean response of 36 genotypes of corn tested in 4 locations in Brazil. The results were presented by joint plots to identify the best genotypes for their adaptability and performance in the set of attributes.


---

## Body

## 1. Introduction

The multienvironment trials are conducted to estimate the genetic stability, to evaluate the performance of the genotypes in different environmental conditions, and to quantify and interpret the genotype × environment interaction (GEI), in such a way to select the best genotypes that will be recombined and planted in other years and other environments in the next selection cycle. Some statistical models are used to evaluate the behavior of the genotypes; the most recently proposed methods are based on the singular value decomposition of the GEI matrix.In [1], the two-way interaction using the principal components analysis (PCA) are explained and considered in their model only one multiplicative term. An extension of these models was made in [2], which applies the PCA to decompose the two-way interaction in several multiplicative terms. This model was called as additive main effects and multiplicative interaction effects model (AMMI) by Gauch [3].Another class of the Linear-bilinear model used in multienvironment trials is the sites regression model (SREG) described by [4]. In this case, the main objective is to evaluate the response of the genotype in each environment. The basic difference in relation to AMMI models is that the effect of genotype is introduced into the residual interaction. This model was used in [5] for clustering of environmental data with no cross-interaction.Gabriel [6] describes the least-square adjustment of the AMMI model, first estimating the additive effects of the model and then making the decomposition in singular values of the matrix (Z) of interaction residuals. The results can be interpreted by graphics called biplot [7] that reflects in reduced dimensions the most important aspects of the GEI, in the case the AMMI model, or the performance of genotypes in different environments, for the case of the model SREG.The researchers in the multienvironmental trials evaluate multiple attributes and must select the best genotypes taking in to account its adaptability and performance in the set of attributes. With the AMMI and SREG models, it is possible to study the genotypes in a single attribute. However, if we want a response with multiple attributes, it is necessary to apply a statistical technique that allows working with three-dimensional structure data.In [8], a three mode principal components analysis (PCA3) is proposed, in order to find the least-square estimates in the model of Tucker [9]. This procedure was used to study the GEI with multiples attributes [10–13]. The method also was applied by [14, 15] to study the three-way interaction in agricultural experiments where genotypes are tested in different locations over several years.In this work, we used the three mode principal component analysis to explain, in reduced dimensions, the most important aspects of multiattribute performance of the 36 genotypes of maize tested in 4 regions of Brazil. First, we studied the GEI using the AMMI model with multiple attributes as described in [15, 16] and finally studied the genotypes in relation to the mean environment for multiple attributes as described in [17].


## 2. Materials and Methods

Data are from experiments with 36 genotypes of maize (first mode,i=1,…,36), evaluated in 4 locations of Brazil (S1: Rio Branco-AC, S2: Milagres-CE, S3: Linhares-ES, cohesive soil Oxisol, S4: Linhares-ES, soil-LVA-Distrophic LVD11; second mode, j=1,…,4). In each of the experiments, 10 attributes were evaluated: X1: flowering (days), X2: plant height (cm), X3: corn-cob length (cm), X4: percentage of plants laid, X5: percentage of plants broken, X6: final stand, X7: number of spikes, X8: percentage of patients corn-cobs, X9: weight of spikes (kg/ha), and X10: weight of grains (kg/ha) (third mode, k=1,…,10).The data were obtained using a simple lattice design in duplicate 6 × 6 with 4 replications. Each plot consisted of two rows of 5 m, used entirely. The spacing between rows was 1 m and between hills it was 0.40 m. Each plot consisted of 26 hills and each hole with 3 seeds. The seeds used in the tests were sent by the National Center for Research in Maize and Sorghum—CNPMS of Empresa Brasileira de Pesquisa Agropecuária—EMBRAPA, already counted, and contained in envelopes coded with the number of genotype and number of the parcel.


### 2.1. Statistical Processing

Each attributek(k=1,…,10) was calculated according to model the estimates corresponding to interactions: (yijk.-yi.k.-y.jk.+y..k.) for AMMI model or (yijk.-y.jk.) for SREG model, where yijk. represents the mean corresponding to attribute k for the genotype i in location j; yi.k. represents the mean corresponding to the attribute k for the genotype i; y.jk. represents the mean corresponding to the attribute k in location j; y..k. represents the overall mean of the data for the attribute k. Two tridimensional arrays (Z) were constructed with the values obtained using the AMMI model or the SREG model, that is, for each model, the value zijk in the array Z represents the residue to the genotype j in the location j, for the attribute k. Each array was applied to three mode principal components analysis. Therefore, this methodology can be used for any dataset.


### 2.2. Three Mode Principal Component Analysis

The method consists of finding the least-square estimators for the Tucker model:(1)zijk.=∑p=1P∑q=1Q∑r=1Raipbjqckrgpqr,
where aip, bjq, and ckr are the elements of the principal components matrix (A, B, and C) associated with each respective mode, and in this example, A is the principal components matrix for the genotype mode, B is the principal component matrix for the location mode, and C is the principal component matrix for the attribute. G is a core array of the three way, where the element gpqr represents the relationship between the pth component of the first mode with the qth component of the second mode and with the rth component of the third mode. P is the rank of the matrix Z1;2⊂3 that contains the residuals of the interaction with all the variables, calculated as an AMMI model or as an SREG model. In this matrix, the rows represent the levels of the genotypes and the columns represent combinations of the levels of locations and attributes. Similarly, Q and R are defined as the rank of matrices Z2;1⊂3 and Z3;1⊂2.To select the number of principal components retained for each mode, the algorithm proposed by [18] was used, based on successive eliminations until finding the optimal solution.


### 2.3. Biplot Representation for Three Marker Matrices

A biplot is a simultaneous representation of rows and columns of a matrix (two-way array) in terms of directions and projections [7]. When we have several matrices, or when we are working with data contained in a three-way array, it is necessary to project on the principal components of one of the modes. Suppose we are going to project onto the components of the second mode (locations) and Q1 represents the number of principal components retained in this mode; then we have to obtain the matrices A×Gq×B′,q=1,…,Q1, and require one biplot for each matrix. Gq is the part of G associated with the component q of the second mode. This procedure is described in [8] and commonly called a joint plot.


## 2.1. Statistical Processing

Each attributek(k=1,…,10) was calculated according to model the estimates corresponding to interactions: (yijk.-yi.k.-y.jk.+y..k.) for AMMI model or (yijk.-y.jk.) for SREG model, where yijk. represents the mean corresponding to attribute k for the genotype i in location j; yi.k. represents the mean corresponding to the attribute k for the genotype i; y.jk. represents the mean corresponding to the attribute k in location j; y..k. represents the overall mean of the data for the attribute k. Two tridimensional arrays (Z) were constructed with the values obtained using the AMMI model or the SREG model, that is, for each model, the value zijk in the array Z represents the residue to the genotype j in the location j, for the attribute k. Each array was applied to three mode principal components analysis. Therefore, this methodology can be used for any dataset.


## 2.2. Three Mode Principal Component Analysis

The method consists of finding the least-square estimators for the Tucker model:(1)zijk.=∑p=1P∑q=1Q∑r=1Raipbjqckrgpqr,
where aip, bjq, and ckr are the elements of the principal components matrix (A, B, and C) associated with each respective mode, and in this example, A is the principal components matrix for the genotype mode, B is the principal component matrix for the location mode, and C is the principal component matrix for the attribute. G is a core array of the three way, where the element gpqr represents the relationship between the pth component of the first mode with the qth component of the second mode and with the rth component of the third mode. P is the rank of the matrix Z1;2⊂3 that contains the residuals of the interaction with all the variables, calculated as an AMMI model or as an SREG model. In this matrix, the rows represent the levels of the genotypes and the columns represent combinations of the levels of locations and attributes. Similarly, Q and R are defined as the rank of matrices Z2;1⊂3 and Z3;1⊂2.To select the number of principal components retained for each mode, the algorithm proposed by [18] was used, based on successive eliminations until finding the optimal solution.


## 2.3. Biplot Representation for Three Marker Matrices

A biplot is a simultaneous representation of rows and columns of a matrix (two-way array) in terms of directions and projections [7]. When we have several matrices, or when we are working with data contained in a three-way array, it is necessary to project on the principal components of one of the modes. Suppose we are going to project onto the components of the second mode (locations) and Q1 represents the number of principal components retained in this mode; then we have to obtain the matrices A×Gq×B′,q=1,…,Q1, and require one biplot for each matrix. Gq is the part of G associated with the component q of the second mode. This procedure is described in [8] and commonly called a joint plot.


## 3. Results and Discusions

### 3.1. AMMI Multiattribute

Table1 shows the results of the 6 × 2 × 3 solution obtained by applying the three mode principal component analysis to the three-way array with the residuals of interaction obtained for the AMMI model with 10 attributes simultaneously. This solution explained 44.92% of the total variability, therefore, this percentage is significant in relation to ten attributes measured in the research. The components associated with the genotype mode explain, respectively, 12.85%, 11.70%, 6.37%, 6.22%, 4.93%, and 2.84%, the components of the location mode explain, respectively, 25.49% and 19, 44%, and the components of the attribute mode explain, respectively, 22.6%, 11.6%, and 11.0% of the total variability.Table 1
6 × 2 × 3 solution for AMMI model with multiples attributes.
(a)
MatrixA
p1p2p3p4p5p6G10.16−0.1−0.20.30.010.29G20.030.120.080.02−0.10.16G30.12−0−0.2−0−0−0.2G40.1−0.20.330.05−0.3−0.2G50.260.19−0.20.12−00.12G6−0.1−0.3−0.10.06−0.1−0G70−0.1−0.30.230.030.08G80.080.020.20.17−00.05G90.16−0.30.080.1−0.20.21G10−00.020.08−00.060.09G110.15−00.19−0.10.220.22G120.25−0.30.04−0.10.01−0.1G130.08−0.10.06−0.1−0.2−0G140.190.280.16−0.20.120.07G15−0.2−0.30.340.130.58−0.1G16−0.300.11−0.1−0.10.09G17−0.20.220.11−00.020.09G18−0.1−0.2−0.1−0.5−0.2−0G190.020.160.330.25−0.3−0.4G20−0.20.080.08−0.1−0.10.02G21−0.10.14−00.0200.15G220.08−0.10.03−0.30.03−0.1G230.110.060−0.3−0.10.19G24−0.30.110.020.15−0.2−0G250.01−0.2−0.20.170.030.21G260.180.130.02−00.28−0.1G270.10.06−0.10.120.020.07G28−00.21−0.100.15−0.3G290.030.03−00.070.05−0.3G300.090.1−0.1−00.010.09G310.02−0.1−0.3−0.30.05−0.4G32−0.5−0−0.20.10.09−0.1G33−00.180.060.040.04−0.1G3400.130−0.10.22−0.2G35−0.2−0.1−0.10−0.1−0.1G36−00.15−0.10.08−0.1−0.2(b)
MatrixB
q1q2S10.42−0.8S20.560.64S3−0.40.11S4−0.60(c)
MatrixC
r1r2r3X10.12−0.10.23X2−0.10.060.19X3−0.1−00.12X40.25−0.6−0.6X50.010.24−0.5X6−0.5−0−0X7−0.3−0.70.39X80.240.070.27X9−0.50.04−0.2X10−0.50.01−0.2(d)
ArrayG̲
q1r1q1r2q1r3q2r1q2r2q2r3P1−5.4−0.9−4.1−5.3−0.71.48P26.06−2−0.2−5.4−0.30.63P3−1−5.50.920.87−1.7−1.7P41.43−0.9−5.11.911.76−1.5P5−0.8−1.11.5−0.94.93−0.5P60.2−1.2−0.21.120.473.74Component 1 of the local mode is characterized by all the sites with relatively high weights (MatrixB, Table 1). This component is described by contrast between sites S1 and S2 (Rio Branco-AC and Milages-CE) versus sites S3 and S4 (Linhares-ES oxisol cohesive soil and Linhares-ES LVA-dystrophic-LVD11 soil).Figure1 shows the most important genotypes × locations interactions for the different attributes evaluated, mainly for the attribute X4 (percentage of plants laid). Genotype 19 interacts positively in locations S1 and S2 and negatively in locations S3 and S4 for attribute X4. Genotypes 18 and 31 have a behavior completely contrary to genotype 19, that is, they interact positively in the locations S3 and S4 and interact negatively in locations S1 and S2 for the attribute X4. In regard to variable X2 (plant height), the behavior of these genotypes is totally contrary to variable X4, for example, the genotype 19 interacts negatively in locations S1 and S2 and positively in the locations S3 and S4 for the attribute X2.Figure 1
Joint plot conditioned on the first component the location mode for AMMI model.Figure1 also shows the contrast between the genotypes 9 and 12 of the genotypes 17, 24, and 32. The first group interacts positively in the attributes X6, X9, and X10 for locations S1 and S2 and negatively in these attributes for the locations S3 and S4. Genotypes 17, 24, and 32 have positive interaction in the attributes X1 and X8 in the locations S1 and S2 but have negative interactions in locations S3 and S4 in these attributes.Figure2 shows the projection biplot onto the second component of locations, which is dominated by contrast of the sites S1 and S2 with respective weights (−0.75 and 0.64) (Matrix B, Table 1). In this graphic, it is possible to study the genotypes × locations interaction in different attributes for the locations S1 and S2. For example, genotype 15 interacts positively in the local S1 for attribute X7 and interacts negatively in local S2 for this attribute. There is also the contrast of genotypes 4 and 18 in relation to genotype 26. The genotype 26 interacts positively with the attribute X4 in local S1 and interacts negatively in local S2 for this attribute, then in S2 there are a few percentage of plant laid, and in S1 many percentage of plant laid. The contrast between the genotypes 5 and 14 in relation to the genotypes 6 and 35 for attributes X6, X9, and X10 can be observed. The first group had positive interactions with S2 and negative with S1 and showed that in the S2 there are high values of final stand, weight of spikes, and weights of grains.Figure 2
Joint plot conditioned on the second component the location mode for AMMI model.


### 3.2. SREG Multiattribute

Table2 shows the results of the 5 × 2 × 3 solution obtained by applying the three mode principal components analysis to the three-way array containing the residuals of interaction for the SREG model with 10 attributes simultaneously. This solution explained 53.70% of the total variability. The components associated with the genotype mode explain, respectively, 27.69%, 13.34%, 6.96%, 3.77%, and 1.94%, the components of the local mode explain 46.23% and 7.48%, and the components of the mode attribute explain 27.96%, 17.31%, and 8.42%. In this analysis the first component of the local mode is characterized by all the sites with relatively high weights (Matrix B, Table 2).Table 2
5 × 2 × 3 solution for SREG model with multiples attributes.
(a)
MatrixA
p1p2p3p4p5G10.07−00.28−0.40G20−0.2−0.10.12−0.1G30.140.11−0.1−0.10.07G4−0.1−0.50.02−0.10.03G50.01−0.1−0−0.10.1G60.02−0.10.1−00.01G7−0.1−0.20.15−0.10.28G80.090.25−0.2−0.10.11G90.060.13−0.2−0.30.06G100.3−00.020.2−0G110.230.080−0.1−0.2G120.07−0−0.2−0.30.05G13−0.30.04−0.1−0.20.03G14−0.10.18−0.3−0.1−0.1G15−0.30.05−0.30.03−0.1G16−0.30.25−00.07−0.1G17−0.20.160.170.13−0.1G180.110.030.18−0−0.3G19−00.160.14−0.20.22G20−0.1−0.20.070.27−0.2G21−0.2−00.070.030.15G22−0.1−0.2−0.1−0−0.5G23−0−0.20.040.04−0.2G24−0.1−00.080.370.32G250.01−0.2−0.100.32G260.13−0.1−0.20.050.02G270.02−0.1−0.20.10.23G280.240.14−0.10.260.04G290.020.03−0.10.010.03G300.08−0.1000G310.330.070.21−0.1−0.1G32−0.20.290.210.230.03G33−0.10.030.18−0.1−0.1G340.280.09−0.10.21−0G350.110.080.150.160.05G36−0.10.150.34−0.20(b)
MatrixB
q1S10.54S20.45S30.46S40.53(c)
MatrixC
r1r2r3X10.45−0.20.21X20.45−0.20.05X30.5−0.3−0.1X40.21−0−0.3X5−0.030.08−0.9X6−0.21−0.30.06X7−0.17−0.2−0X8−0.170.370.31X9−0.29−0.50.01X10−0.28−0.50(d)
ArrayG̲
q1r1q1r2q1r3q2r1q2r2q2r3P119.5−2.9−0.41.312.42−0.9P22.2513.2−1.5−2−2.2−1.1P30.15−2−8.6−1.5−4.4−0.6P40.62−1.23.53−2.6−5.2−2.4P50.470.140.630.34−24.83Figure3 studies the behavior of the genotypes in relation to the mean environmental for different attributes. Genotypes 10, 11, 28, 31, and 34 are above mean environmental for the attributes X1, X2, and X3 in the 4 locations and, in contrast, the genotypes 13, 14, 15, 16, 17, and 32 have a performance below the mean environmental in these attributes. Genotype 4 has good performance in the variables X9 and X10 to the 4 sites in contrast to genotype 8 which is below the mean environmental in these attributes in the same locations. An opposite behavior occurs in the variable X8 in which genotype 4 is below mean environmental and genotypes 8, 14, and 32 have good performance.Figure 3
Joint plot conditioned on the first component the location mode for SREG model.The second component of the local mode is characterized by localS2 and S4 with respective weights of −0.73 and 0.61 (Table 2, Matrix B). This component reflects the contrast between these two locations. In the corresponding joint plot (Figure 4), it can be observed that genotypes 1, 4, 7, 9, 12, and 25 have good behavior with the variable X8 for the local S4, but are below the mean environmental in local S2 for this same attribute. For the same genotypes, the opposite happens, they have good performance in variable X5 for the local S2, but they are below mean environmental for the same attribute to the local S4.Figure 4
Joint plot conditioned on the second component the location mode for SREG model.The fundamental difference between AMMI models with multiple attributes and SREG models with multiple attributes is that in the first case study, the two-way interaction considering all the attributes, while in the SREG model, it is possible to study the multivariate behavior of genotypes in relation to the mean environmental.


## 3.1. AMMI Multiattribute

Table1 shows the results of the 6 × 2 × 3 solution obtained by applying the three mode principal component analysis to the three-way array with the residuals of interaction obtained for the AMMI model with 10 attributes simultaneously. This solution explained 44.92% of the total variability, therefore, this percentage is significant in relation to ten attributes measured in the research. The components associated with the genotype mode explain, respectively, 12.85%, 11.70%, 6.37%, 6.22%, 4.93%, and 2.84%, the components of the location mode explain, respectively, 25.49% and 19, 44%, and the components of the attribute mode explain, respectively, 22.6%, 11.6%, and 11.0% of the total variability.Table 1
6 × 2 × 3 solution for AMMI model with multiples attributes.
(a)
MatrixA
p1p2p3p4p5p6G10.16−0.1−0.20.30.010.29G20.030.120.080.02−0.10.16G30.12−0−0.2−0−0−0.2G40.1−0.20.330.05−0.3−0.2G50.260.19−0.20.12−00.12G6−0.1−0.3−0.10.06−0.1−0G70−0.1−0.30.230.030.08G80.080.020.20.17−00.05G90.16−0.30.080.1−0.20.21G10−00.020.08−00.060.09G110.15−00.19−0.10.220.22G120.25−0.30.04−0.10.01−0.1G130.08−0.10.06−0.1−0.2−0G140.190.280.16−0.20.120.07G15−0.2−0.30.340.130.58−0.1G16−0.300.11−0.1−0.10.09G17−0.20.220.11−00.020.09G18−0.1−0.2−0.1−0.5−0.2−0G190.020.160.330.25−0.3−0.4G20−0.20.080.08−0.1−0.10.02G21−0.10.14−00.0200.15G220.08−0.10.03−0.30.03−0.1G230.110.060−0.3−0.10.19G24−0.30.110.020.15−0.2−0G250.01−0.2−0.20.170.030.21G260.180.130.02−00.28−0.1G270.10.06−0.10.120.020.07G28−00.21−0.100.15−0.3G290.030.03−00.070.05−0.3G300.090.1−0.1−00.010.09G310.02−0.1−0.3−0.30.05−0.4G32−0.5−0−0.20.10.09−0.1G33−00.180.060.040.04−0.1G3400.130−0.10.22−0.2G35−0.2−0.1−0.10−0.1−0.1G36−00.15−0.10.08−0.1−0.2(b)
MatrixB
q1q2S10.42−0.8S20.560.64S3−0.40.11S4−0.60(c)
MatrixC
r1r2r3X10.12−0.10.23X2−0.10.060.19X3−0.1−00.12X40.25−0.6−0.6X50.010.24−0.5X6−0.5−0−0X7−0.3−0.70.39X80.240.070.27X9−0.50.04−0.2X10−0.50.01−0.2(d)
ArrayG̲
q1r1q1r2q1r3q2r1q2r2q2r3P1−5.4−0.9−4.1−5.3−0.71.48P26.06−2−0.2−5.4−0.30.63P3−1−5.50.920.87−1.7−1.7P41.43−0.9−5.11.911.76−1.5P5−0.8−1.11.5−0.94.93−0.5P60.2−1.2−0.21.120.473.74Component 1 of the local mode is characterized by all the sites with relatively high weights (MatrixB, Table 1). This component is described by contrast between sites S1 and S2 (Rio Branco-AC and Milages-CE) versus sites S3 and S4 (Linhares-ES oxisol cohesive soil and Linhares-ES LVA-dystrophic-LVD11 soil).Figure1 shows the most important genotypes × locations interactions for the different attributes evaluated, mainly for the attribute X4 (percentage of plants laid). Genotype 19 interacts positively in locations S1 and S2 and negatively in locations S3 and S4 for attribute X4. Genotypes 18 and 31 have a behavior completely contrary to genotype 19, that is, they interact positively in the locations S3 and S4 and interact negatively in locations S1 and S2 for the attribute X4. In regard to variable X2 (plant height), the behavior of these genotypes is totally contrary to variable X4, for example, the genotype 19 interacts negatively in locations S1 and S2 and positively in the locations S3 and S4 for the attribute X2.Figure 1
Joint plot conditioned on the first component the location mode for AMMI model.Figure1 also shows the contrast between the genotypes 9 and 12 of the genotypes 17, 24, and 32. The first group interacts positively in the attributes X6, X9, and X10 for locations S1 and S2 and negatively in these attributes for the locations S3 and S4. Genotypes 17, 24, and 32 have positive interaction in the attributes X1 and X8 in the locations S1 and S2 but have negative interactions in locations S3 and S4 in these attributes.Figure2 shows the projection biplot onto the second component of locations, which is dominated by contrast of the sites S1 and S2 with respective weights (−0.75 and 0.64) (Matrix B, Table 1). In this graphic, it is possible to study the genotypes × locations interaction in different attributes for the locations S1 and S2. For example, genotype 15 interacts positively in the local S1 for attribute X7 and interacts negatively in local S2 for this attribute. There is also the contrast of genotypes 4 and 18 in relation to genotype 26. The genotype 26 interacts positively with the attribute X4 in local S1 and interacts negatively in local S2 for this attribute, then in S2 there are a few percentage of plant laid, and in S1 many percentage of plant laid. The contrast between the genotypes 5 and 14 in relation to the genotypes 6 and 35 for attributes X6, X9, and X10 can be observed. The first group had positive interactions with S2 and negative with S1 and showed that in the S2 there are high values of final stand, weight of spikes, and weights of grains.Figure 2
Joint plot conditioned on the second component the location mode for AMMI model.


## 3.2. SREG Multiattribute

Table2 shows the results of the 5 × 2 × 3 solution obtained by applying the three mode principal components analysis to the three-way array containing the residuals of interaction for the SREG model with 10 attributes simultaneously. This solution explained 53.70% of the total variability. The components associated with the genotype mode explain, respectively, 27.69%, 13.34%, 6.96%, 3.77%, and 1.94%, the components of the local mode explain 46.23% and 7.48%, and the components of the mode attribute explain 27.96%, 17.31%, and 8.42%. In this analysis the first component of the local mode is characterized by all the sites with relatively high weights (Matrix B, Table 2).Table 2
5 × 2 × 3 solution for SREG model with multiples attributes.
(a)
MatrixA
p1p2p3p4p5G10.07−00.28−0.40G20−0.2−0.10.12−0.1G30.140.11−0.1−0.10.07G4−0.1−0.50.02−0.10.03G50.01−0.1−0−0.10.1G60.02−0.10.1−00.01G7−0.1−0.20.15−0.10.28G80.090.25−0.2−0.10.11G90.060.13−0.2−0.30.06G100.3−00.020.2−0G110.230.080−0.1−0.2G120.07−0−0.2−0.30.05G13−0.30.04−0.1−0.20.03G14−0.10.18−0.3−0.1−0.1G15−0.30.05−0.30.03−0.1G16−0.30.25−00.07−0.1G17−0.20.160.170.13−0.1G180.110.030.18−0−0.3G19−00.160.14−0.20.22G20−0.1−0.20.070.27−0.2G21−0.2−00.070.030.15G22−0.1−0.2−0.1−0−0.5G23−0−0.20.040.04−0.2G24−0.1−00.080.370.32G250.01−0.2−0.100.32G260.13−0.1−0.20.050.02G270.02−0.1−0.20.10.23G280.240.14−0.10.260.04G290.020.03−0.10.010.03G300.08−0.1000G310.330.070.21−0.1−0.1G32−0.20.290.210.230.03G33−0.10.030.18−0.1−0.1G340.280.09−0.10.21−0G350.110.080.150.160.05G36−0.10.150.34−0.20(b)
MatrixB
q1S10.54S20.45S30.46S40.53(c)
MatrixC
r1r2r3X10.45−0.20.21X20.45−0.20.05X30.5−0.3−0.1X40.21−0−0.3X5−0.030.08−0.9X6−0.21−0.30.06X7−0.17−0.2−0X8−0.170.370.31X9−0.29−0.50.01X10−0.28−0.50(d)
ArrayG̲
q1r1q1r2q1r3q2r1q2r2q2r3P119.5−2.9−0.41.312.42−0.9P22.2513.2−1.5−2−2.2−1.1P30.15−2−8.6−1.5−4.4−0.6P40.62−1.23.53−2.6−5.2−2.4P50.470.140.630.34−24.83Figure3 studies the behavior of the genotypes in relation to the mean environmental for different attributes. Genotypes 10, 11, 28, 31, and 34 are above mean environmental for the attributes X1, X2, and X3 in the 4 locations and, in contrast, the genotypes 13, 14, 15, 16, 17, and 32 have a performance below the mean environmental in these attributes. Genotype 4 has good performance in the variables X9 and X10 to the 4 sites in contrast to genotype 8 which is below the mean environmental in these attributes in the same locations. An opposite behavior occurs in the variable X8 in which genotype 4 is below mean environmental and genotypes 8, 14, and 32 have good performance.Figure 3
Joint plot conditioned on the first component the location mode for SREG model.The second component of the local mode is characterized by localS2 and S4 with respective weights of −0.73 and 0.61 (Table 2, Matrix B). This component reflects the contrast between these two locations. In the corresponding joint plot (Figure 4), it can be observed that genotypes 1, 4, 7, 9, 12, and 25 have good behavior with the variable X8 for the local S4, but are below the mean environmental in local S2 for this same attribute. For the same genotypes, the opposite happens, they have good performance in variable X5 for the local S2, but they are below mean environmental for the same attribute to the local S4.Figure 4
Joint plot conditioned on the second component the location mode for SREG model.The fundamental difference between AMMI models with multiple attributes and SREG models with multiple attributes is that in the first case study, the two-way interaction considering all the attributes, while in the SREG model, it is possible to study the multivariate behavior of genotypes in relation to the mean environmental.


## 4. Conclusions

Based on three mode principal components analysis, it is possible to work simultaneously with multiple attributes in the adjustment of AMMI and SREG models, so we can do a multivariate study of the most important aspects of the genotype × environment interaction and the response of genotypes in different environmental conditions.When we have many combinations of levels of factors, the conditioned joint plot is a powerful tool to represent three arrays of markers and gives us information on the most important aspects of the response of genotypes. It was possible to detect the combinations of genotypes × locations × attributes causing the significant interaction and the differential response, that is, identify the best genotypes in view of its adaptability and performance.



---

*Source: 215843-2011-12-05.xml*