predict6_9_unique: Produce model estimated (p=6, topology=9, unique portion)...

View source: R/model_functions.R

predict6_9_uniqueR Documentation

Produce model estimated (p=6, topology=9, unique portion) y-coordinates of the kmer spectra given the kmer size, repetitiveness, average polyploid kmer coverage, bias, and x-coordinates of the kmer spectra.

Description

Produce model estimated (p=6, topology=9, unique portion) y-coordinates of the kmer spectra given the kmer size, repetitiveness, average polyploid kmer coverage, bias, and x-coordinates of the kmer spectra.

Usage

predict6_9_unique(raaaabb, raaabcc, raabcdd, raabcde, rabcdef, k, d,
  kmercov, bias, x)

Arguments

raaaabb, raaabcc, raabcdd, raabcde, rabcdef

Numerics corresponding to the nucleotide heterozygosities aaaabb, aaabcc, aabcdd, aabcde, and abcdef respectively.

k

An integer corresponding to the kmer length.

d

A numeric corresponding to the repetitiveness.

kmercov

A numeric corresponding to the estimated average kmer coverage of the polyploid genome.

bias

A numeric corresponding to the overdispersion of the negative binomial distribution.

x

An integer vector of the x-coordinates of the histogram (after filtering out low coverage errors and high coverage kmers).

Value

A numeric vector of the model estimated y-coordinates of the kmer spectra.


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