> fd$fr.gpd<-gpd.fit(fd$fr,3.5,npy=250)
能得出参数估计的值,但是有错误提示:
错误于`$<-.data.frame`(`*tmp*`, "fr.gpd", value = list(trans = FALSE, :
replacement has 16 rows, data has 975
接下来画诊断图,
> gpd.diag(fd$fr.gpd)
出现错误提示:错误于if (z$trans) { : 变元长度为零
最重要的是这次没有结果,图出不来
恳请大侠指点