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

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

optim函数命令初始参数怎么设置

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：极大似然估计 伽马分布 极大似然 似然估计 Optim

1.你先确定gamma分布参数的初始值para；
2.定义一个函数，例如f(x)，其是gamma分布相应的似然函数，取对数，再取负值，因为optim计算的是最小值；
3.optim(para,f).即可得到你要的估计值。
注：初始值的选择很重要。这个要根据你的实际背景。

lww1993 发表于 2014-11-18 19:27
1.你先确定gamma分布参数的初始值para；
2.定义一个函数，例如f(x)，其是gamma分布相应的似然函数，取对数 ...

谢谢大神！我还有个问题，这个optim对于有位置参数的即3个参数的伽马分布也可以拟合吗？那mle函数也是定义个函数求极大似然估计的，我应该选择哪个函数呢？非常感谢您！

mle中是调用了optim来计算相应的参数的。所以两个方法是都可以用的。

optim是可以估计gamma的参数的。

lww1993 发表于 2014-11-19 14:43
mle中是调用了optim来计算相应的参数的。所以两个方法是都可以用的。

optim是可以估计gamma的参数的。

非常感谢您！我在做单样本Kolmogorov-Smirnov 检验时，出现结，这个数据重复应该怎么处理呢？
Warning message:
In ks.test(x, "pnorm", mean(x), sd(x)) :
  ties should not be present for the Kolmogorov-Smirnov test

qinglingxuanfen 发表于 2014-11-20 15:44
非常感谢您！我在做单样本Kolmogorov-Smirnov 检验时，出现结，这个数据重复应该怎么处理呢？
Warning m ...

既然你要检验的是正态分布，你用shapiro.test吧。shapiro.test(x)就行

lww1993 发表于 2014-11-21 10:58
既然你要检验的是正态分布，你用shapiro.test吧。shapiro.test(x)就行

不是，我要做的是伽马分布拟合检验

qinglingxuanfen 发表于 2014-11-21 14:32
不是，我要做的是伽马分布拟合检验

只是用pnorm试了下是否是正态的，拟合伽马同样出现了上述警告信息出现了ties，这个应该怎么办啊？

qinglingxuanfen 发表于 2014-11-21 14:36
只是用pnorm试了下是否是正态的，拟合伽马同样出现了上述警告信息出现了ties，这个应该怎么办啊？

你的重复量大吗？不大的话我觉得可以删除后再做。

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

我的数据很多，怎么找重复数据呢？有什么命令可以直接提出重复数据及所在位置编号吗