scale_kinship: Scale kinship matrix
In rqtl/qtl2geno: Treatment of Marker Genotypes for QTL Experiments
Description
Usage
Arguments
Details
Value
Examples
Description
Scale kinship matrix to be like a correlation matrix.
Usage
1
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 / sqrt(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
1
2
3
4
5
grav2 <- read_cross2(system.file("extdata", "grav2.zip", package="qtl2geno"))
map <- insert_pseudomarkers(grav2$gmap, step=1)
probs <- calc_genoprob(grav2, map, error_prob=0.002)
K <- calc_kinship(probs)
Ka <- scale_kinship(K)
rqtl/qtl2geno documentation built on May 28, 2019, 2:36 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Description Usage Arguments Details Value Examples
Description
Scale kinship matrix to be like a correlation matrix.
Usage
1 | 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
|
Details
We take c_ij = k_ij / sqrt(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
1 2 3 4 5 | grav2 <- read_cross2(system.file("extdata", "grav2.zip", package="qtl2geno"))
map <- insert_pseudomarkers(grav2$gmap, step=1)
probs <- calc_genoprob(grav2, map, error_prob=0.002)
K <- calc_kinship(probs)
Ka <- scale_kinship(K)
|
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.