search for: mloglik

Hi, I would like estimate a model for function of production's Coob-Douglas using maximum likelihood. The model is log(Y)= beta[1]+beta[2]*log(L)+beta[3]*log(K). I tried estimate this model using the tools nlm ( ) and optim ( ) using the log-likelihood function below: > mloglik <- function (beta, Y, L, K) { + n <- length(Y) + sum ( (log(Y)- beta[1]-beta[2]*log(L)-beta[3]*log(K))^2)/2*beta[4]^2 + n/2*log(2*pi)+ n*log(beta[4]) + } Then I did estimates the parameters using nlm ( ) and optim ( ), but the estimates were very bad. I used these codes: > mlem &lt...