HelloI am trying to estimate the parameter of a function by the Maximum
Likelihood Estimation method.If the function is the difference between two
integrals: C<-function(n){integrand3<-function(x) {((2-x)^n)*(exp(ax-2))}cc<-
integrate (integrand3,0,2)print(cc)}
D<-function(n){integrand4<-function(x) {((5-ax)^n)}cc<- integrate
(integrand4,0,2)print(cc)}
f(n) = C(n) - D(n)
I need to estimate parameter (a). loglikelihood function is in the form of
sum[F(k) log(f(k))]=lnL. I am wondering how to introduce my logliklihood
function to the loop below. Can anybody help me for correcting the following
loop? if there are some other packages better than this , please let me
know.Thank you,Diba
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(C(n)-D(n)))
re <- maxLik (loglik, start=.5)
summary(re)
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
