eval_model: Evaluate distinct model forms, in order to resolve ambiguity...

In PauloMoraMartinez/Genomescope2.0_: Reference-free profiling of genomes

Evaluate distinct model forms, in order to resolve ambiguity of which peak is the homozygous peak

Description

Evaluate distinct model forms, in order to resolve ambiguity of which peak is the homozygous peak

Usage

eval_model(kmer_hist_orig, nls0, nls1, p, round, foldername, arguments)

Arguments

kmer_hist_orig	A data frame of the original histogram data (starting at 1 and with last position removed).
nls0, nls1	The nlsLM model objects to evaluate and compare.
p	An integer corresponding to the ploidy.
round	An integer corresponding to the iteration number (0, 1, 2, 3) for the fitting process.
foldername	A character vector corresponding to the name of the output directory.
arguments	A data frame of the user-specified inputs.

Value

A list (nls, nlsscore) where nls is the nlsLM model object (with some additional components) and nlsscore is the score (model RSSE) corresponding to the best fit (of the p forms).

PauloMoraMartinez/Genomescope2.0_ documentation built on April 5, 2022, 12:02 a.m.