obj_parent_reparam: A reparameterization of 'obj_parent'

Description Usage Arguments Author(s)

View source: R/RcppExports.R

Description

A reparameterization of obj_parent

Usage

1
2
obj_parent_reparam(pcounts, psize, ploidy, pgeno, s, ell, r, weight = 1,
  outlier = FALSE, out_prop = 0.01)

Arguments

pcounts

The number of counts of the reference allele we observe in a parent.

psize

The total number of reads we observe in a parent.

ploidy

An integer. The ploidy of the species. This is assumed to be the same for each individual.

pgeno

The genotype of the parent.

s

Same as log(bias_val) in obj_offspring.

ell

We have seq_error = expit(ell) from obj_offspring.

r

Same as log((1.0 - od_param) / od_param) from obj_offspring.

weight

A double between 0 and 1. This should only not be 1 when doing the em algorithm and the weight is the probabiility that a point being ok.

outlier

A logical. Should we include an outlier model (TRUE) or not (FALSE)? Defaults to FALSE.

out_prop

The proportion of points that are outliers. Defaults (quite arbitrarily) to 0.01.

Author(s)

David Gerard


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