HiI have a probability mass function similar to?pr(N=n)=
integral(((2-x)^n)*(exp(ax-2))) - integral (((5-ax)^n)), ? both integrals are
defined over the interval(0,2) with respect to x.?I am going to estimate the
parameter (a) with method of maximum likelihood estimation. The loglikelihood
?is : ? ? ? ? ? ? ? [F log(pr(n))]=lnL??where F is the vector of observations
and (n) is the vector of input for ?the defined pmf?.?Can anybody suggest me the
fastest way of getting the MLE?I have tried this
program:n<-c(0,1,2,3,4,5,6,7,8)F<-c(0,0,1,3,5,7,8,11,10)loglik<-
function(a) {sum(F*log(pr(n)))}re<- maxlik(loglik, start=.5)summary(re)
I don't know how to define the probability mass function ( pr(n) ) in the
written program.I would appreciate for any help.Best Regards
[[alternative HTML version deleted]]