使用glm.nb函数并使用AIC函数来比较模型可靠性。
示例代码:
# 假设x和y分别为预测变量和响应变量 model1 <- glm.nb(y ~ x, data = mydata) # 拟合模型1 model2 <- glm.nb(y ~ x + z, data = mydata) # 拟合模型2 AIC(model1, model2) # 比较两个模型的AIC值
注:mydata是包含x、y、z变量的数据框。