「遺漏值，其中爲TRUE/FALSE需要」GA包中的錯誤

當使用GA軟件包運行遺傳算法時，我收到「缺少需要TRUE/FALSE的值」錯誤。在下面我的代碼，like是用於健身的功能和X是長度的數值向量6.「遺漏值，其中爲TRUE/FALSE需要」GA包中的錯誤

library(GA) 
y<-c(46,38,49,55) 
W<-matrix(c(0,1,0,0,1,0,0,0,0,0,0,1,0,0,1,0),ncol=4) 
ei<-eigen(W)$values 
data<- matrix(c(14691,14150,14607,15140,5337,5386,5207,5178),ncol=2) 
     like=function(x) 
      { 
       co<-c(x[2],x[3]) 
       e<-x[4]*W%*%y-data%*%co-x[1] 
       L<-30*log(2)-15*log(2*pi)+log(prod(1-x[4]*Re(ei)))-15*log(x[5])+sum(log(pnorm(-x[6]*e/sqrt(x[5])))) 
       -0.5/x[5]*t(e)%*%e 
       return(L) 
      } 
    GA <- ga(type = "real-valued", 
      fitness = like, 
      min=c(-100,0,-50,-10,0.001,0), 
      max=c(100,100,0,10,100,10), 
      popSize = 100, maxiter = 200, 
      optim = TRUE)

結果和錯誤是如下：

GA | iter = 1 
Mean =  -Inf | Best = -2.870387 
Error in if (fmin > (sfactor * fave - fmax)/(sfactor - 1)) { : 
    missing value where TRUE/FALSE needed

來源

2016-12-23 ShuangXin

這是一個神祕的錯誤，但這個問題的一個很好的提示是在迭代1總結信息：

Mean = -Inf | Best = -2.870387

您函數的平均客觀值是負無窮大，這意味着某些迭代正在返回客觀值負無窮大。調試代碼的位標識，讓這個問題的樣本數據：

x <- c(-3.2, 5.7, -26.9, -6.1, 37, 9) 
like(x) 
# [1] -Inf

既然你最大化，你可以只用return(pmax(L, -1e6))更換return(L)封頂返回函數值在一個合理的下限，例如。那麼你會得到明智的結果：

GA | iter = 200

Mean = 102.2197 |最佳= 106.0328 |最終本地搜索= 106.0328

來源

2016-12-23 19:27:38 josliber

我改變sqrt（x [5]）與一個新的x [5]在定義函數like和它的作品。我認爲負無窮是由log（0）引起的。感謝你的幫助。 – ShuangXin

「遺漏值，其中爲TRUE/FALSE需要」GA包中的錯誤

回答

相關問題