极值理论中如果尺度参数alpha，位置参数beta和形状参数随着时间改变的非其次泊松过程，怎么用R估计参数？

有现成的函数吗？还是要自己变成求解最大似然函数？

求教育啊