问题描述
我正在使用包和函数 lpsymphony::lpsymphony_solve_LP(obj,mat,dir,rhs,types,max)
来解决我的线性规划问题。
现在我想证明参数 mat
中的限制很重要,因此我只用限制的一个子集解决了这个问题 - 只有 mat
的一些行用于获得解决方案。每次我重新执行时,解决方案都会改变。
问题:为什么没有重现性,我该如何解决这个问题?!
但是现在让我们看看我的代码:
# objective function
obj = rep(1,times = (196))
# other arguments
...
# first solution
set.seed(455)
res_231=lpsymphony_solve_LP(obj,mat[1:231,],dir[1:231],rhs[1:231],types = types,max = TRUE)
matrix(c(res_231$solution[29:56]),ncol=4,byrow=TRUE)
## [,1] [,2] [,3] [,4]
## [1,] 0 0 0 0
## [2,] 0 0 0 1
## [3,] 0 0 0 0
## [4,] 1 1 0 0
## [5,] 0 0 1 0
## [6,] 0 0 0 0
## [7,] 0 0 0 0
# second solution
set.seed(455)
res_231=lpsymphony_solve_LP(obj,] 0 1 0 0
## [5,] 1 0 0 0 # difference in [5,1] compared to first solution
## [6,] 0 0 1 0 # difference in [6,3] compared to first solution
## [7,] 0 0 0 0
以及会议信息的摘录
session_info()
## version R version 4.1.0 (2021-05-18)
## ...
## lpsymphony * 1.20.0 2021-05-19 [1] Bioconductor (R 4.1.0)
非常感谢!
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)