alpha_cond_params: Compute conditional distribution of sires effect

View source: R/unbalanced_EM.R

alpha_cond_paramsR Documentation

Compute conditional distribution of sires effect

Description

Given the model

y[ij] | alpha[i] ~ N(alpha[i], Sigma[E]) alpha[i] ~ N(0, Sigma[A]),

computes the parameters of the Gaussian conditional distribution

alpha[i] | y[i].

Usage

alpha_cond_params(y_data, Sigma_E, Sigma_A, mu = rep(0, nrow(Sigma_E)))

Arguments

y_data

Observed data that inherits "fullsibdata"

Sigma_E

Between-individuals covariance matrix

Sigma_A

Between-sires covariance matrix

mu

Global mean

Value

A list of multivariate Gaussian parameters for the conditional distribution of the sires random effect. Each entry is a list with entries mean and cov, which are the mean vector and covariance matrix respectively.


damian-t-p/halfsibdesign documentation built on March 14, 2023, 4:55 a.m.