run_unc_nls: Run Unconstrained NLS Regression

In GrindeLab/STEAMcpp_mirror:

Description

Estimate the number of generations since admixture (g) from the observed local ancestry correlation curves using non-linear least squares regression. More details in Grinde et al. (TBD).

Usage

run_unc_nls(lacorr, k1, k2, start.a = 0, start.b = 1, start.g = 10)

Arguments

lacorr	local ancestry correlation; data frame with 3 columns: theta = recomb frac; corr = correlation; anc = ancestry components (e.g., 1_1, 1_2, 1_3, 2_2, 2_3, 3_3)
k1, k2	ancestry component indices (number between 1 and K)
start.a	starting value for intercept; default = 0
start.b	starting value for slope; default = 1
start.g	starting value for g; default = 10

Value

A single number indicating the estimated number of generations since admixture.

Examples

run_unc_nls(lacorr = example_corr, k1 = 1, k2 = 1)
run_unc_nls(lacorr = example_corr, k1 = 1, k2 = 2)

GrindeLab/STEAMcpp_mirror documentation built on Dec. 31, 2020, 12:05 p.m.