Displaying 2 results from an estimated 2 matches for "newdis2".
Did you mean:
newdisk
2012 Jul 11
0
declaring negative log likelihood of a distribution
...return(-KBLL1)
}*
#here a and b are the parameters to be estimated and x=binomial
values, fre=frequencies and n=binomial trials
#and since the inner series is a convergent infinite series, a large
number like 1000 is defined as the maximum
#with for loop
*>Loglik.newdis2<-function(x,a,b,fre,n,imax) {
term<-0
for (i in 0:imax) {
term=term+(((-1)**i)*(choose(b-1,i))*(beta(x+a+a*i,n-x+1)))
}
dens=a*b*choose(n,x)*term
KBLL2<-sum(fre*log(dens))
return(-KBLL2)
} *
#Now to estimate the param...
2012 Jul 17
0
Maximum Likelihood estimation of KB distribution
...782/kb.png
For a given data I need to estimate the paramters (alpha and beta) of this
distribution(Known as Kumaraswamy binomial Distribution, A Binomial Like
Distribution). For that, in order to use *optim()*, I first declared the
Negative Log-likelihood of this distribution as follows;
*Loglik.newdis2<-function(x,a,b,n,imax) {
term<-0
for (i in 0:imax) {
term=term+(((-1)**i)*(choose(b-1,i))*(beta(x+a+a*i,n-x+1)))
}
dens=a*b*choose(n,x)*term
KBLL2<-sum(log(dens))
return(-KBLL2)
} *
since there is an infinite convergent s...