Estimate breeding values - Adapted function from QTLpoly
breeding_values(qtl_info, probs, selected_mks, blups, beta.hat, pos)
data.frame containing breeding values
data.frame with: LG - linkage group ID; Pos - position in linkage map (cM); Pheno - phenotype ID; Pos_lower - lower position of confidence interval; Pos_upper - upper position of the confidence interval; Pval - QTL p-value; h2 - herdability
data.frame with first column (named `ind`) as individuals ID and next columns named with markers ID and containing the genotype probability at each marker
data.frame with: LG - linkage group ID; mk - marker ID; pos - position in linkage map (cM)
data.frame with: haplo - haplotype ID; pheno - phenotype ID; qtl - QTL ID; u.hat - QTL estimated BLUPs
data.frame with: pheno - phenotype ID; beta.hat - estimated beta
selected QTL position (cM)