R 使用一般线性模型（GLM）后的事后测试

我想看看3组（klasse）和他们的攻击性行为（通过计算攻击性互动的数量）之间是否存在显著差异。情况与此无关，但必须加以考虑。我用一个一般的线性模型比较了两组之间的亲缘关系，并用泊松族作为计数数据

这是Rstudio中的数据输入：

data.frame(agressiekrab=c(8,5,1,10,6,12,4,17,1,1,5,9,11,2,3,0,21,17,4,1,10,4,14,15,22,8,19,0,6,16,4,10), klasse=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3), situatie=c(1,1,1,1,1,1,2,2,2,2,2,3,3,3,3,3,1,1,1,2,2,2,1,1,1,1,1,2,2,2,2,2))

structure(list(agressiekrab = c(8, 5, 1, 10, 6, 12, 4, 17, 1, 
1, 5, 9, 11, 2, 3, 0, 21, 17, 4, 1, 10, 4, 14, 15, 22, 8, 19, 
0, 6, 16, 4, 10), klasse = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3
), situatie = c(1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 3, 3, 3, 3, 
3, 1, 1, 1, 2, 2, 2, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2)), class = "data.frame", row.names = c(NA, 
-32L))

agmodel1<-glm(agressiekrab~klasse*situatie,poisson) 
summary(agmodel1)

这是非参数数据，所以不是正态分布 这是我在Rstudio中使用的代码：

data.frame(agressiekrab=c(8,5,1,10,6,12,4,17,1,1,5,9,11,2,3,0,21,17,4,1,10,4,14,15,22,8,19,0,6,16,4,10), klasse=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3), situatie=c(1,1,1,1,1,1,2,2,2,2,2,3,3,3,3,3,1,1,1,2,2,2,1,1,1,1,1,2,2,2,2,2))

structure(list(agressiekrab = c(8, 5, 1, 10, 6, 12, 4, 17, 1, 
1, 5, 9, 11, 2, 3, 0, 21, 17, 4, 1, 10, 4, 14, 15, 22, 8, 19, 
0, 6, 16, 4, 10), klasse = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3
), situatie = c(1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 3, 3, 3, 3, 
3, 1, 1, 1, 2, 2, 2, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2)), class = "data.frame", row.names = c(NA, 
-32L))

agmodel1<-glm(agressiekrab~klasse*situatie,poisson) 
summary(agmodel1)

现在我想使用一个事后测试来了解de组之间的差异（klasse）。我尝试了配对威尔科克斯测试，但它没有给我两对两的比较。我不知道如何解释结果。 我还尝试了TukeyHSD测试，但也没有成功（我认为这是不对的，因为数据是非参数的）

tukey方法给出了一个错误：

TukeyHSD(aov(agressiekrab~klasse*situatie))    
  Error in TukeyHSD.aov(aov(agressiekrab ~ klasse * situatie)) : 
  no factors in the fitted model

---

你能用R代码提供你的例子吗，这样人们可以更容易地重现一切？如果TukeyHSD出现错误，错误是什么？您是如何使用该函数的？一般来说，如果你让我们看看你到底在做什么会更容易。好吧，我在我生成的结果中添加了所有代码。很好，谢谢。如果您还提供data.frame，它将变得更容易。您添加了一张图片，但是类似于

data.frame（agresiekrab=c（8,5,1），klasse=c（1,1,1），situatie=c（1,1,1））的东西会更容易些。您可以在data.frame上使用
dput
函数来获得可以粘贴到post中的输出。data.frame现在在post中