MixBetaBinomialFitStep: MixBetaBinomialFitStep

In gimelbrantlab/Qllelic: Detects differential allelic imbalance from RNA-seq data

MixBetaBinomialFitStep

Description

One step of fitting beta-binomial mixture distribution of AI in particular bin.

Usage

MixBetaBinomialFitStep(initials_old, coverage, observations)

Arguments

initials_old	Initials for EM step: initials = c(w1, alpha1, alpha2), weight of first component and alphas for both beta-binomial distributions in a mixture
coverage	A number, that represents the coverage bin
observations	A vector of "maternal counts" in the bin

Value

Re-fitted initials for next EM step.

gimelbrantlab/Qllelic documentation built on Dec. 30, 2024, 4:02 p.m.