lpsolve不在R中打印lp模型

Question

我是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

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

Answer 1

使用write.lp打印出更大的LP。

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

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

Answer 2

我可以用一个简单的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看来，打印函数的所有附加参数都将被忽略：

用法

## S3方法用于类'lpExtPtr'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