尽管在 R - 编程之家

如何解决尽管在 R

我正在使用包和函数 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)

非常感谢！