Abstract

Complete replicate block designs are fully efficient for treatment effects and are the designs of choice for many agricultural field experiments. For experiments with a large number of treatments, however, they may not provide good control of variability over the whole experimental area. Nested incomplete block designs with a single level of nesting can then improve ‘within-block’ homogeneity for moderate sized experiments. For very large designs, however, a single level of nesting may not be adequate and this paper discusses multi-level nesting with hierarchies of nested blocks. Multi-level nested block designs provide a range of block sizes which can improve ‘within-block’ homogeneity over a range of scales of measurement. We discuss design and analysis of multi-level block designs for hierarchies of nested blocks including designs with crossed block factors. We describe an R language package for multi-level block design and we exemplify the design and analysis of multi-level block designs by a simulation study of block designs for cereal variety trials in the UK. Finally, we re-analyse a single large row-and-column field trial for 272 spring barley varieties in 16 rows and 34 columns assuming an additional set of multi-level nested column blocks superimposed on the existing design. For each example, a multi-level mixed blocks analysis is compared with a spatial analysis based on hierarchical generalized additive (HGAM) models. We discuss the combined analysis of random blocks and HGAM smoothers in the same model.

Highlights

  • Comparative experiments in agriculture often involve the estimation of treatment effects against a background of high plot variability

  • The analysis shows the number of model terms, the mean estimated degrees of freedom (EDF), the mean AIC and the mean standard error of pairwise treatment differences (SED) averaged over simulations combination, the table shows the number of marginal model terms, the mean EDF, the mean AIC and the mean SED, calculated as described in the previous example

  • The mean pairwise SED for H3.1 was about 0.249 compared with about 0.260 for B3.3, which suggests that the best fitting HGAM model gave an improvement in precision of about 4% compared with the best fitting multi-level blocks model

Read more

Summary

INTRODUCTION

Comparative experiments in agriculture often involve the estimation of treatment effects against a background of high plot variability. Effective control of plot variability is essential for good treatment estimation and the most common method of control is the randomized complete blocks design. One commonly used method is the mixed model which assumes a mixture of fixed treatment effects and random block effects, see Pinheiro and Bates (2000) and Piepho et al (2003). This paper examines multi-level block designs for block and spatial effects in field experiments. We first make a simulation study of 3 replicates of 48 varieties in complete randomized blocks with four levels of nesting using simulation data based on a spatial correlation model from a study of 244 UK cereal variety trials by Patterson and Hunter (1983). Edmondson analysis so the assumption of fixed block effects at the design stage does not preclude the choice of an alternative model at the analysis stage. See Gilmour and Trinca (2006) Section 5.1.1. for further discussion of fixed versus random effects block models

MULTI-LEVEL BLOCK DESIGNS
OPTIMIZATION
OPTIMIZATION CRITERION
MIXED MODELS AND HGAM MODELS
MODEL FORMULATION
MODEL COMPARISON
MEAN PAIRWISE SED
SIMULATED EXAMPLE OF A NESTED BLOCKS DESIGN
SPATIAL DATA MODEL
SIMULATED DATA
EXAMPLE NESTED BLOCKS DESIGN
RANDOM BLOCKS ANALYSIS
INTERPRETATION
HGAM ANALYSIS
COMPARISON OF MIXED MODEL VERSUS HGAM ANALYSIS
SIMULATED EXAMPLE OF NESTING WITHIN A SET OF CROSSED BLOCKS
EXAMPLE CROSSED BLOCKS DESIGN
COMPARISON OF NESTED BLOCKS VERSUS HGAM ANALYSIS
A LARGE ROW-AND-COLUMN FIELD TRIAL EXAMPLE
SUPERIMPOSED COLUMN BLOCKS
HGAM MODELS
COMPARISON OF NESTED BLOCKS MODELS WITH HGAM MODELS
RESIDUAL COPLOTS
WIREFRAME PLOTS OF FIELD VARIABILITY
COMPARISONS OF MODEL PREDICTIONS
6.10. COMMENT ON DESIGN LAYOUT
Findings
DISCUSSION
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call