powered by
This function tunes the final fitting for heterozygous genomes by adjusting the delta values for heterozygous and homozygous regions.
error_minimize2(tooptimize, h_het, h_hom, h_target)
A numeric value representing the minimized difference.
A numeric vector containing the scale factors to optimize.
A numeric vector representing the raw fitting for the heterozygous region.
A numeric vector representing the raw fitting for the homozygous region.
A numeric vector representing the target k-mer frequency.
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)
Run the code above in your browser using DataLab