EstMGFBetaBin: Estimating the shape parameters a and b for Beta-Binomial...

Description Usage Arguments Details Value References See Also Examples

View source: R/Beta.R

Description

The functions will estimate the shape parameters using the maximum log likelihood method and moment generating function method for the beta-binomial distribution when the binomial random variables and corresponding frequencies are given

Usage

1
EstMGFBetaBin(x,freq)

Arguments

x

vector of binomial random variables

freq

vector of frequencies

Details

a,b > 0

x = 0,1,2,...

freq ≥ 0

NOTE : If input parameters are not in given domain conditions necessary error messages will be provided to go further

Value

The output of EstMGFBetaBin will produce a list format consisting

a shape parameter of beta distribution representing for alpha

b shape parameter of beta distribution representing for beta

References

Young-Xu, Y. & Chan, K.A., 2008. Pooling overdispersed binomial data to estimate event rate. BMC medical research methodology, 8(1), p.58.

Available at: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=2538541&tool=pmcentrez&rendertype=abstract .

Trenkler, G., 1996. Continuous univariate distributions. Computational Statistics & Data Analysis, 21(1), p.119.

Available at: http://linkinghub.elsevier.com/retrieve/pii/0167947396900158 .

Hughes, G., 1993. Using the Beta-Binomial Distribution to Describe Aggregated Patterns of Disease Incidence. Phytopathology, 83(9), p.759.

Available at: http://www.apsnet.org/publications/phytopathology/backissues/Documents/1993Abstracts/Phyto_83_759.htm

See Also

mle2

Examples

1
2
3
4
5
6
7
8
9
No.D.D=0:7        #assigning the random variables
Obs.fre.1=c(47,54,43,40,40,41,39,95)   #assigning the corresponding frequencies
#estimating the parameters using maximum log likelihood value and assigning it
parameters=suppressWarnings(bbmle::mle2(EstMLEBetaBin,start = list(a=0.1,b=0.1),
data = list(x=No.D.D,freq=Obs.fre.1)))
bbmle::coef(parameters)   #extracting the parameters

#estimating the parameters using moment generating function methods
EstMGFBetaBin(No.D.D,Obs.fre.1)

Amalan-ConStat/R-fitODBOD documentation built on Oct. 1, 2018, 7:13 p.m.