Question

我正在尝试从optimr软件包中运行optimr函数。我的代码是

library(optimr)

#____function to compute Sum of Square Errors 
SSE=function(par,T) {
    options(warn = -1)

    nr = dim(T)[1]   # number of redds (rows)
    S  = numeric(nr) # redd survival vector
    TC = par[1]    
    bt = par[2] 

    #--- temperature deviations dT=(Ti-TC)
    dT = sweep(T,1,TC,FUN='-') 
    dT = ifelse(is.na(dT)|dT<0,0,dT)

    h = bt*dT 

    for(i in 1:nr) S[i] = 1-prod(exp(-h[i,]))


    SSE = sum((S-mean(S))^2)  # same as var(S)*(nr-1)

   options(warn = 0)
}

#___ Data
Temp = matrix(round(rnorm(2000,12,1.5),2), ncol=20)  #fictious data used   
hist(Temp,col='grey60')

#____ minimize SSE via optimizer
o = optimr(par=c(8,.01),fn=SSE,T=Temp,method='Nelder-Mead') 
print(cbind(Tcrit=o$par[1],bt=o$par[2]),digits=3)

问题在于，无论我给命令提供什么初始参数值，它都将返回相同的初始值。我知道，如果起始值位于平坦的表面上，则这些值不会改变，但是无论我提供什么值，都将返回相同的值。有什么建议吗？

使用优化来最小化功能

0 个答案: