R glm.nb 拟合负二项式广义线性模型

说明

系统函数 glm() 的修改包括对负二项式广义线性模型的附加参数 theta 的估计。

用法

glm.nb(formula, data, weights, subset, na.action,
       start = NULL, etastart, mustart,
       control = glm.control(...), method = "glm.fit",
       model = TRUE, x = FALSE, y = TRUE, contrasts = NULL, ...,
       init.theta, link = log)

参数

formula , data , weights , subset , na.action , start , etastart , mustart , control , method , model , x , y , contrasts , ...	glm() 函数的参数。请注意，这些不包括family 和offset(但可以使用offset())。
init.theta	theta 参数的可选初始值。如果省略，则在使用泊松 GLM 进行初始拟合后使用矩估计器。
link	链接函数。当前必须是 log 、 sqrt 或 identity 之一。

细节

使用交替迭代过程。对于给定的 theta，GLM 使用与 glm() 所使用的相同过程进行拟合。对于固定方式，theta 参数是使用分数和信息迭代来估计的。两者交替进行，直至两者收敛。 (估计 theta 时的交替次数和迭代次数由 glm.control 的 maxit 参数控制。)

设置trace > 0 跟踪交替迭代过程。设置 trace > 1 跟踪 glm 拟合，设置 trace > 2 跟踪 theta 的估计。

值

继承自 glm 和 lm 类的拟合模型对象 negbin 。该对象类似于 glm 的输出，但包含三个附加组件，即 theta 用于 θ 的 ML 估计，SE.theta 用于其近似标准误差(使用观察到的信息而不是预期信息)，twologlik 用于两次对数似然函数。

例子

quine.nb1 <- glm.nb(Days ~ Sex/(Age + Eth*Lrn), data = quine)
quine.nb2 <- update(quine.nb1, . ~ . + Sex:Age:Lrn)
quine.nb3 <- update(quine.nb2, Days ~ .^4)
anova(quine.nb1, quine.nb2, quine.nb3)

参考

Venables, W. N. and Ripley, B. D. (2002) Modern Applied Statistics with S. Fourth edition. Springer.

也可以看看

glm , negative.binomial , anova.negbin , summary.negbin , theta.md

有一个simulate方法。