请选择 进入手机版 | 继续访问电脑版
482 1

[问答] 一个最优化问题 [推广有奖]

  • 8关注
  • 2粉丝

讲师

33%

还不是VIP/贵宾

-

威望
0
论坛币
310 个
通用积分
0
学术水平
1 点
热心指数
1 点
信用等级
0 点
经验
15111 点
帖子
218
精华
0
在线时间
379 小时
注册时间
2012-10-31
最后登录
2020-8-31

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
请问各位,以下程序到底哪里出了问题,函数ln1 可以计算,但是最优化的时候就会出问题

ln1=function(lam1,h,data,beta){
  x1=data$x1; y1=data$y1; n1=length(y1)
  
  x2=data$x2; y2=data$y2; n2=length(y2)
  
  theta=n1/(n1+n2)
  
  z1=x1*as.vector(dnorm(y1,x1%*%beta,h)*(y1-x1%*%beta))
  
  part=1+z1%*%lam1/theta
  
  index=which(1+z1%*%lam1/theta<=1/n1)
  part[index]=-log(n1)-1.5+2*n1*part[index]-(n1*part[index])^2/2
  
  nln=-2*sum(log(part+1e-50))
  
  if( is.na(nln)|is.nan(nln)|is.infinite(nln)|nln>1e10) nln=1e10
  
  nln  
}


x1=as.matrix(rnorm(100,0,1)); x2=as.matrix(rnorm(100,0,1))
x1=cbind(1,x1); x2=cbind(1,x2)
y1=rnorm(100,0,1); y2=rnorm(100,0,1)
data=list(x1=x1,x2=x2,y1=y1,y2=y2)
beta=c(1,2); h=1; lam1=c(0.1,0.1)
ln1(lam1,h,data,beta)
lam=c(0,0)
nlminb(start=lam,ln1,h=h,data=data,beta=beta)

二维码

扫码加我 拉你入群

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

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

关键词:最优化问题 最优化 function length matrix

现在的论坛怎么老出现乱码或者特殊符号发不出来

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-3-29 21:43