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
if (FALSE) {
#estimating the parameters using maximum log likelihood value and assigning it
parameters <- EstMLEKumBin(x=No.D.D,freq=Obs.fre.1,a=10.1,b=1.1,it=10000)
bbmle::coef(parameters) #extracting the parameters
aKumBin <- bbmle::coef(parameters)[1] #assigning the estimated a
bKumBin <- bbmle::coef(parameters)[2] #assigning the estimated b
itKumBin <- bbmle::coef(parameters)[3] #assigning the estimated iterations
#fitting when the random variable,frequencies,shape parameter values are given.
results <- fitKumBin(No.D.D,Obs.fre.1,aKumBin,bKumBin,itKumBin*100)
results
#extracting the expected frequencies
fitted(results)
#extracting the residuals
residuals(results)
}
Run the code above in your browser using DataLab