zscore: Z-score of relatedness

In genostats/gaston.pop: Gaston package extension for genetic population applications

Description

Compute the Z-score for given relatedness

Usage

 zscore(app_pairs, m, sd) 

Arguments

app_pairs	A vector of relatedness values
m	Reference mean used in Z-score
sd	Reference standard deviation used in Z-score

Details

The Z-score is defined as Z=\frac{m_{pairs}- m}{sd} where m_{pairs} is the mean of relatedness values.

Value

Value of Z-score.

Author(s)

Claire Dandine-Roulland

See Also

app.pairs, app.unpairs

Examples

# load chr2 data set (~10k SNPs in low LD)
x <- read.bed.matrix( system.file("extdata", "chr2.bed", package="gaston") )

# Compute Rousset Relationship Matrix
K <- Rousset(x)$Rousset

# Compute Z-score relatedness between individuals from similar population compared to all pairs
zscore(app.pairs(K, x@ped$population), m=app.unpairs(K, x@ped$population, mean), sd=app.unpairs(K, x@ped$population, sd))

genostats/gaston.pop documentation built on Jan. 17, 2022, 10:58 a.m.