The purpose of this study was to conduct selection, genetic parameter estimation, and prediction of genetic values for 18 S 1 families of guava trees using mixed model methodology and simultaneous selection of traits by means of the additive selection index, multiplicative selection index, and mean rank adapted from Mulamba. All families analyzed were obtained by means of self-fertilization of superior genotypes (full siblings) from the genetic breeding program of guava trees at the Universidade Estadual do Norte Fluminense. An experimental randomized block design with 18 S 1 families, three replicates, and ten plants per plot was used. A total of 540 genotypes (individual plants) of guava tree were evaluated. Genetic parameter estimation and selection of the best genotypes based on the genetic value were performed using the statistical procedure, from the Selegen-REML/BLUP program. The analyses of the additive selection index, multiplicative selection index, and the sum of rank adapted from Mulamba were also performed under the Selegen program. During the evaluation by the individual BLUPs, families 1, 12, 4, 6, and 8 contributed to most of the genotypes selected for the traits under evaluation, suggesting their significant potential to generate high quality and high yield genotypes. In the selection indexes via mixed models, the multiplicative index showed higher values for genetic gains (74 %), followed by the mean rank index adapted from Mulamba (19 %), and the additive index (2 %).