fitBeta: Fit Beta Mixture Model The function fits a two-component Beta...

Description Usage Arguments Value Examples

View source: R/Beta.R

Description

Fit Beta Mixture Model The function fits a two-component Beta mixture model. This function directly maximize the log likelihood function through optimization.

Usage

1
fitBeta(y, inits = NULL, model = "V")

Arguments

y

A vector representing the RNAseq raw count.

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 dispersion phi while V model doesn't put any constraint. Currently only V model is implemented.

Value

A vector consisting parameter estimates of alpha1, beta1, alpha2, beta2 and pi1, logLik and BIC.

Examples

1
2
3
4
set.seed(100)
y=c(rbeta(100,1,4),rbeta(200,4,1))
fitBeta(y=y)
SIBER(y, model='Beta')

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