Question

我正在尝试在R-中复制Excel Solver，这基本上是一个约束优化问题

我正在尝试最小化每次操作的费用，即总支出/总操作数等于以下功能并带有一些约束。

CPA功能：

(a+b+c+d)/((consta+(Baln(a)))+ (constb+(Bbln(b)))+(constc+(Bcln(c)))+(constd+(Bdln(d)))

其中未知变量为a，b，c，d和const*代表来自回归和B*展位的常数对于回归系数（因此它们是我所拥有的值）。

这是我试图最小化的简化填充函数：

(a+b+c+d)/ (((69.31*ln(a))+(14.885*ln(b))+(21.089*ln(c))+(9.934*ln(d))-(852.93))

约束：

a+b+c+d>=0    
a+b+c+d<=130000(total spend)
a<=119000 (maxa)
a>=272.56(mina)
b<=11000(maxb)
b>=2.04(minb)
c<=2900(maxc)
c>=408.16(minc)
d<=136800(maxd)
d>=55.02(mind)

我正在使用约束优化功能。我的代码如下：

g<-function(a,b,c,d) { (a+b+c+d)/((consta+(Balog(a)))+ (constb+(Bblog(b)))+ (constc+(Bclog(c)))+ (constd+(Bdlog(d)))) }

gb<-function(a) g(a[1], a[2], a[3],a[4])

A<-matrix(c(1,0,0,0,-1,0,0,0,0,1,0,0,0,-1,0,0,0,0,1,0,0,0,-1,0,0,0,0,1,0,0,0,-1,-1,-1,-1,-1,1,1,1,1),4,10)

B<- c(mina, -maxa, minb, -maxb, minc, -maxc, mind, -maxd,-totalspend, 0)

constrOptim(c(273,6,409,56),g,gb,A,B)

当我运行优化函数时，它声明我的参数出错（Error in ui %*% theta : non-conformable arguments）。我认为这是函数的梯度编码错误但我不确定。任何帮助表示赞赏。

ContrOptim函数 - 参数错误

0 个答案: