In CIP-RIU/hidap: High Interactive Data Analysis Platform for Plant breeding

library(knitr)

traits <- params$traits
geno <- params$geno
env <- params$env
rep <- params$rep
data <- params$data
maxp <- params$maxp

1. Model specification and data description

There is data from r nlevels(as.factor(data[, env])) environments and r nlevels(as.factor(data[, geno])) genotypes. In each environment the genotypes were evaluated using a randomized complete block design with r nlevels(as.factor(data[, rep])) blocks. The statistical model is $$ y_{ijk} = \mu + \alpha_i + \beta_j + (\alpha\beta){ij} + \gamma{k(j)} + \epsilon_{ijk} $$ where

• $y_{ijk}$ is the observed response with genotype $i$, environment $j$, and block $k$ nested in environment $j$.
• $\mu$ is the mean response over all genotypes, environments and blocks.
• $\alpha_i$ is the effect for genotype $i$.
• $\beta_j$ is the effect for environment $j$.
• $(\alpha\beta)_{ij}$ is the interaction effect between genotype $i$ and environment $j$.
• $\gamma_{k(j)}$ is the effect of block $k$ nested in environment $j$.
• $\epsilon_{ijk}$ is the error term.

In this model we assume that the errors are independent and have a normal distribution with common variance, that is, $\epsilon_{ijk} \sim N(0,\sigma_{\epsilon}^2)$.

out <- NULL
for (i in 1:length(traits))
  out <- c(out, knit_expand('child_met.Rmd'))

r paste(knit(text = out), collapse = '\n')

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