E. Worku Cluster analysis, correlation, path analysis, RILs, tef 115-131 9 (4) 2025 October, 2025 Forty-nine tef RILs with checks was conducted in Gondar Zuria district, Ethiopia to assess multivariate variability and association among traits using 7×7 simple lattice design. Analysis of variance revealed a highly significant and significant variation among the RILs except for thousand seed weight. Grain yield exhibited highly significant positive genotypic correlation with plant height, biomass yield and harvest index and highly significant positive phenotypic correlation with plant height, main shoot panicle length, number of primary branches per main shoot panicle, main shoot panicle grain yield, biomass yield and harvest index. The path coefficient analyses both at genotypic and phenotypic levels revealed harvest index and above ground biomass yield exert strong direct effect on grain yield; indicating that these traits could be used as selection criteria for yield improvement. Cluster analysis, grouped the tested genotypes into five clusters with an inter- and intra-cluster distance ranging from 12.79 to 85.14 and 7.82 to 15.29, respectively. Also cluster III had larger distance from Cluster II and Cluster V; the divergence of cluster III from these two clusters indicate existence of RILs to be crossed and selected. The first five principal components (PCs) accounted for 70.8% of the total genetic variation. Generally, this considerable genetic variability among the RILs confirmed the possibility of improving tef productivity.