The selection orders and the three-stage procedure. In stage 1, we get the big scale selection order by NGVS. NGVS incorporates the prior knowledge of LD structure and gene expression network with the flexibility of model setting. In stage 2, we find the small scale selection orders within each marker set by GLasso. In stage 3 we combine the two selection orders together using P-values of single QTL regression and the three criteria discussed. Results in the picture were obtained from the adjacent interaction two-effect model (model assumption), full co-expression network of gene Ankrd24, Sfrp1 and Ergic1 (gene expression network) and marker sets formation by chromosomes (marker sets).