我刚开始用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回覆了无限值