Learn R Programming

sommer (version 3.1)

blockerT: Applying transversal postblocking in a field

Description

The blockerT function takes a dataframe that contains columns indicating x and y coordinates and use them to create transversal blocking with different numbers or row width. For example, a square field with no blocks after applying the function will have 2 new columns (forward and backwards) for the transversal blocking desired (angle degree and #rows specified).

Usage

blockerT(dat, d= 60, nr= 5, rows="ROW",
           ranges="RANGE", by=NULL)

Arguments

dat

a dataframe with 2 obligatory columns; rows and ranges which can have different names and can be matched with the following arguments.

d

angle degree used to draw the transversal blocks.

nr

number of rows width of the block.

rows

the name of the numeric column that indicates the x coordinate direction in the field.

ranges

the name of the numeric column that indicates the y coordinate other direction in the field.

by

optional argument to indicate the name of the column of the dataframe x that indicates the environments so the field is filled by environment.

Value

$fin

a new dataframe with a 2 new columns for transversal blocking forward and backwards.

References

Covarrubias-Pazaran G (2016) Genome assisted prediction of quantitative traits using the R package sommer. PLoS ONE 11(6): doi:10.1371/journal.pone.0156744

Fikret Isik. 2009. Analysis of Diallel Mating Designs. North Carolina State University, Raleigh, USA.

See Also

The core functions of the package mmer and mmer2

Examples

Run this code
# NOT RUN {
  
####=========================================####
#### For CRAN time limitations most lines in the 
#### examples are silenced with one '#' mark, 
#### remove them and run the examples using
#### command + shift + C |OR| control + shift + C
####=========================================####
data(CPdata)
#### look at the data
head(CPpheno)
#### fill the design
gg <- fill.design(x=CPpheno, rows="Row",ranges="Col")
head(gg)
#### apply the postblocking 
gg2 <- blockerT(dat=gg, d= 60, nr= 5, rows="Col",ranges="Row")
head(gg2)
#### see the new blocking
# lattice::levelplot(TBLOCKB~Row*Col|FIELDINST, data=gg2)
# lattice::levelplot(TBLOCKF~Row*Col|FIELDINST, data=gg2)
#### now you can use them in your mixed models


# }

Run the code above in your browser using DataLab