我尝试使用阶跃函数拟合准泊松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")
答案 0 :(得分:0)
无法在准泊松模型中计算AIC。功能&#34;步骤&#34;使用AIC选择更多相关变量以包含在模型中。有一个&#34;不是纯粹主义者&#34;解决方案计算&#34;准AIC&#34;。
我正在反对同样的问题。看来可以解决:
https://cran.r-project.org/web/packages/bbmle/vignettes/quasi.pdf