br: Block Relaxation for parameters estimation

In radmixture: Calculate Population Stratification

Description

This function is also used for estimating Q and F but faster than EM.

Usage

br(g, q, f, acc, max.iter, tol, model)

Arguments

g	Genotype matrix with dimensions n × p, where n is sample size and p is the number of SNPs.
q	Ancestry coefficient matrix with dimensions n × K, where n is sample size and K is the number of populations.
f	Minor allele frequency matrix with dimensions K × p, where K is the number of populations and p is the number of SNPs.
acc	a logical value indicating whether use quasi-Newton accelerated BR or not.
max.iter	If acc = T, max.iter must be set, the default is 3. max.iter should greater than 1.
tol	Tolerance, if acc = F, tolerance must be set, the default is 1e-4.
model	Choose which model you want to use. Supervised learning or unsupervised learning.

Value

Estimation results of q, f and the loglikelihood value of each iteration.

radmixture documentation built on May 2, 2019, 6:11 a.m.