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方法。