EBSeqHMM (version 1.6.0)

f0: Calculate the prior predictive distribution of the Beta-Negative Binomial model

Description

Calculate the prior predictive distribution of the Beta-Negative Binomial model

Usage

f0(Input, AlphaIn, BetaIn, EmpiricalR, NumOfGroups, log)

Arguments

Input

expression values

AlphaIn,BetaIn,EmpiricalR

The parameters estimated from last iteration of EM.

NumOfGroups

How many transcripts within each Ng group

log

If set as TRUE, the output will in log scale.

Value

output a numeric vector, each element shows the prior predictive probability of one gene/isoform

Details

Function f0() will calculate the Beta-Negative Binomial prior predictive probability for a given set of parameters

Examples

f0(matrix(rnorm(100,100,1),ncol=10), .5, .6,
           matrix(rnorm(100,200,1),ncol=10), 100, TRUE)