calc_XHiY: Calculate XHiY

In gemma2: GEMMA Multivariate Linear Mixed Model

Description

Calculate XHiY

Usage

calc_XHiY(eval, D_l, X, UltVehiY)

Arguments

eval	vector of eigenvalues from the decomposition of the relatedness matrix
D_l	vector of length d_size
X	design matrix
UltVehiY	a matrix

Value

numeric vector

Examples

readr::read_tsv(system.file("extdata",
"mouse100.pheno.txt",
package = "gemma2"),
col_names = FALSE) -> pheno
phe16 <- as.matrix(pheno[, c(1, 6)])
as.matrix(readr::read_tsv(system.file("extdata",
"mouse100.cXX.txt",
package = "gemma2"),
col_names = FALSE)[, 1:100]) -> kinship
eigen2(kinship) -> eout
eout$values -> eval
eout$vectors -> U
UltVehi <- matrix(c(0, -1.76769, -1.334414, 0),
nrow = 2,
byrow = FALSE) # from output of eigen_proc()
calc_XHiY(eval = eval,
D_l = c(0.9452233, 5.9792268),
          X = rep(1, 100) %*% U,
          UltVehiY = UltVehi %*% t(phe16) %*% U
          )

gemma2 documentation built on Oct. 24, 2020, 5:06 p.m.