  In this study, we addressed the problem of genome-wide prediction accounting
for partial correlation of marker effects when the partial correlation
structure, or equivalently, the pattern of zeros of the precision matrix is
unknown. This problem requires estimating the partial correlation structure of
marker effects, that is, learning the pattern of zeros of the corresponding
precision matrix, estimating its non-null entries, and incorporating the
inferred concentration matrix in the prediction of marker allelic effects. To
this end, we developed a set of statistical methods based on Gaussian
concentration graph models (GCGM) and Gaussian directed acyclic graph models
(GDAGM) that adapt the existing theory to perform covariance model selection
(GCGM) or DAG selection (GDAGM) to genome-wide prediction. Bayesian and
frequentist approaches were formulated. Our frequentist formulations combined
some existing methods with the EM algorithm and were termed Glasso-EM,
CONCORD-EM and CSCS-EM, whereas our Bayesian formulations corresponded to
hierarchical models termed Bayes G-Sel and Bayes DAG-Sel. Results from a
simulation study showed that our methods can accurately recover the partial
correlation structure and estimate the precision matrix. Methods CONCORD-EM and
Bayes G-Sel had an outstanding performance in estimating the partial
correlation structure and a method based on CONCORD-EM yielded the most
accurate estimates of the precision matrix. Our methods can be used as
predictive machines and as tools to learn about the covariation of effects of
pairs of loci on a given phenotype conditioned on the effects of all the other
loci considered in the model. Therefore, they are useful tools to learn about
the underlying biology of a given trait because they help to understand
relationships between different regions of the genome in terms of the partial
correlations of their effects on that trait.
