在matlab里如何进行GPD分布的参数估计？

我是论坛的新人，程序运行出了错误，这两天正好交论文，所以麻烦会的人，知道一下，非常感谢！！！
这是我在matlab里运行的程序。采用的最大似然估计
x = trnd(5,2000,1);            
q = quantile(x,.95);           
y = x(x>q) - q;                  
n = numel(y)                    
lowerBound = [-0.5 0.0];
startingGuess = [0.3 1.0];
这里为止都运行出来了。
paramEsts = mle(y, @gpnegloglike, startingGuess, lowerBound);
kHat= paramEsts(1)   
sigmaHat= paramEsts(2)
运行显示错误