现在有一些随机数据,我想看它和哪个已知的分布拟合的最好,其中,带指数截尾的幂律分布(power law with exponential cutoff)是其中一个,它的概率密度函数为:
首先,我随机生成一些随机值来代替我的实际数据:
然后,使用`nls`函数来进行参数评估:
进行到这里,总是会抛出两个错误中的其中一个:
1.奇异值矩阵;在网上搜索了,是初始值的原因的,尝试了几次,有时候能解决;
2.Error in numericDeriv(form[[3L]], names(ind), env) : 在计算模型的时候产生了缺省值或无限值;尝试做过一些数据处理,但是这个问题一直无法解决;
对于这个参数估计的问题,我应该怎么求解,或者是否有更好的方法或者函数来实现带指数截尾的幂律分布的参数估计,先谢谢大家了