使用nls函数求解指数截尾的幂律分布求助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

现在有一些随机数据，我想看它和哪个已知的分布拟合的最好，其中，带指数截尾的幂律分布（power law with exponential cutoff）是其中一个，它的概率密度函数为：

1. f(x;α,λ)=C*x^(−α)*exp(−λ*x)

复制代码

首先，我随机生成一些随机值来代替我的实际数据：

1. data <- rlnorm(1000,0.6,1.23)
   
2. h <- hist(data,breaks=1000,plot=FALSE)
   
3. x <- h$mids
   
4. y <- h$density

复制代码

然后，使用`nls`函数来进行参数评估：

1. nls(y~c*x^(-a)*exp(-b*x),start=list(a=1,b=1,c=1))

复制代码

进行到这里，总是会抛出两个错误中的其中一个：
1.奇异值矩阵；在网上搜索了，是初始值的原因的，尝试了几次，有时候能解决；
2.Error in numericDeriv(form[[3L]], names(ind), env) :   在计算模型的时候产生了缺省值或无限值；尝试做过一些数据处理，但是这个问题一直无法解决；
对于这个参数估计的问题，我应该怎么求解，或者是否有更好的方法或者函数来实现带指数截尾的幂律分布的参数估计，先谢谢大家了

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：函数求解 nls Exponential 概率密度函数 cutoff power 最好