是否可以使用具有R中组合列的变量的泊松分布运行GLM?
我正在研究不同物种的影响,笼子的密度以及鸡蛋放置的日数以及孵化的卵数量,因此我将阴影和未孵化的柱子联系起来。我的数据是计数数据。代码适用于family = binomial
,但我想测试泊松是否是更好的模型。
我的代码如下:
attach(EggV)
density <- as.factor(Density)
day <- as.factor(Day)
Y <- cbind (Hatched, Unhatched)
model.pois <- glm(Y ~ Species + density + day, data = EggV, family = poisson)
但是一旦我运行代码就会给我一个错误:
x [good,,drop = FALSE]中的错误:(下标)逻辑下标太长
如果我只使用变量“Hatched”或“Unhatched”运行相同的代码,它可以工作,但这不足以进行数据分析。