fitKumBin: Fitting the Kumaraswamy Binomial Distribution when binomial...

Description Usage Arguments Details Value References See Also Examples

View source: R/Kumaraswamy.R

Description

The function will fit the Kumaraswamy binomial distribution when random variables, corresponding frequencies and shape parameters are given. It will provide the expected frequencies, chi-squared test statistics value, p value, degree of freedom and over dispersion value so that it can be seen if this distribution fits the data.

Usage

1
fitKumBin(x,obs.freq,a,b,it,print)

Arguments

x

vector of binomial random variables

obs.freq

vector of frequencies

a

single value for shape parameter alpha representing a

b

single value for shape parameter beta representing b

it

number of iterations to converge as a proper probability function replacing infinity

print

logical value for print or not

Details

0 < a,b

x = 0,1,2,...n

obs.freq ≥ 0

it > 0

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

Value

The output of fitKumBin gives a list format consisting

bin.ran.var binomial random variables

obs.freq corresponding observed frequencies

exp.freq corresponding expected frequencies

statistic chi-squared test statistics

df degree of freedom

p.value probability value by chi-squared test statistic

over.dis.para over dispersion value.

References

Li, X. H., Huang, Y. Y., & Zhao, X. Y. (2011). The Kumaraswamy Binomial Distribution. Chinese Journal of Applied Probability and Statistics, 27(5), 511-521.

See Also

mle2

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
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(EstMLEKumBin,start = list(a=10.1,b=1.1,it=10000),
          data = list(x=No.D.D,freq=Obs.fre.1)))
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.
fitKumBin(No.D.D,Obs.fre.1,aKumBin,bKumBin,itKumBin*100)

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