 # problem with evaluating variable

Hi,
Im stuck with this code, wasted like 5h trying to find solution, thats why Im here asking you for advice. My task is to
1.create loglikelihood function for Poisson distribution (from function with 2 variables)
2. with maxNR (package maxLik) estimate the lambda value
I think function is calculated how it should be, still result looks a bit odd (less concerning)
the worst part is that I cant evaluate the lambda with maxNR becouse i have some red text and I dont know why. I will be very grateful for any response becouse I feel helpless.

``````library("maxLik")
library("rgl")
x = c(4,4,2,4,5,3,1,1,8,3,2,3,1,4,3,4,3,1,2,2,2,1,2,2,2,3,2,3,2,1,1,2,5,2,3,3,2,2,1,2,4,1,1,
1,1,3,1,6,4,5,1,1,3,1,3,2,2,1,1,2)
N = length(x)

loglik= function(k,lambda){
val = N*(-lambda)+sum(k)*log(lambda)-sum(log(factorial(k))+N*(-lambda)))
return(val)
}
k = seq(from=0, to=7 , by= 0.1)
(g=factorial(k))
lambda= seq(from=0, to=8 , by= 0.1)
fval = matrix(0,nrow=length(k),ncol=length(lambda))

for(i in 1:length(k)) {
for(j in 1:length(lambda)) {
fval[i,j] = loglik(k[i],lambda[j])
}
}
persp3d(k,lambda,exp(fval),col="skyblue")
image(k,lambda,fval)
image(k,lambda,exp(fval)

lnL= function(parametry){
k=param
lambda=param
val = N*(-lambda)+sum(k)*log(lambda)-sum(log(factorial(k))+N*(-lambda))
return(val)
}
Estimate = maxNR(fn=lnL,start = 0,0)
summary(Estimate)
``````

you have one more closing parenthesis than than there is a need for after the final lambda

This topic was automatically closed 21 days after the last reply. New replies are no longer allowed.

If you have a query related to it or one of the replies, start a new topic and refer back with a link.