updog_maximize: Estimate parental genotypes using offspring genotypes.

Description Usage Arguments Author(s)

View source: R/updog.R

Description

Estimate parental genotypes using offspring genotypes.

Usage

1
2
3
4
updog_maximize(ocounts, osize, qarray, r1vec, r2vec, pk, pival = 0.99,
  alpha = 0.1, beta = 0.1, est_fudge = TRUE, tol = 10^-4,
  itermax = 1000, update_geno = FALSE, update_pi = FALSE,
  update_beta = FALSE)

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 three-way array of proportions. The (i, j, k)th element is the probability of an offspring having k - 1 reference alleles given that parent 1 has i - 1 refrerence alleles and parent 2 has j - 1 reference alleles. Each dimension of the array is ploidy + 1.

r1vec

A vector of prior probabilities whose ith element is the prior probability that parent 1 has i -1 copies of allele A.

r2vec

A vector of prior probabilities whose ith element is the prior probability that parent 2 has i -1 copies of allele A.

pk

The parameter space. Should all be between 0 and 1.

pival

A fudge factor. The probability of not begin a mistake. Estimated if est_fudge = TRUE. Set to 0.99. Only used if est_fudge = TRUE.

alpha

The alpha parameter in the beta distribution of the outliers. Only used if est_fudge = TRUE.

beta

The beta parameter in the beta distribution of the outliers. Only used if est_fudge = TRUE.

est_fudge

A logical. Should we estimate the fudge factor TRUE or not FALSE?

tol

A positive numeric. The tolerance for the stopping criterion.

itermax

A positive integer. The maximum number of iterations to run in the EM.

update_geno

A logical. Update the parental genotypes?

update_pi

A logical. Update the mixing proporiton?

update_beta

A logical. Update the beta distribution?

Author(s)

David Gerard


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