Introduction
The goal of a plant breeding program is to produce genotypes which are, in some sense, optimum for the conditions under which they will be grown. Many models have been developed to measure the stability parameters and partitioning of variation due to GxE interactions (Finlay and Wilkinson, 1963; Eberhart and Russel, 1966). The model proposed by Eberhart and Russel (1966) is considered more appropriate to interpret the stability statistics and is more commonly used for stability studies in crops. Lin et al. (1986) has established three concepts of stability. Type 1 stability measures (genotype mean square=S^{2}_{i} and genotypic coefficient of variation=CV_{i}) are those which measure the variation within a genotype across environment. These statistics do not depend upon the other genotypes which might be included in the trials. Thus they provide very broad based inference and are commonly avoided for making the final decision. Type 2 stability measures (ecovalence=W^{2}_{i} and Shukla’s stability variance (σ^{2}_{i}) which basically measure the deviation of the individual genotype from the location means of all genotypes in test. Type 3 (regression slope = b_{i}) stability is calculated by the residual mean square from the regression of individual cultivar yields on an environmental index (Eberhart and Russell, 1966). In this method, the slope of regression provides an indication of regions of adaptability as well as stability. It also indicates the cultivar response to the predictable component of the environment. Naazar et al. (2002) indicated that GxE interaction was highly significant. A top yielding genotype SLM046 was found a stable cultivar for grain yield. Genotypes Regent, Cobra and A.W. were found suitable for favorable environments, Whereas PF7045/91 and Eureka could be recommended for poor environments. The objective of the present study was to evaluate and identify the promising genotypes which could be considered adaptable under broad environmental conditions and can be grouped together.
Materials and Methods
The data comes from normal duration replicated NUWYT trials for 20002001. Two seeding dates (Normal and Late) were used for trials at 30 locations in Pakistan. Twenty candidate genotypes submitted by various wheat breeders of the country were sown with a local check at each location. It was recommended that each plot consists of 6 rows, 5 m long and 30 cm apart and each genotype is grown in four replicates. A randomization plan was given by wheat program NARC to lay out experiments at different locations. Data was collected at maturity and sent to wheat program, NARC for further evaluation.
The methods of analysis used for this data set were done in the following sequence. The combined analysis of variance of yield data over all environments, using GenotypeEnvironment interaction data for stability analysis with conventional and unconventional approaches. This was supplemented by graphical representation of the data using GEBI software which uses cluster analysis to form groups of genotypes which are similar to one another based on response pattern towards grain yield.
Results and Discussion
The combined analysis of variance revealed significant GxE indicating the influence
of environments on the yield performance of genotypes (Gomes and Gomes, 1984)
(Table 1).
Since in the combined analysis GxE interaction was highly significant therefore a stability analysis based on location index was in order. However rather than depending only on the Eberhart and Russell (1966) approach a more comprehensive set of stability analysis as summarized by Lin et al. (1986) was employed. The results of stability analysis are summarized in Table 2.
The more frequently used method of stability analysis involves comparing the
competing genotypes with respect to other genotypes by regressing yields attained
by each genotype at different locations over an environmental index which is
based on average yield of all genotypes at each location. The resulting slope
for each genotype can then be considered as a measure of stability using a unit
slope to be stable standard. The deviations from the regression line are also
considered as they provide a good measure of fit. However Finlay and Wilkinson
(1963) pointed out that the slope provides an indication of regions of adaptability
as well as stability. Genotypes that have a slope significantly greater than
one are specifically adapted to high yielding environments.
Table 1: 
Pooled analysis of variance of the gain yield data of wheat
genotypes 

** Statistically significant at 1% level 
Table 2: 
Stability Statistics for twenty genotype of wheat grown at
thirty locations in Pakistan, 20002001 

On the other hand, genotypes with a slope less than one are insensitive to
change in environment and are, therefore better adapted to poor environments.
The competing genotypes in the seeding trials on the basis of regression model
can be divided into three distinct groups.
The first group which is categorized as stable group has b values ranging from
0.97 to 1.04 include 7 genotypes V7004, V97024, 92T009, 97B2236, Inqalab,
V97112 and IBW96405. However some of these genotypes are weak on other measures
of stability which makes it dangerous to recommend these on the basis of b value
only. The genotypes that can be safely termed as stable are V7004 and V97024.
The other genotypes in the stable group have high values of ecovalence so we
should use these with some caution.

Fig. 1: 
Dendogram showing clustering of genotypes 

Fig. 2: 
Dendogram showing clustering of different locations of experiment 
The other group consists of 6 genotypes which have values significantly greater
than one are DN16, V97046, PR17, D97603, 97B2210 and local check. The last
group of seven genotypes comprises of those genotypes which by and large are
insensitive to environmental change. These genotypes are V7005, SD1200/14,
V8964, 91BT0105, V97052, SI91195 and PR73 (Table 2).
The last technique that is used to process this data set is multivariate method of cluster analysis (Sneath and Sokal, 1973). The data of genotype x environment (G X E) tables of yields, clustering is used to simplify the data set by grouping the genotypes, over all environments, with similar response patterns of all yields. In a similar fashion growing of the environments, over all genotypes, with similar response pattern for all yields (Byth et al., 1976). The method used for hierarchical classification requires a measure of association (proximity measure) among the individuals and a fusion strategy. The proximity measure provides a measure of the distance or closeness in multidimensional space.

Fig. 3: 
Performance plots showing performance of different groups
of genotypes 
The pictorials include dendograms and performance plots. Two major contributors accounted for the overall variability in the genotype x environment data. Genotypes accounted for about 15.98% and GxE interaction contribution was 84.02%. The dendogram of genotypes is given in Fig. 1. Clearly five groups at fusion level of 15 are formed. Similarly environments are be grouped into seven clusters (Fig. 2).
Fig. 3. shows that consistently no group performed well over all environment groups. However performance of some groups was much better than the others. It is therefore necessary that performance plots are given due consideration while deciding on better genotypes. The genotypes in group14 performed consistently poor on all seven environmental group positions, but group15 and group9 genotypes performed poorly on only four out of seven. In this respect the genotypes in groups 12 and 13 performed much better as at only two environment group locations they performed at below par and at five environment group locations they performed better than average. This makes eight genotypes and local check as the better genotypes than the others in the group for these trials and can be recommended for wider adaptation. The names of these genotypes are V97024, V97046, V97112, V97052, D97603, IBW96405, PR70, 97B2210 and LCHECK. When compared with the recommendations made for stability using parametric approach only three genotypes namely V97024, V97112 and IBW96405 satisfied criterion given by Eberhart and Russell (1966), but only V97024 was declared as stable considering other parameters of stability.