BBprob {BinaryEPPM}

Calculation of vector of probabilities for the beta binomial distribution.

Description

Given a vector of parameters and a scalar of the number of trials the function returns a vector of probabilities.

Usage

BBprob(twoparameter, nt)

Arguments

twoparameter

A vector of the parameters of the beta binomial distribution.

nt

The number of trials.

Value

Vector of probabilities

Author(s)

David M. Smith <[e-mail]>

References

Smith D (1982). Algorithm AS189. Maximum Likelihood Estimation of the Parameters of the Beta Binomial Distribution. Applied Statistics, 32, 196-204.

Williams D (1996). "Overdispersion in Logistic Linear Models." In B Mrgan (ed.), Statistics in Toxicology, pp75-84. Oxford Science Publications.

Examples

twoparameter <- c(0.96477815,0.7561417)
names(twoparameter) <- c('p','theta')
nt <- 37
BBprob(twoparameter,nt)