traits <- params$traits
geno <- params$geno
env <- params$env
rep <- params$rep
dfr <- params$dfr
means <- params$means
dgg <- params$dgg
units <- params$units
sf <- params$sf

pbindex <- pesekbaker(traits, geno, env, rep, dfr, means, dgg, units, sf)

Desired genetic gains

You have computed the Pesek-Baker index for traits r traits with the following desired genetic gains in actual units:

pbindex$Desired.Genetic.Gains

Standard deviation

The estimated genotypic standard deviations for these traits are:

pbindex$Standard.Deviations

The Pesek-Baker index

The index coefficients for the Pesek-Baker index are:

pbindex$Index.Coefficients

With this, the Pesek-Baker index for each genotype is shown below. As you see, genotype r (pbindex$Pesek.Baker.Index[sort(pbindex$Pesek.Baker.Index$PB.Index, decreasing = TRUE, index.return = TRUE)$ix, ])[1,1] is the one with the highest value.

pbindex$Pesek.Baker.Index[sort(pbindex$Pesek.Baker.Index$PB.Index, decreasing = TRUE,
                               index.return = TRUE)$ix, ]

The response to selection

For a selection fraction of r sf, the responses to selection in actual units are:

pbindex$Response.to.Selection

and in standardized units:

pbindex$Std.Response.to.Selection


CIP-RIU/hidap documentation built on April 30, 2021, 9:21 p.m.