MBttest (version 1.0.0)

betaparametab: Estimation of Beta Parameters a And b

Description

parameters alpha(a) and beta (b) in betat distribution are estimated by using modified Baggerly et al(2003)'s iterative optimal method.

Usage

betaparametab(xn, w, P, V)

Arguments

xn

column vector, a set of library sizes.

w

column vector, a set of weights

P

proportion of counts of a gene or an isoform

V

variance for proportions of counts of a gene or an isoform over m replicate libraries in a condition

Value

return parameters a and b.

References

Baggerly KA, Deng L, Morris JS, Aldaz CM (2003) Differential expression in SAGE: accounting for normal between-library variation. Bioinformatics 19: 1477-1483.

See Also

betaparametVP, betaparametw

Examples

XX<-c(2000,2000,2000)
p<-0.15
V=0.004
w<-c(0.3,0.3,0.3)
betaparametab(xn=XX,w=w,P=p,V=V)
#[1] 1.145868 6.493254