Computer and Biometric Services Unit International Center for Agricultural Research in the Dry Areas (ICARDA) P.O. Box 5466, Aleppo, Syria Online BioComputing Service ......... Date(D/M/Y) : 7 / 7 / 2004 Time : 9 hr 33 min Stability analysis of Multi-location multi-year trials using data on means ........This program computes stability parameters of Type 4 : C.S. Lin and M.R. Binns (1988). A method of analyzing cultivar x location x year experiments: a new stability parameter. Theo. Appl. Genet 76:425-430) ...................................................................................................... Notations GenoLev........ = Genotype number GenoMean....... = Genotype mean Slop........... = Slope on location index SeSlop......... = Standard error of Slop Prb1........... = Prob[t>|Slop-1|/SE] significance level DeviMS......... = Deviation mean-square PrDeviMS....... = Prob[Chisq >DeviMS/PoolMS] significance level MSYr_Loc....... = mean-square due to year within locations PrMS_Yr........ = Prob[Chisquare>MSYr_Loc*DF/PoolMS] TMP_CV1........ = (temporal) CVs of adjusted-means averaged over locations TMP_CV2........ = square-root of (average of (squares of temporal CVs of adjusted-means) ) ....................................................................................................... 1. Experiments and parameters of the designs ................................................... No. of locations/sites........... = 4 No. of years..................... = 4 No. of genotypes................. = 25 No. of plots in each experiments = 400 Variables analyzed................ : MnGYield,MnBYield ........................................................................... Variable = MnGYield ............................................................... ............................................................................ 2. Genotype - location means ............................................................... GenoLevF MeanL[1] MeanL[2] MeanL[3] MeanL[4] GenoMean 1 1192 1345 1517 969.5 1256 2 1253 1440 1433 920.0 1261 3 1180 1393 1450 998.8 1255 4 1218 1378 1391 931.8 1230 5 1101 1468 1614 838.2 1255 6 1165 1494 1411 1166.5 1309 7 1207 1637 1419 928.8 1298 8 1248 1386 1517 1063.5 1304 9 1303 1536 1575 928.2 1336 10 1143 1482 1451 890.8 1242 11 1161 1515 1521 958.2 1289 12 1251 1401 1652 997.2 1325 13 1179 1490 1549 941.0 1290 14 1282 1444 1507 1114.2 1337 15 1216 1402 1410 1043.0 1268 16 1122 1419 1452 976.8 1242 17 1182 1434 1360 1031.2 1252 18 974 1505 1468 833.2 1195 19 1136 1338 1356 969.5 1200 20 1179 1358 1317 1087.8 1235 21 1248 1427 1388 1075.2 1284 22 1038 1198 1392 967.8 1149 23 1205 1324 1214 749.5 1123 24 920 1152 1306 1027.2 1101 25 1064 1191 1418 877.8 1138 .................................................................................... 3. Pooled error mean square ............................................................... Pooled Error Mean-square (rep.-mean based) 16009 .................................................................................... 4. Stability indices (slope, deviation mean-sq, year (within location) mean-sq) and ranks ............................................................... GenoLevF GenoMean RankMean Slop SeSlop Prb1 DeviMS PrDeviMS RankDeviMS MSYr_Loc PrMS_Yr RankMSYr _Loc 1 1256 12 1.022 0.1793 0.4562 4708 0.3084 16 528404 0 17 2 1261 11 1.071 0.1821 0.3666 4855 0.2972 17 476191 0 9 3 1255 14 0.936 0.0301 0.0842 133 0.9674 3 453254 0 6 4 1230 19 0.946 0.1451 0.3720 3080 0.4632 14 548236 0 18 5 1255 13 1.583 0.1291 0.0229 2439 0.5436 11 523408 0 16 6 1309 4 0.685 0.2394 0.1596 8393 0.1228 18 567503 0 21 7 1298 6 1.277 0.3489 0.2555 17816 0.0117 24 431746 0 3 8 1304 5 0.859 0.1319 0.1981 2545 0.5295 12 504234 0 13 9 1336 2 1.324 0.1635 0.0931 3914 0.3761 15 497087 0 11 10 1242 17 1.257 0.1014 0.0632 1506 0.6865 10 523362 0 15 11 1289 8 1.251 0.0798 0.0440 932 0.7923 6 520521 0 14 12 1325 3 1.179 0.2743 0.2901 11010 0.0639 22 578630 0 22 13 1290 7 1.284 0.0225 0.0031 74 0.9816 2 500610 0 12 14 1337 1 0.794 0.0472 0.0244 326 0.9218 5 597771 0 24 15 1268 10 0.789 0.0388 0.0161 221 0.9464 4 550061 0 19 16 1242 16 1.039 0.0806 0.3375 951 0.7885 7 589457 0 23 17 1252 15 0.797 0.1368 0.1384 2740 0.5043 13 438872 0 5 18 1195 21 1.505 0.2498 0.0902 9134 0.1020 20 471670 0 8 19 1200 20 0.828 0.0188 0.0058 52 0.9872 1 402545 0 1 20 1235 18 0.552 0.0846 0.0169 1047 0.7698 8 464608 0 7 21 1284 9 0.708 0.0975 0.0478 1392 0.7062 9 551432 0 20 22 1149 22 0.782 0.2424 0.2319 8598 0.1167 19 437453 0 4 23 1123 24 0.999 0.4076 0.4991 24314 0.0023 25 406431 0 2 24 1101 25 0.568 0.3467 0.1697 17596 0.0123 23 671447 0 25 25 1138 23 0.965 0.2507 0.4503 9202 0.1003 21 477775 0 10 .................................................................................... 5. Stability indices (items in (4) and temporal %CV1,%CV2) and ranks ............................................................... GenoLevF GenoMean RankMean Slop DeviMS RankDeviMS MSYr_Loc RankMSYr_Loc TMP_CV1 RankCV1 TMP_CV2 RankCV2 1 1256 12 1.022 4708 16 528404 17 58.53 19 58.87 19 2 1261 11 1.071 4855 17 476191 9 54.44 7 54.57 7 3 1255 14 0.936 133 3 453254 6 52.60 2 53.40 3 4 1230 19 0.946 3080 14 548236 18 59.89 21 60.69 21 5 1255 13 1.583 2439 11 523408 16 59.44 20 60.50 20 6 1309 4 0.685 8393 18 567503 21 57.95 15 58.23 15 7 1298 6 1.277 17816 24 431746 3 50.48 1 50.86 1 8 1304 5 0.859 2545 12 504234 13 54.74 8 54.85 8 9 1336 2 1.324 3914 15 497087 11 53.22 5 53.62 5 10 1242 17 1.257 1506 10 523362 15 60.94 23 62.69 23 11 1289 8 1.251 932 6 520521 14 56.70 12 57.09 10 12 1325 3 1.179 11010 22 578630 22 57.49 14 57.62 14 13 1290 7 1.284 74 2 500610 12 56.53 11 57.33 12 14 1337 1 0.794 326 5 597771 24 58.08 17 58.46 17 15 1268 10 0.789 221 4 550061 19 58.30 18 58.84 18 16 1242 16 1.039 951 7 589457 23 63.49 24 64.46 24 17 1252 15 0.797 2740 13 438872 5 53.20 4 53.56 4 18 1195 21 1.505 9134 20 471670 8 52.82 3 53.35 2 19 1200 20 0.828 52 1 402545 1 53.34 6 53.72 6 20 1235 18 0.552 1047 8 464608 7 55.13 9 55.35 9 21 1284 9 0.708 1392 9 551432 20 57.25 13 57.26 11 22 1149 22 0.782 8598 19 437453 4 57.98 16 58.29 16 23 1123 24 0.999 24314 25 406431 2 55.59 10 57.38 13 24 1101 25 0.568 17596 23 671447 25 70.80 25 74.94 25 25 1138 23 0.965 9202 21 477775 10 60.26 22 61.31 22 .................................................................................... 6. Correlations between stability indices ............................................................... Corr GenoMean 1.0000 Slop 0.2262 1.0000 DeviMS -0.4896 0.0185 1.0000 MSYr_Loc 0.1895 -0.1724 -0.0940 1.0000 TMP_CV1 -0.4045 -0.2597 0.1109 0.7679 1.0000 TMP_CV2 -0.4615 -0.2420 0.1763 0.7312 0.9926 1.0000 GenoMean Slop DeviMS MSYr_Loc TMP_CV1 TMP_CV2 .................................................................................... 7. Graphs of stability indices with mean ............................................................... -+---------+---------+---------+---------+---------+---------+---------+---------+---------+---------+- I I I I I I 1.65 I I I I I I I 5 I I I I I 1.50 I 18 I I I I I I I I I I I 1.35 I I I 9 I I I I 13 7 I I 10 11 I I I 1.20 I I I 12 I I I I I I I I 2 I 1.05 I 16 I I 1 I I 23 I I 25 I I 4 I I 3 I 0.90 I I I I I 8 I I 19 I I 17 15 14 I I 22 I 0.75 I I I I I 21 I I 6 I I I I I 0.60 I I I 24 I I 20 I I I I I I I 0.45 I I -+---------+---------+---------+---------+---------+---------+---------+---------+---------+---------+- 1100.0 1125.0 1150.0 1175.0 1200.0 1225.0 1250.0 1275.0 1300.0 1325.0 1350.0 Slop v. GenoMean using factor GenoLevF -+---------+---------+---------+---------+---------+---------+---------+---------+---------+---------+- I I I I I 23 I 24000.0 I I I I I I I I I I I I 21000.0 I I I I I I I I I I I I 18000.0 I 7 I I 24 I I I I I I I I I 15000.0 I I I I I I I I I I I I 12000.0 I I I I I 12 I I I I I I I 9000.0 I 25 18 I I 22 6 I I I I I I I I I 6000.0 I I I I I 2 I I 1 I I 9 I I I 3000.0 I 4 I I 5 8 I I I I 10 21 I I 20 16 11 I I 14 I 0.0 I 19 3 15 13 I -+---------+---------+---------+---------+---------+---------+---------+---------+---------+---------+- 1100.0 1125.0 1150.0 1175.0 1200.0 1225.0 1250.0 1275.0 1300.0 1325.0 1350.0 DeviMS v. GenoMean using factor GenoLevF Points coinciding with 5 17 -+---------+---------+---------+---------+---------+---------+---------+---------+---------+---------+- I I I I I I 720000.0 I I I I I I I I I I I I 680000.0 I I I 24 I I I I I I I I I 640000.0 I I I I I I I I I I I I 600000.0 I 14 I I I I 16 I I 12 I I I I 6 I 560000.0 I I I 15 21 I I 4 I I I I I I 10 1 I 520000.0 I 11 I I I I 8 I I 13 9 I I I I I 480000.0 I 25 I I 18 2 I I 20 I I I I 3 I I I 440000.0 I 22 17 I I 7 I I I I I I I I 23 I 400000.0 I 19 I -+---------+---------+---------+---------+---------+---------+---------+---------+---------+---------+- 1100.0 1125.0 1150.0 1175.0 1200.0 1225.0 1250.0 1275.0 1300.0 1325.0 1350.0 MSYr_Loc v. GenoMean using factor GenoLevF Points coinciding with 1 5 -+---------+---------+---------+---------+---------+---------+---------+---------+---------+---------+- I I I 24 I I I 70.0 I I I I I I I I I I I I 67.5 I I I I I I I I I I I I 65.0 I I I I I I I I I 16 I I I 62.5 I I I I I I I I I 10 I I 25 I 60.0 I 4 I I 5 I I I I I I 1 15 I I 22 6 14 I 57.5 I 12 I I 21 I I 11 I I I I I I 23 I 55.0 I 20 I I 2 8 I I I I I I 19 17 9 I I 18 I 52.5 I 3 I I I I I I I I I I 7 I 50.0 I I -+---------+---------+---------+---------+---------+---------+---------+---------+---------+---------+- 1100.0 1125.0 1150.0 1175.0 1200.0 1225.0 1250.0 1275.0 1300.0 1325.0 1350.0 TMP_CV1 v. GenoMean using factor GenoLevF Points coinciding with 11 13 ........................................................................... Variable = MnBYield ............................................................... ............................................................................ 2. Genotype - location means ............................................................... GenoLevF MeanL[1] MeanL[2] MeanL[3] MeanL[4] GenoMean 1 5365 6145 6140 4660 5578 2 5261 6172 6798 4727 5740 3 5115 5827 6066 4858 5467 4 4999 5868 6478 4429 5444 5 4920 5990 7067 4880 5714 6 5280 6517 6455 5351 5901 7 5004 6126 7080 4251 5615 8 5008 5623 6573 4978 5545 9 5038 6111 6534 4676 5590 10 5222 6500 6322 4719 5691 11 4595 5502 6777 4749 5406 12 4784 5798 6814 4729 5531 13 5384 6229 6437 5177 5807 14 5290 6678 6408 4882 5815 15 5075 5894 7042 4852 5716 16 4796 5786 6676 4884 5536 17 5306 6110 6526 4584 5631 18 4942 6444 7016 4602 5751 19 5548 6122 6628 4512 5703 20 5385 5926 7378 4954 5911 21 5716 6662 6543 4744 5916 22 4886 5852 6618 3756 5278 23 5318 6406 6106 4242 5518 24 5366 5920 6152 5460 5724 25 5649 5416 6286 4756 5527 .................................................................................... 3. Pooled error mean square ............................................................... Pooled Error Mean-square (rep.-mean based) 156539 .................................................................................... 4. Stability indices (slope, deviation mean-sq, year (within location) mean-sq) and ranks ............................................................... GenoLevF GenoMean RankMean Slop SeSlop Prb1 DeviMS PrDeviMS RankDeviMS MSYr_Loc PrMS_Yr RankMSYr _Loc 1 5578 16 0.801 0.1921 0.2043 78685 0.1339 11 1754666 0 1 2 5740 7 1.095 0.0265 0.0350 1502 0.9623 1 2911140 0 11 3 5467 22 0.676 0.0455 0.0095 4408 0.8935 3 2303284 0 5 4 5444 23 1.076 0.0394 0.0962 3310 0.9189 2 3889640 0 22 5 5714 10 1.190 0.2229 0.2415 105988 0.0666 18 2513469 0 7 6 5901 3 0.741 0.2184 0.1787 101756 0.0743 17 2956142 0 13 7 5615 14 1.471 0.0748 0.0122 11948 0.7369 5 3329335 0 17 8 5545 17 0.835 0.2083 0.2561 92516 0.0940 14 3309748 0 16 9 5590 15 1.036 0.0560 0.2943 6699 0.8427 4 3378326 0 18 10 5691 12 0.959 0.2445 0.4413 127529 0.0384 19 2812872 0 10 11 5406 24 1.091 0.3206 0.4016 219250 0.0037 24 3592393 0 20 12 5531 19 1.133 0.2059 0.2917 90463 0.0991 13 4068850 0 24 13 5807 5 0.727 0.0756 0.0344 12185 0.7325 6 3222497 0 14 14 5815 4 0.949 0.2751 0.4345 161409 0.0162 20 3541312 0 19 15 5716 9 1.137 0.2097 0.2903 93816 0.0910 15 3682667 0 21 16 5536 18 1.000 0.2168 0.4996 100224 0.0772 16 3249502 0 15 17 5631 13 1.012 0.1061 0.4606 23992 0.5417 8 2634407 0 9 18 5751 6 1.370 0.1074 0.0374 24606 0.5333 9 2214580 0 4 19 5703 11 1.037 0.2054 0.4371 90009 0.1003 12 2024144 0 2 20 5911 2 1.182 0.2919 0.2980 181762 0.0096 22 3954167 0 23 21 5916 1 0.970 0.2886 0.4637 177608 0.0107 21 2553985 0 8 22 5278 25 1.446 0.1761 0.0634 66143 0.1845 10 2926124 0 12 23 5518 21 1.025 0.3623 0.4754 279968 0.0008 25 2415420 0 6 24 5724 8 0.421 0.0970 0.0135 20071 0.5988 7 5042684 0 25 25 5527 20 0.618 0.2996 0.1651 191428 0.0075 23 2175136 0 3 .................................................................................... 5. Stability indices (items in (4) and temporal %CV1,%CV2) and ranks ............................................................... GenoLevF GenoMean RankMean Slop DeviMS RankDeviMS MSYr_Loc RankMSYr_Loc TMP_CV1 RankCV1 TMP_CV2 RankCV2 1 5578 16 0.801 78685 11 1754666 1 22.99 2 24.18 1 2 5740 7 1.095 1502 1 2911140 11 27.27 8 29.11 9 3 5467 22 0.676 4408 3 2303284 5 26.23 6 27.75 6 4 5444 23 1.076 3310 2 3889640 22 34.26 22 36.24 22 5 5714 10 1.190 105988 18 2513469 7 27.90 10 28.55 8 6 5901 3 0.741 101756 17 2956142 13 29.34 11 29.53 11 7 5615 14 1.471 11948 5 3329335 17 29.73 12 32.02 15 8 5545 17 0.835 92516 14 3309748 16 32.81 21 33.86 21 9 5590 15 1.036 6699 4 3378326 18 31.64 17 32.47 16 10 5691 12 0.959 127529 19 2812872 10 29.86 13 30.47 12 11 5406 24 1.091 219250 24 3592393 20 34.63 23 37.10 24 12 5531 19 1.133 90463 13 4068850 24 36.44 24 36.61 23 13 5807 5 0.727 12185 6 3222497 14 30.78 14 31.01 13 14 5815 4 0.949 161409 20 3541312 19 31.50 16 32.95 17 15 5716 9 1.137 93816 15 3682667 21 32.10 18 33.16 19 16 5536 18 1.000 100224 16 3249502 15 32.14 19 33.59 20 17 5631 13 1.012 23992 8 2634407 9 26.80 7 28.04 7 18 5751 6 1.370 24606 9 2214580 4 24.15 3 25.36 3 19 5703 11 1.037 90009 12 2024144 2 22.44 1 24.28 2 20 5911 2 1.182 181762 22 3954167 23 32.23 20 33.04 18 21 5916 1 0.970 177608 21 2553985 8 25.16 5 26.20 5 22 5278 25 1.446 66143 10 2926124 12 30.86 15 31.91 14 23 5518 21 1.025 279968 25 2415420 6 27.58 9 29.13 10 24 5724 8 0.421 20071 7 5042684 25 38.63 25 40.25 25 25 5527 20 0.618 191428 23 2175136 3 24.91 4 25.63 4 .................................................................................... 6. Correlations between stability indices ............................................................... Corr GenoMean 1.0000 Slop -0.1592 1.0000 DeviMS 0.0483 -0.0256 1.0000 MSYr_Loc 0.0521 -0.0616 -0.1263 1.0000 TMP_CV1 -0.1638 -0.0580 -0.0554 0.9478 1.0000 TMP_CV2 -0.2117 -0.0267 -0.0610 0.9461 0.9892 1.0000 GenoMean Slop DeviMS MSYr_Loc TMP_CV1 TMP_CV2 .................................................................................... 7. Graphs of stability indices with mean ............................................................... -+---------+---------+---------+---------+---------+---------+---------+---------+---------+---------+- I I I I I I 1.50 I I I 7 I I 22 I I I I I I 18 I 1.35 I I I I I I I I I I I I 1.20 I 5 I I 2 I I I I 12 15 I I 11 2 I I 4 I 1.05 I I I 23 9 19 I I 16 17 I I 2I I 10 14 I I I 0.90 I I I I I I I 8 I I 1 I I I 0.75 I 6 I I 13 I I I I 3 I I I I 25 I 0.60 I I I I I I I I I I I I 0.45 I I I 24 I I I I I I I I I 0.30 I I -+---------+---------+---------+---------+---------+---------+---------+---------+---------+---------+- 5120.0 5200.0 5280.0 5360.0 5440.0 5520.0 5600.0 5680.0 5760.0 5840.0 5920.0 Slop v. GenoMean using factor GenoLevF -+---------+---------+---------+---------+---------+---------+---------+---------+---------+---------+- I I I I I I 320000.0 I I I I I I I I I I I I 280000.0 I 23 I I I I I I I I I I I 240000.0 I I I I I I I 11 I I I I I 200000.0 I I I 25 I I I I 22I I I I I 160000.0 I 14 I I I I I I I I I I 10 I 120000.0 I I I I I 5 I I 16 6 I I 128 15 I I I 80000.0 I 1 I I I I 22 I I I I I I I 40000.0 I I I I I 17 18 I I 24 I I 7 13 I I 3 9 I 0.0 I 4 2 I -+---------+---------+---------+---------+---------+---------+---------+---------+---------+---------+- 5120.0 5200.0 5280.0 5360.0 5440.0 5520.0 5600.0 5680.0 5760.0 5840.0 5920.0 DeviMS v. GenoMean using factor GenoLevF Points coinciding with 15 19 -+---------+---------+---------+---------+---------+---------+---------+---------+---------+---------+- I I I I I I 5500000.0 I I I I I I I I I I I 24 I 5000000.0 I I I I I I I I I I I I 4500000.0 I I I I I I I I I I I 12 I 4000000.0 I I I 4 2 I I I I I I 15 I I 11 I 3500000.0 I 14 I I 9 I I 8 7 I I 16 13 I I I I I 3000000.0 I I I 22 2 6 I I 10 I I I I 17 I I 2I 2500000.0 I 5 I I 23 I I 3 I I 18 I I 25 I I I 2000000.0 I 19 I I I I I I 1 I I I I I 1500000.0 I I -+---------+---------+---------+---------+---------+---------+---------+---------+---------+---------+- 5120.0 5200.0 5280.0 5360.0 5440.0 5520.0 5600.0 5680.0 5760.0 5840.0 5920.0 MSYr_Loc v. GenoMean using factor GenoLevF -+---------+---------+---------+---------+---------+---------+---------+---------+---------+---------+- I I I 24 I I I 38.0 I I I I I I I I I I I 12 I 36.0 I I I I I I I I I 11 I I 4 I 34.0 I I I I I I I I I 8 I I 2 I 32.0 I 16 15 I I 9 I I 14 I I 22 I I 13 I I I 30.0 I 10 I I 7 I I 6 I I I I I I I 28.0 I 5 I I 23 I I 2 I I I I 17 I I 3 I 26.0 I I I I I I I 25 2I I I I I 24.0 I 18 I I I I I I 1 I I I I 19 I 22.0 I I -+---------+---------+---------+---------+---------+---------+---------+---------+---------+---------+- 5120.0 5200.0 5280.0 5360.0 5440.0 5520.0 5600.0 5680.0 5760.0 5840.0 5920.0 TMP_CV1 v. GenoMean using factor GenoLevF