calc_Bhat: Calculate estimated allele effects, B matrix

In fboehm/qtl2pleio: Testing Pleiotropy in Multiparental Populations

Description

Calculate estimated allele effects, B matrix

Usage

calc_Bhat(X, Sigma_inv, Y)

Arguments

X	dn by df block-diagonal design matrix that incorporates genetic info for d markers. Note that we can use the same marker data twice.
Sigma_inv	dn by dn inverse covariance matrix, often composed as the inverse of K \otimes V_g + I_n \otimes V_e
Y	dn by 1 matrix, ie, a column vector, of d phenotypes' measurements

Value

a df by 1 matrix of GLS-estimated allele effects

Examples

X1 <- as.matrix(rbinom(n = 100, size = 1, prob = 1 / 2))
X <- gemma2::stagger_mats(X1, X1)
Sigma_inv <- diag(200)
Y <- runif(200)
calc_Bhat(X, Sigma_inv, Y)

fboehm/qtl2pleio documentation built on July 17, 2021, 11:34 p.m.