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