准泊松GLM的阶跃函数

时间:2014-07-22 14:39:08

标签: r glm gam

我尝试使用阶跃函数拟合准泊松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")

1 个答案:

答案 0 :(得分:0)

无法在准泊松模型中计算AIC。功能&#34;步骤&#34;使用AIC选择更多相关变量以包含在模型中。有一个&#34;不是纯粹主义者&#34;解决方案计算&#34;准AIC&#34;。

我正在反对同样的问题。看来可以解决:

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