Selection according to weighted breeding values

Usage

selectWBV(pop, n, QTNeffects)

Arguments

pop: (population object) Population of individuals to select (see:population)
n: (integer) number of individuals to select
QTNeffects: (numeric vector) effect of the genetic markers

Value

character vector of the selected individuals' names

References

Jannink, Jean-Luc. “Dynamics of Long-Term Genomic Selection.” Genetics Selection Evolution 42, no. 1 (December 2010). https://doi.org/10.1186/1297-9686-42-35.

Examples

mySpec <- specie$new(specName = "Statisticae exempli",
                     nChr = 10,
                     lchr = 1e6,
                     lchrCm = 100)
#> A new species has emerged: Statisticae exempli !
#> 
SNPs <- SNPinfo$new(SNPcoord = exampleData$snpCoord, specie = mySpec)
initPop <- createPop(geno = exampleData$genotypes,
                     SNPinfo = SNPs,
                     popName = "Initial population")
#> Create population: Initialisation...
#> Create population: Create individuals...
#> 
1%
2%
3%
4%
5%
6%
7%
8%
9%
10%
11%
12%
13%
14%
15%
16%
17%
18%
19%
20%
21%
22%
23%
24%
25%
26%
27%
28%
29%
30%
31%
32%
33%
34%
35%
36%
37%
38%
39%
40%
41%
42%
43%
44%
45%
46%
47%
48%
49%
50%
51%
52%
53%
54%
55%
56%
57%
58%
59%
60%
61%
62%
63%
64%
65%
66%
67%
68%
69%
70%
71%
72%
73%
74%
75%
76%
77%
78%
79%
80%
81%
82%
83%
84%
85%
86%
87%
88%
89%
90%
91%
92%
93%
94%
95%
96%
97%
98%
99%
100%
#> Create population: Create population object...
#> Create population: Add individuals...
#> 
1%
2%
3%
4%
5%
6%
7%
8%
9%
10%
11%
12%
13%
14%
15%
16%
17%
18%
19%
20%
21%
22%
23%
24%
25%
26%
27%
28%
29%
30%
31%
32%
33%
34%
35%
36%
37%
38%
39%
40%
41%
42%
43%
44%
45%
46%
47%
48%
49%
50%
51%
52%
53%
54%
55%
56%
57%
58%
59%
60%
61%
62%
63%
64%
65%
66%
67%
68%
69%
70%
71%
72%
73%
74%
75%
76%
77%
78%
79%
80%
81%
82%
83%
84%
85%
86%
87%
88%
89%
90%
91%
92%
93%
94%
95%
96%
97%
98%
99%
100%
#> A new population created: Initial population !
myTrait <- trait$new(name = "myTrait",
                     qtn = names(exampleData$snpEffects),
                     qtnEff = exampleData$snpEffects)
selectWBV(pop = initPop,
          QTNeffects = myTrait$qtnEff,
          n = 10)
#>  [1] "Coll0068" "Coll0079" "Coll0008" "Coll0020" "Coll0036" "Coll0049"
#>  [7] "Coll0016" "Coll0002" "Coll0083" "Coll0074"