source("https://raw.githubusercontent.com/Flavjack/inti/master/pkgdown/favicon/docs.r") knitr::opts_chunk$set(echo = TRUE)
Broad-sense heritability ($H^2$) is defined as the proportion of phenotypic variance that is attributable to an overall genetic variance for the genotype [@schmidt2019Estimating]. There are usually additional interpretations associated with $H^2$: (i) It is equivalent to the coefficient of determination of a linear regression of the unobservable genotypic value on the observed phenotype; (ii) It is also the squared correlation between predicted phenotypic value and genotypic value; and (iii) It represents the proportion of the selection differential ($S$) that can be realized as the response to selection ($R$) [@falconer2005Introduction].
There are two main reasons why heritability on an entry-mean basis is of interest in plant breeding [@schmidt2019Heritability]:
$$\Delta G=H^2S$$ Where:
In practice, most trials are conducted in a multienvironment trial (MET) presente unbalanced data as not all cultivars are tested at each environment or simply when plot data is lost or when the number of replicates at each location varies between genotypes [@schmidt2019Estimating]. However, the standard method for estimating heritability implicitly assumes balanced data, independent genotype effects, and homogeneous variances.
According @schmidt2019Heritability, the variance components could be calculated in two ways:
For the two stage approach, in the first stage each experiment is analyzed individually according their experiment design (Lattice, CRBD, etc) [@zystro2018Alternative]. And for the second stage environments are denotes a year-by-location interaction. This approach assumes a single variance for genotype-by-environment interactions (GxE), even when multiple locations were tested across multiple years [@buntaran2020Cross].
$$y_{ikt}=\mu\ +\ G_i+E_t+GxE_{it}+\varepsilon_{ikt}$$
$$\sigma_p^2=\sigma_g^2+\frac{\sigma_{g\cdot e}^2}{n_e}+\frac{\sigma_{\varepsilon}^2}{n_e\cdot n_r}$$
For the one stage approach only one model is used for the MET analysis. The environmental effects are included via separate year, and location main interaction effects.
$$y_{ikt}=\mu+G_i+Y_m+E_q+YxE_{mq}+GxY_{im}+GxE_{iq}+GxYxE_{imq}+\varepsilon_{ikmq}$$
$$\sigma_p^2=\sigma_g^2+\frac{\sigma_{g\cdot e}^2}{n_e}+\frac{\sigma_{g\cdot y}^2}{n_y}+\frac{\sigma_{g\cdot y\cdot e}^2}{n_y\cdot n_e}+\ \frac{\sigma_{\epsilon}^2}{n_e\cdot n_y\cdot n_r}$$
tibble( Standart = "$H^2=\\frac{\\sigma_g^2}{\\sigma_p^2}=\\frac{\\Delta G}{S}$", Cullis = "$H_{Cullis}^2=1-\\frac{\\overline{V}_{\\Delta..}^{^{BLUP}}}{2\\cdot\\sigma_g^2}$", Piepho = "$H_{Piepho}^2=\\frac{\\sigma_g^2}{\\sigma_g^2+\\frac{\\overline{V}_{\\Delta..}^{BLUE}}{2}}$" ) %>% kable(caption = "Differentes heritability calculation")
For calculate the standard heritability in MET experiments the number of location and replication should be include manually in the function H2cal()
. In the case of difference number of replication in each experiments, take the maximum value (often done in practice) [@schmidt2019Estimating].
For remove the outliers the function implemented is the Method 4 used for @bernal-vasquez2016Outliera: Bonferroni-Holm using re-scaled MAD for standardizing residuals (BH-MADR).
library(inti)
dt <- potato hr <- H2cal(data = dt , trait = "stemdw" , gen.name = "geno" , rep.n = 5 , fixed.model = "0 + (1|bloque) + geno" , random.model = "1 + (1|bloque) + (1|geno)" , emmeans = TRUE , plot_diag = TRUE , outliers.rm = TRUE )
hr$model %>% summary()
hr$tabsmr %>% kable(caption = "Variance component table")
hr$blues %>% kable(caption = "BLUEs")
hr$blups %>% kable(caption = "BLUPs")
hr$outliers$fixed %>% kable(caption = "Outliers fixed model")
hr$outliers$random %>% kable(caption = "Outliers random model")
::: {#refs} :::
if(!file.exists("files/pkgs.bib")){write_bib(c(.packages()),'files/pkgs.bib')}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.