error_minimize2 {findGSEP}

Tuning Final Fitting for Heterozygous Genomes

Description

This function tunes the final fitting for heterozygous genomes by adjusting the delta values for heterozygous and homozygous regions.

Usage

error_minimize2(tooptimize, h_het, h_hom, h_target)

Arguments

tooptimize

A numeric vector containing the scale factors to optimize.

h_het

A numeric vector representing the raw fitting for the heterozygous region.

h_hom

A numeric vector representing the raw fitting for the homozygous region.

h_target

A numeric vector representing the target k-mer frequency.

Value

A numeric value representing the minimized difference.

Examples

tooptimize <- c(0.5)
h_het <- rnorm(100)
h_hom <- rnorm(100)
h_target <- rnorm(100)
diff <- error_minimize2(tooptimize, h_het, h_hom, h_target)
print(diff)