在做完这一步后出现Warning且保存表格数值都为NA 但是其他几组数据却能够顺利跑出
for(i in 1:ncol(byD1)) {
u.left.thresh = findthresh(stre, n - k)
u.right.thresh = findthresh(stre, k)
u.gpd = gpd.tail(stre, upper = u.right.thresh, lower = u.left.thresh, plot = F)
bzD1[, i] = gpd.2q(byD1[, i], u.gpd)
}
Warning messages:
Maximum Likelihood Method for the GPD did not converge in: gpd.ml(sample = excess, location = 0)
Problem in .Fortran(if(!cmplx) "dqr" else "zqr",: subroutine dqr: Missing values in argument 1
Use traceback() to see the call stack
有没有S-Plus高手可以解答一下 谢谢


雷达卡


京公网安备 11010802022788号







