h_cullis_spt | R Documentation |
Calculate Cullis heritabilities from SpATS objects
h_cullis_spt(model)
model |
an object of class SpATS as produced by SpATS() |
A data frame. The data frame has the following components
trait
: Character string with the trait being analyzed
H2Cullis
: Generalized heritability proposed by Cullis (2006)
H2Oakey
: Generalized heritability proposed by Oakey (2006)
reBLUP_avg
: Average BLUP reliability
vdBLUP_avg
: Average pairwise prediction error variance of genotype effects
PEV_avg
: Average predictive error variance (PEV) of genotype effects
var_G
: Genotypic Variance
Johan Aparicio
Cullis, B. R., Smith, A. B., & Coombes, N. E. (2006). On the design of early generation variety trials with correlated data. Journal of agricultural, biological, and environmental statistics, 11, 381-393.
Oakey, H., A. Verbyla, W. Pitchford, B. Cullis, and H. Kuchel (2006). Joint modeling of additive and non-additive genetic line effects in single field trials. Theoretical and Applied Genetics, 113, 809 - 819.
library(SpATS)
library(agriutilities)
data(wheatdata)
wheatdata$R <- as.factor(wheatdata$row)
wheatdata$C <- as.factor(wheatdata$col)
m1 <- SpATS(
response = "yield",
spatial = ~ PSANOVA(col, row, nseg = c(10, 20), nest.div = 2),
genotype = "geno",
genotype.as.random = TRUE,
fixed = ~ colcode + rowcode,
random = ~ R + C,
data = wheatdata,
control = list(tolerance = 1e-03, monitoring = 0)
)
h_cullis_spt(m1)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.