run_unc_nls: Run Unconstrained NLS Regression

View source: R/run_unc_nls.R

run_unc_nlsR Documentation

Run Unconstrained NLS Regression

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)


kegrinde/STEAM documentation built on June 30, 2024, 11:27 p.m.