alphaBeta_to_meanVar: fit a 2-component beta mixture

In nickytong/SIBER: Systematic Identification of Bimodally Expressed Genes Using RNAseq Data

Description

can impose the two component to have equal variance (model='E') or unequal variance (model='V')

Usage

alphaBeta_to_meanVar(alpha, beta)

Arguments

y	ration values between 0 and 1
inits	Initial value to fit the mixture model. A vector with elements alpha1, beta1, alpha2, beta2 and pi1.
model	Character specifying E or V model. E model fits the mixture model with equal variance while V model doesn't put any constraint.

Value

A vector consisting parameter estimates of mu1, mu2, phi1, phi2, pi1, logLik and BIC. For 0-inflated model, mu1=phi1=0.

nickytong/SIBER documentation built on May 23, 2019, 5:08 p.m.