fitAseNullMulti: fit QuASAR for a multiple samples

In piquelab/QuASAR: Quantitative Allele Specific Analysis of Reads

Description

joint geontype across samples and estimate sample error

Usage

fitAseNullMulti(ref, alt, eps = rep(0.1, ncol(ref)), log.gmat, max.it = 100,
  tol = 1e-16, fixGprior = TRUE, verbose = TRUE)

Arguments

ref	matrix of reference read count
alt	matrix alternate read count
eps	vector of starting values
log.gmat	log of genotype priors
max.it	maximum number of iterations for algorithm
tol	tolerance to assess convergance
fixGprior	logical to fix genotype priors across all steps of the algorithm
verbose	logical turn on reporting during algorithm iterations

Value

list of genotypes, log genotypes, vector oferror estimates, log-likelihood, and the sum of log likelihoods

piquelab/QuASAR documentation built on May 25, 2019, 7:14 a.m.