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)
You have computed the Pesek-Baker index for traits r traits
with the following desired genetic gains in actual units:
pbindex$Desired.Genetic.Gains
The estimated genotypic standard deviations for these traits are:
pbindex$Standard.Deviations
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, ]
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
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.