我最近拂去了一个脚本,该脚本从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。
答案 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
下构建的版本。