Abstract Root system architecture (RSA) is a crucial factor in resource acquisition and plant productivity. Roots are difficult to phenotype in the field, thus new tools for predicting phenotype from genotype are particularly valuable for plant breeders aiming to improve RSA. This study identifies quantitative trait loci (QTLs) for RSA and agronomic traits in a rice (Oryza sativa) recombinant inbred line (RIL) population derived from parents with contrasting RSA traits (PI312777 × Katy). The lines were phenotyped for agronomic traits in the field, and separately grown as seedlings on agar plates which were imaged to extract RSA trait measurements. QTLs were discovered from conventional linkage analysis and from a machine learning approach using a Bayesian network (BN) consisting of genome-wide SNP data and phenotypic data. The genomic prediction abilities (GPAs) of multi-QTL models and the BN analysis were compared with the several standard genomic prediction (GP) methods. We found GPAs were improved using multitrait (BN) compared to single trait GP in traits with low to moderate heritability. Two groups of individuals were selected based on GPs and a modified rank sum index (GSRI) indicating their divergence across multiple RSA traits. Selections made on GPs did result in differences between the group means for numerous RSA. The ranking accuracy across RSA traits among the individual selected RILs ranged from 0.14 for root volume to 0.59 for lateral root tips. We conclude that the multitrait GP model using BN can in some cases improve the GPA of RSA and agronomic traits, and the GSRI approach is useful to simultaneously select for a desired set of RSA traits in a segregating population.