betaparametab {MBttest} | R Documentation

## Estimation of Beta Parameters alpha and beta

### Description

parameters alpha(a) and beta (b) in beta distribution are estimated by using an iteration algorithm.

### 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.

### Author(s)

Yuan-De Tan tanyuande@gmail.com

### 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

[Package MBttest version 1.8.1 Index]