NegLLBetaBin {fitODBOD}

Negative Log Likelihood value of Beta-Binomial Distribution

Description

This function will calculate the Negative Log Likelihood value when the vector of binomial random variables and vector of corresponding frequencies are given with the shape parameters a and b.

Usage

NegLLBetaBin(x,freq,a,b)

Arguments

x

vector of binomial random variables.

freq

vector of frequencies.

a

single value for shape parameter alpha representing as a.

b

single value for shape parameter beta representing as b.

Details

0 < a,b

freq \ge 0

x = 0,1,2,...

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

Value

The output of NegLLBetaBin will produce a single numeric value.

References

Young-Xu Y, Chan KA (2008). “Pooling overdispersed binomial data to estimate event rate.” BMC medical research methodology, 8, 1–12. Trenkler G (1996). “Continuous univariate distributions.” Computational Statistics and Data Analysis, 21(1), 119–119. HUGHES G, MADDEN L (1993). “Using the beta-binomial distribution to describe aggegated patterns of disease incidence.” Phytopathology, 83(7), 759–763.

Examples

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

NegLLBetaBin(No.D.D,Obs.fre.1,.3,.4)   #acquiring the negative log likelihood value