dmvnorm结果出现-Inf

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我是在对200条曲线进行平均曲线的拟合，用optim进行优化，总是报出optim() non-finite finite-difference value 。一步步看完优化的函数，发现，优化的函数中用到了dmvnorm，dmvnorm出现了-Inf导致的。
我的dmonorm
Lfy<-function(par,areadf,weightdf){
  ind_N<-dim(weightdf)[2]; Lfy_ls<-c()
  for(k in 1:ind_N){
    area<-areadf[,k];weight<-weightdf[,k]
    mu<-NP_get_mu1(par,weight); #得到5个数字组成的向量
    SAD<-SAD1(par,weight);#得到5x5的对称矩阵
    library(mvtnorm);
    Lfy_ind<-dmvnorm(area,mu,SAD,log=TRUE)
    Lfy_ls<-c( Lfy_ls,Lfy_ind)#由ind_N个数据组成的向量，由于会出现个别的-Inf导致return也是Inf
    }
  Lfy<--sum(Lfy_ls)
  return(Lfy)
}



我想知道，dmvnorm为什么会算出-Inf？


非常感谢各位的解惑！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Norm NOR ORM difference finite