Consider the statistical model for *Complete Randomized design (CRD)*
y_{ij} = \mu + \tau_i + \epsilon_{ij}, \quad i=1,\ldots ,k, \quad j=1,\ldots, n_i,
where $y_{ij}$ is a random variable that represents the response obtained on the $j$th
observation of the ith treatment, and \epsilon_{ij}...