updog_iterate: Iterate between indiividual level estimates and joint...

Description Usage Arguments Author(s)

View source: R/bad_updog_ideas.R

Description

Iterate between indiividual level estimates and joint averaging.

Usage

1
2
updog_iterate(ocounts, osize, qarray, r1vec, r2vec, seq_error = 0.01,
  itermax = 1000, tol = 10^-4)

Arguments

ocounts

A vector of non-negative integers. The ith element is the number of reads of the reference allele in the ith child.

osize

A vector of positive integers. The ith element is the total number of reads for the ith child.

qarray

An array of proportions. Each dimension size is equal to the ploidy plus 1.

r1vec

The posterior probabilities of the genotypes of parent 1 conditional only on the sequence data from parent 1.

r2vec

The posterior probabilities of the genotypes of parent 2 conditional only on the sequence data from parent 2.

seq_error

A non-negative numeric. This is the known sequencing error rate. The default is to estimate this using data that are all approximately the reference allele.

itermax

The number of times to iterate between the two approaches.

tol

A positive numeric. The tolderance for the stopping criteria.

Author(s)

David Gerard


dcgerard/updogAlpha documentation built on May 14, 2019, 3:10 a.m.