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
estimate <- EstMLEBetaBin(No.D.D,Obs.fre.1,a=0.1,b=0.1)
bbmle::coef(estimate) #extracting the parameters
#estimating the parameters using moment generating function methods
results <- EstMGFBetaBin(No.D.D,Obs.fre.1)
# extract the estimated parameters and summary
coef(results)
summary(results)
AIC(results) #show the AIC value
Run the code above in your browser using DataLab