本人最近在做极值统计方面的文章,想要用极值理论中的广义帕累托分布(GPD),对保险损失分布进行拟合。按照一般文献中经常使用的平均剩余函数图、GPD指数图等,对门限值(阈值)进行大致的选取,利用R的极值软件包对比例参数(Scale)、形状参数(Shape)进行了极大似然估计(MLE),但是总是通不过K-S拟合优度的检验,P-value非常小。
很多极值的教材,也都说通过看图来确定的门限值都是很主观的。我也知道存在通过一些非图示的方法来确定门限值的,请教高手,有没有合适的可以在R软件上实现确定门限值相对精确些的方法?