MLE for Cobb Douglas的功能

Question

我试图估算MLE的cobb道格拉斯生产函数。我创建了我的日志功能：

myDD<-function(alpha,gamma,delta,sigma,theta){
+ n<-NROW(lk)
+ e<-ly+theta*yy-log(gamma)-alpha*(1-delta)*lk-alpha*delta*ll
+ DD<-((-n/2)*log(2*pi))-((n/2)*log(sigma))-((t(e)%*%e)/(2*sigma))
+ return(-DD)
+ }

并尝试运行命令：

myMLE<-mle2(myDD,start = list(alpha=1,gamma=1,delta=1,sigma=1,theta=0.001))

但是它说我没找到lk。我无法理解为什么以及在哪里指定lk。