flexdog where we now allow an outlier distribution.E-step in flexdog where we now allow an outlier distribution.
get_wik_mat_out(probk_vec, out_prop, refvec, sizevec, ploidy, seq, bias, od)The vector of current prior probabilities of each genotype.
The probability of being an outlier.
A vector of counts of reads of the reference allele.
A vector of total counts.
The ploidy of the species. Assumed to be the same for each individual.
The starting value of the sequencing error rate.
The starting value of the bias.
The starting value of the overdispersion parameter.
Same as get_wik_mat but the last column is
for the outlier class.
flexdog for the full EM algorithm.
get_wik_mat for the equivalent function
without outliers. doutdist for the outlier
density function.