急求！怎么用极大似然估计伽马分布的两个参数，然后进行对数据进行拟合并检验

lww1993 发表于 2014-11-21 17:08
你的重复量大吗？不大的话我觉得可以删除后再做。

我用optim，不管怎么调初始参数，结果总会有警告信息说产生NANS，可是数据中确实有大量的有正有负，这个怎么处理？

qinglingxuanfen 发表于 2014-11-24 09:54
我的数据很多，怎么找重复数据呢？有什么命令可以直接提出重复数据及所在位置编号吗

你可以从第一个数开始，寻找与第一个数相同的数据，再从第二个数开始，寻找。如此直到所有数遍历。

qinglingxuanfen 发表于 2014-11-24 09:57
我用optim，不管怎么调初始参数，结果总会有警告信息说产生NANS，可是数据中确实有大量的有正有负，这个怎 ...

你既然是gamma分布，就不能出现数字为负的。你需要换其它的分布。

lww1993 发表于 2014-11-24 11:43
你可以从第一个数开始，寻找与第一个数相同的数据，再从第二个数开始，寻找。如此直到所有数遍历。

这样是要编个循环程序吧

lww1993 发表于 2014-11-24 11:46
你既然是gamma分布，就不能出现数字为负的。你需要换其它的分布。

因为有大量正值与负值，我拟合的有位置参数的伽马分布，我设置的初始参数是可以保证log（x-c）中x-c是正的，应该是迭代中产生了三次负值，可是就用的optim啊，求极大似然，函数也是没有错的，这个怎么办啊？

qinglingxuanfen 发表于 2014-11-24 14:52
这样是要编个循环程序吧

对 我没有找到什么好用的函数。

qinglingxuanfen 发表于 2014-11-24 15:02
因为有大量正值与负值，我拟合的有位置参数的伽马分布，我设置的初始参数是可以保证log（x-c）中x-c是正的 ...

为何不先做相应的平移，保证所有的数都是正数，再做gamma分布拟合呢？

lww1993 发表于 2014-11-24 16:30
为何不先做相应的平移，保证所有的数都是正数，再做gamma分布拟合呢？

我把所有的值都加了3，都变成了正值，拟合普通的两个参数的伽马分布，得到的形状参数和尺度参数与之前拟合的3个参数的分布中两个参数没有加减或者比例关系，这怎么办？是不是如果x服从伽马分布，x+3不一定服从伽马分布呢？

qinglingxuanfen 发表于 2014-11-24 17:49
我把所有的值都加了3，都变成了正值，拟合普通的两个参数的伽马分布，得到的形状参数和尺度参数与之前拟合 ...

正态分布满足线性性质，可是若x服从伽马，y=ax+b，则y-b服从伽马分布，这个系数可以，但我现在要给原始数据统一加一个数，采用新数据拟合的话，前后两者的分布不一致了，没办法代回去了

qinglingxuanfen 发表于 2014-11-24 19:01
正态分布满足线性性质，可是若x服从伽马，y=ax+b，则y-b服从伽马分布，这个系数可以，但我现在要给原始数 ...

而且我尝试了constrOptim函数，但是ui和ci针对这个问题怎么设定呢？