r - 具有潜在警告/错误和使用 tryCatch 的 R 全局优化问题
问题描述
简单地说,我有一个函数f(x, t2)
,我想找到关于的x
积分最大化的值。我选择 pso 算法进行优化。可执行代码如下f(x, t2)
t2
library(pso)
xl=0; xu=2000; n=1; t2l=100; t2u=2000; t1=1
g<-function(x, t2) t1*x/(t2+x)
h<-function(z) 1/z^n
gdot<-function(x, t2){
c(x/(t2+x),-t1*x/(t2+x)^2)
}
logdetHinv<-function(dp, dw, t2){
gmat=mapply(function(x) gdot(x,t2),dp)
D0=gmat%*%diag(dw)%*%t(gmat)
D1=gmat%*%diag(1/h(g(dp,t2)))%*%diag(dw)%*%t(gmat)
2*log(det(D1))-log(det(D0))
}
obj<-function(x){
dp=x[1:2]; dw=c(x[3],1-x[3])
fitness_value=-integrate(Vectorize(function(t2) logdetHinv(dp, dw, t2)*1/(t2u-t2l)), t2l, t2u)$value
return(ifelse(dw[2]>0, fitness_value, fitness_value+1e3))
}
x <- psoptim(rep(1,3), fn = obj, lower = c(rep(xl,2),0.1), upper = c(rep(xu,2), 0.9))$par
x
因为全局优化涉及一些随机过程,它有时会报告正确的结果
> x
[1] 2000.0000 754.4146 0.5000
其他时候它报告错误
Error in integrate(Vectorize(function(t2) logdetHinv(dp, dw, t2) * 1/(t2u - :
non-finite function value
In addition: There were 11 warnings (use warnings() to see them)
> warnings()
Warning messages:
1: In log(det(D1)) : NaNs produced
2: In log(det(D0)) : NaNs produced
3: In log(det(D1)) : NaNs produced
4: In log(det(D0)) : NaNs produced
我想该算法会尝试记录 中的一些负值logdetHinv
,它会返回NaN
一条警告消息,而不是错误,最后会导致 中的错误integrate
。
我想避免这样的值,也许有tryCatch
,比如如果函数中有警告logdetHinv
,它返回一个非常小的值,但不是NaN
,所以不会导致错误integrate
,并且psoptim
在最大化目标函数时不太可能选择这样的值(最小化-integrate(logdetHinv)
)。我不熟悉tryCatch
这种复杂的情况。我应该把它放在哪里tryCatch
?谢谢。
此外,我想知道 R 中是否有一些调试技术可以让我知道D0/D1
在这种情况下是什么随机值 () 导致了错误。我猜它是对数中的一些负值,但它不应该,因为对数内部是正定矩阵的行列式。在回溯模式下,在浏览中,如果我键入D0
对象“D0”将找不到。
解决方案
在这种情况下,我不会使用在测试中通常比在您的主代码中更合适的 tryCatch。为什么不简单地测试函数中的决定因素?像这样的东西应该工作:
logdetHinv<-function(dp, dw, t2){
gmat=mapply(function(x) gdot(x,t2),dp)
D0=gmat%*%diag(dw)%*%t(gmat)
D1=gmat%*%diag(1/h(g(dp,t2)))%*%diag(dw)%*%t(gmat)
detD1 <- max(0.01, det(D1))
detD0 <- max(0.01, det(D0))
2*log(detD1)-log(detD0)
}
推荐阅读
- jquery - 如何在 jquery 事件中调用 vue 方法?
- sql-server - 在同一个表中插入 2 行的比率
- angular - 注销 MatDialog 在移动浏览器中不起作用
- java - ANTLR4:Java - TokenStreamRewriter getText() 方法的性能问题
- mark.js - Mark.js 突出显示包含在 html 标记中的文本
- apache - 如何在 Plack/Apache 中拦截回复
- reactjs - 用 react-image-gallery 反应垂直轮播
- python - sklearn线性回归中没有学习率和迭代次数
- .net - 我的 .net 项目中的错误程序集引用 [已解决]
- batch-file - 如何将 2 个批处理命令合并为一行