为什么' .QP_qpgen2'在quadprog中找不到?

时间:2017-07-11 15:21:31

标签: r quadprog

我最近拂去了一个脚本,该脚本从quadprog包中调用solve.QP(我目前的版本为1.5-5)。现在它会生成错误" object' .QP_qpgen2'找不到"。我不明白为什么。

这个对象不是由我创建的,而是由quadprog中的solve.QP函数创建的。 在Github Quadprog.R上有代码(第117行):

 res1 <- .Fortran(.QP_qpgen2,
               as.double(Dmat), dvec=as.double(dvec),
               as.integer(n), as.integer(n),
               sol=as.double(sol), lagr=as.double(lagr),
               crval=as.double(crval),
               as.double(Amat), as.double(bvec), as.integer(n),
               as.integer(q), as.integer(meq),
               iact=as.integer(iact), nact=as.integer(nact),
               iter=as.integer(iter), work=as.double(work),
               ierr=as.integer(factorized))

错误可以从解决问题的documentation代码中生成:

##
## Assume we want to minimize: -(0 5 0) %*% b + 1/2 b^T b
## under the constraints:      A^T b >= b0
## with b0 = (-8,2,0)^T
## and      (-4  2  0) 
##      A = (-3  1 -2)
##          ( 0  0  1)
## we can use solve.QP as follows:
##
Dmat       <- matrix(0,3,3)
diag(Dmat) <- 1
dvec       <- c(0,5,0)
Amat       <- matrix(c(-4,-3,0,2,1,0,0,-2,1),3,3)
bvec       <- c(-8,2,0)
solve.QP(Dmat,dvec,Amat,bvec=bvec)

如果有帮助,我正在使用R v3.4.1。

2 个答案:

答案 0 :(得分:5)

正如我的评论所述,R 3.4有一种注册外部例程的新方法。 Quadprog依赖于Fortran例程。 要解决这个问题,您需要使用当前的Rtools从R 3.4中的源代码构建软件包。您需要安装Rtools并进行设置(谷歌搜索将为您提供如何为您使用的任何系统设置Rtools的指南)。然后,转到quadprog包的CRAN页面并下载源文件quadprog.tar.gz。最后,运行命令

install.packages("PATH_TO_FILE/quadprog_1.5-5.tar.gz", repos = NULL, type="source", INSTALL_opts = "--merge-multiarch")

或者,你可以等几天。我敢肯定,CRAN上的软件包很快就会更新。

答案 1 :(得分:1)

我不断收到有关此问题的电子邮件:

使用packageDescription("quadprog")查看已安装的软件包的R版本。

如果版本为R 3.3.x(或更早版本),请使用update.packages(checkBuilt=TRUE)将您的版本更新为在R 3.4.x下构建的版本。