我正在使用 R 中包 MASS 中的函数 glm.nb()。如何找到正确的 init.theta 值？在我看到的关于这个主题的其他问题中，提问者似乎已经有了他们的 init.theta。
我知道 theta 是一个色散参数。当我不输入 init.theta 时，摘要向我显示了高得离谱的 theta 和标准错误值。代码运行但给了我一条错误消息，说我达到了 theta.ml 迭代限制。
我进行了修改并添加了 init.theta = 1，这给了我较低的 theta、标准误差和 p 值。我仍然不明白这个论点是什么意思。
模拟我正在使用的代码：
glm.nb(measurement ~ year, data = df, init.theta = 1, link = log)