scale_kinship {qtl2}

Scale kinship matrix

Description

Scale kinship matrix to be like a correlation matrix.

Usage

scale_kinship(kinship)

Arguments

kinship

A kinship matrix, or a list of such in the case of the "leave one chromosome out" method, as calculated by calc_kinship().

Details

We take c_ij = k_ij / √(k_ii k_jj)

Value

A matrix or list of matrices, as with the input, but with the matrices scaled to be like correlation matrices.

Examples

grav2 <- read_cross2(system.file("extdata", "grav2.zip", package="qtl2"))
map <- insert_pseudomarkers(grav2$gmap, step=1)
probs <- calc_genoprob(grav2, map, error_prob=0.002)
K <- calc_kinship(probs)
Ka <- scale_kinship(K)