Displaying 1 result from an estimated 1 matches for "estimates1".
Did you mean:
estimates
2012 Jul 11
0
declaring negative log likelihood of a distribution
...}
dens=a*b*choose(n,x)*term
KBLL2<-sum(fre*log(dens))
return(-KBLL2)
} *
#Now to estimate the parameters, I used
http://rgm2.lab.nig.ac.jp/RGM2/func.php?rd_id=bbmle:mle2 mle2(from the
package bbmle) , which is a wrapper around the optim function
*>estimates1=mle2(Loglik.newdis1, start=list(a=1,b=1),
data=list(x=values,fre=frequency, n=7, imax=2000)) *
Error in optim(par = c(1, 1), fn = function (p) :
non-finite finite-difference value [1]
In addition: There were 50 or more warnings (use warnings() to see the first
50)
/#this gives the above error...