2

我正在尝试使用 step 函数拟合准泊松 GLM 和 GAM 来执行向后选择。我使用了以下代码,但 R 返回了一个错误,因为它无法计算 AIC。使用泊松分布的相同模型有效,但不是准泊松,有没有办法改变函数使其有效?

full.mod03q<-glm(Total.Positivo_2003~X.baix03 + Tot_PrevDesnut04
                 + Tot_HosDesid04 + Despesa_06 + prpnurb03 + prpnkid03
                 + dev03 + watertoilet03 + elec03 + dollpercap03
                 + abspov03 + gini03 + Importado.Outro.Municipio.Da.Uf_2003 
                 + Importado.Uf.Pais_2003 + elevation + floodplain + BIO1 
                 + BIO6 + BIO12 + deforestation_prop,
                 family="quasipoisson", offset=log(pop_2003), 
                 data=ama1@data)
step(full.mod03q, direction="backward")
4

1 回答 1

0

AIC 不能在准泊松模型中计算。函数“step”使用 AIC 选择更多相关变量以包含在模型中。有一个计算“准AIC”的“非纯粹”解决方案。

我正在与同样的问题作斗争。它似乎可以在以下位置解决:

https://cran.r-project.org/web/packages/bbmle/vignettes/quasi.pdf

于 2017-01-27T13:58:41.580 回答