一个有点复杂的函数求极值

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

目标函数是copula函数密度、两个边缘分布密度函数的乘积，三个密度函数的表达式都比较复杂,求其最大值。。求指点，用什么方法做比较好？代码如下，fn是目标函数，heq是等式约束。。u、v是边缘分布概率，parq、parw是边缘分布参数。。

自己试验过一般的非线性规划求极值，包括alabama::auglag函数、Rsolnp::solnp函数。但效果不理想。。

1. library(lmomco)
   
2. library(copula)
   
3. fn <- function(x){
   
4.   u <- par2cdf(x[1],parq)
   
5.   v <- par2cdf(x[2],parw)
   
6.   f <- -1*dCopula(c(u,v),cop)*par2pdf(x[1],parq)*par2pdf(x[2],parw)
   
7.   return(f)
   
8. }
   
9. 
   
10. heq = function(x){
    
11.   h <- rep(NA, 1);
    
12.   u = par2cdf(x[1],parq);
    
13.   v = par2cdf(x[2],parw);
    
14.   h[1]=pCopula(c(u,v),cop)-0.99
    
15.   return(h)
    
16. }

复制代码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Alabama Library Copula BRARY 非线性规划

这是个constrained opt问题，尝试一下这个package：nloptr
遇到有什么问题可以继续回帖