使用Optim（）來估計的函數

Question

我想

loglike<-function (albe,m,n,res,resv){ 
t=n*m 
l=0 
for (i in 1:t) { 
    k=1 
    for (j in 1:m-1) { 

    k=k+albe[j]+albe[j+m]*resv[i]^2 
    } 
    l=l-k 
} 
for (i in 1:m-1) { 
    k=0 
    for (j in 1:n) { 
    k=k+albe[i]+albe[i+m]*res[j,]^2 
    } 
    l=l+k 
} 
loglike=l 
return(loglike) 
} 

因此，我使用Optim來估計在下面的函數的參數的參數（）

x<-runif(146) 
obj<-optim(fn=loglike, par=x, m=73, n=101,res=m$matrix, resv=m$vetor) 

RES是矩陣n×m的和RESV是長度爲nxm的矢量。

有人可以幫助我，給我下面的錯誤

Error in optim(fn = loglike, par = x, m = 73, n = 101, res = m$matrix, : 
    objective function in optim evaluates to length 0 not 1 

Answer 1

它不應該是僅僅

loglike 

代替

return(loglike) 

？