calc_Bhat {qtl2pleio} | R Documentation |
Calculate estimated allele effects, B matrix
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 |
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)
[Package qtl2pleio version 1.4.3 Index]