depmix 函数来拟合两个状态伽马分布

问题描述

我在 R 中使用 depmixS4 包。我有一个看起来像伽马分布的数据，我假设有两种状态。我想将两态伽马分布拟合到 R 中的数据中。以下是我的代码：

mod <- depmix(freq ~ 1,data = mod.data,nstates = 2,family = gamma()) # use gamma 
fit.mod <- fit(mod)

但是，似乎我有一个错误，因为我没有在 family = gamma() 部分传递参数。如果我只使用 family = gaussian() 就可以了。有人可以帮我解决这个问题吗？谢谢！

解决方法

使用 gamma 分布可能会提供不合适的起始值问题，因此在使用 gamma 系列时一定要提供合适的值。