【新手求教】Error in optim(init, gpd.lik,optim回覆了无限值）

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我刚开始用R求数据GPD分布，按照史道济《极值理论及应用》方法学习
数据有缺失，导入数据时应该删除了吧
temp<-read.table("EV.txt",header=T,na.string="")
得到平均寿命图，大致判定阈值0为正斜率
然后，求阈值范围内的参数估计，就出现问题了
>gpd.fitrange(T1.Idr,-2,0.5,nint=50)
Error in optim(init, gpd.lik, hessian = TRUE, method = method, control = list(maxit = maxit,  :
  optim回覆了无限值
> T1.gpd<-gpd.fit(T1.Idr,0,npy=100)
Error in optim(init, gpd.lik, hessian = TRUE, method = method, control = list(maxit = maxit,  :
  optim回覆了无限值

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：数据分析 GPD

已解决，过程如下：
> T1.Idr<-log(T1[2:1035])-log(T1[1:1034])
Warning messages:
1: In log(T1[2:1035]) : 产生了NaNs  
我发现数据T1数组有负数，对数无法计算，导致新变量T1.Idr包含异值，
所以后面新变量无法正确操作。
把负数据删除，之后就有结果了。

z1075328886 发表于 2020-4-16 21:30
已解决，过程如下：
> T1.Idr

你好，同样的问题，具体您是如何解决的呢