score.calc: Calculate -log10(p) for single-SNP GWAS

In KosukeHamazaki/RAINBOW: Genome-Wide Association Study with SNP-Set Methods

Description

Calculate -log10(p) of each SNP by the Wald test.

Usage

score.calc(
  M.now,
  ZETA.now,
  y,
  X.now,
  Hinv,
  P3D = TRUE,
  optimizer = "nlminb",
  eigen.G = NULL,
  min.MAF = 0.02,
  count = TRUE
)

Arguments

M.now	A n.sample x n.mark genotype matrix where n.sample is sample size and n.mark is the number of markers.
ZETA.now	A list of variance (relationship) matrix (K; m \times m) and its design matrix (Z; n \times m) of random effects. You can use only one kernel matrix. For example, ZETA = list(A = list(Z = Z, K = K)) Please set names of list "Z" and "K"!
y	A n \times 1 vector. A vector of phenotypic values should be used. NA is allowed.
X.now	A n \times p matrix. You should assign mean vector (rep(1, n)) and covariates. NA is not allowed.
Hinv	The inverse of H = ZKZ' + λ I where λ = σ^2_e / σ^2_u.
P3D	When P3D = TRUE, variance components are estimated by REML only once, without any markers in the model. When P3D = FALSE, variance components are estimated by REML for each marker separately.
optimizer	The function used in the optimization process. We offer "optim", "optimx", and "nlminb" functions.
eigen.G	A list with $values Eigen values $vectors Eigen vectors The result of the eigen decompsition of G = ZKZ'. You can use "spectralG.cpp" function in RAINBOW. If this argument is NULL, the eigen decomposition will be performed in this function. We recommend you assign the result of the eigen decomposition beforehand for time saving.
min.MAF	Specifies the minimum minor allele frequency (MAF). If a marker has a MAF less than min.MAF, it is assigned a zero score.
count	When count is TRUE, you can know how far RGWAS has ended with percent display.

Value

-log10(p) for each marker

References

Kennedy, B.W., Quinton, M. and van Arendonk, J.A. (1992) Estimation of effects of single genes on quantitative traits. J Anim Sci. 70(7): 2000-2012.

Kang, H.M. et al. (2008) Efficient Control of Population Structure in Model Organism Association Mapping. Genetics. 178(3): 1709-1723.

Kang, H.M. et al. (2010) Variance component model to account for sample structure in genome-wide association studies. Nat Genet. 42(4): 348-354.

Zhang, Z. et al. (2010) Mixed linear model approach adapted for genome-wide association studies. Nat Genet. 42(4): 355-360.

KosukeHamazaki/RAINBOW documentation built on Dec. 12, 2020, 8:35 p.m.