h_cullis: Cullis heritability for lme4 models
In AparicioJohan/agriutilities: Utilities for Data Analysis in Agriculture

h_cullis

R Documentation

Cullis heritability for lme4 models

Description

Cullis heritability for lme4 models

Usage

h_cullis(model, genotype, re_MME = FALSE)

Arguments

`model`	Object of class `lmer`.
`genotype`	A character string indicating the column in data that contains genotypes.
`re_MME`	A logical value to ask if we want to reconstruct the mixed models equations to estimate the Cullis heritability. (`FALSE` by default)

Value

A numerical value of the Cullis heritability estimate. If re_MME is TRUE, a list with matrices of the mixed models equations is returned.

Author(s)

Paul Schmidt, Johan Aparicio.

Examples


library(lme4)
library(agridat)
library(agriutilities)
dat <- john.alpha
g.ran <- lmer(
  formula = yield ~ rep + (1 | gen) + (1 | rep:block),
  data = dat
)
h_cullis(model = g.ran, genotype = "gen")

AparicioJohan/agriutilities documentation built on Jan. 20, 2025, 7:53 a.m.