MVP.K.VanRaden: Calculate Kinship matrix by VanRaden method

In rMVP: Memory-Efficient, Visualize-Enhanced, Parallel-Accelerated GWAS Tool

Calculate Kinship matrix by VanRaden method

Description

Build date: Dec 12, 2016 Last update: Dec 12, 2019

Usage

MVP.K.VanRaden(
  M,
  maxLine = 5000,
  ind_idx = NULL,
  cpu = 1,
  verbose = TRUE,
  checkNA = TRUE
)

Arguments

M	Genotype, m * n, m is marker size, n is population size
maxLine	the number of markers handled at a time, smaller value would reduce the memory cost
ind_idx	the index of effective genotyped individuals used in analysis
cpu	the number of cpu
verbose	whether to print detail.
checkNA	whether to check NA in genotype.

Value

K, n * n matrix

Examples

genoPath <- system.file("extdata", "06_mvp-impute", "mvp.imp.geno.desc", package = "rMVP")
genotype <- attach.big.matrix(genoPath)
print(dim(genotype))

K <- MVP.K.VanRaden(genotype, cpu=1)

rMVP documentation built on Sept. 11, 2024, 8:56 p.m.