我试图在R中运行生物量数据(还原生物量和生殖生物量与营养生物量的比率)的GLM作为栖息地类型("hab")的函数,收集年份数据("年"),以及数据收集("网站").我的数据看起来很适合Gamma分布,但我有8个生物量零观测值(约800个观测值),因此模型不会运行.处理这个问题的最佳方法是什么？使用的另一个错误分布是什么？或者在我的零观察中添加一个非常小的值(例如.0000001)是否可行？
我的模型是:
reproductive_biomass<-glm(repro.biomass~hab*year + site, data=biom, family = Gamma(link = "log"))