calc_XHiY: Calculate XHiY

Description Usage Arguments Value Examples

View source: R/calc_XHiY.R

Description

Calculate XHiY

Usage

1
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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
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.