4

我是 R 中 LP 建模的新手。我正在使用 lpSolveAPI。当我尝试一个带有两个决策变量的小示例并打印模型时，它会打印完整的模型。

library(lpSolveAPI)
lprec <- make.lp(nrow=0,ncol=2,verbose="full")

set.objfn(lprec,c(6,5))
lp.control(lprec,sense="max")

add.constraint(lprec,c(1,1),"<=",5)
add.constraint(lprec,c(3,2),"<=",12)

set.bounds(lprec,lower=c(0,0),columns = c(1,2))


RowNames <- c("A","B")
ColNames <- c("R1","R2")
dimnames(lprec) <- list(RowNames, ColNames)

print(lprec)

#   Model name: 

#            R1    R2        
#Maximize     6     5        
#A            1     1  <=   5
#B            3     2  <=  12
#Kind       Std   Std        
#Type      Real  Real        
#Upper      Inf   Inf        
#Lower        0     0

但是当我尝试使用 25 个决策变量的模型并添加一些约束后，如果我尝试打印模型，那么它只会说：

Model name: 
a linear program with 25 decision variables and 5 constraints

请建议如何显示更大的模型。

4

2 回答 2

2

用于write.lp打印较大的 LP。

 write.lp(lprec, filename="test.lp")

将生成一个文本文件，您可以使用任何文本编辑器对其进行检查。

于 2015-08-05T18:16:40.653 回答
0

我可以用一个简单的 LP 重现它：

library(lpSolveAPI)
lprec <- make.lp(nrow=0,ncol=25,verbose="full")
add.constraint(lprec, rep(1, 25), "<=", 1)
add.constraint(lprec, c(1, rep(0, 24)), "<=", 5)
print(lprec)
# Model name: 
#   a linear program with 25 decision variables and 2 constraints

?print.lpExtPtr看来 print 函数的所有附加参数都将被忽略：

用法

## 类 'lpExtPtr' 的 S3 方法 print(x, ...)

论据

x 一个 lpSolve 线性规划模型对象。

... 其他参数被忽略。

因此，您最好的选择可能是提取单个信息并输出它们。例如：

# Rows and columns
nr <- dim(lprec)[1]
nc <- dim(lprec)[2]

# Constraint matrix
sapply(1:(dim(lprec)[2]), function(x) {
  ret <- rep(0, nr)
  c <- get.column(lprec, x)
  ret[c$nzrow] <- c$column
  ret
})
#      [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19]
# [1,]    1    1    1    1    1    1    1    1    1     1     1     1     1     1     1     1     1     1     1
# [2,]    1    0    0    0    0    0    0    0    0     0     0     0     0     0     0     0     0     0     0
#      [,20] [,21] [,22] [,23] [,24] [,25]
# [1,]     1     1     1     1     1     1
# [2,]     0     0     0     0     0     0

# Constraint directions
get.constr.type(lprec)
# [1] "<=" "<="

# Right-hand sides
get.constr.value(lprec)
# [1] 1 5
于 2015-08-05T06:17:57.243 回答