具有多项式项的ANCOVA的事后测试

Question

我有一个在多项式协变量和因子之间具有相互作用项的模型。我进行了ANCOVA测试，并且交互作用意义重大。现在，我想执行事后测试，但是我不知道该怎么做。我想测试不同物种的生长曲线是否存在显着差异（我有3种）。 我的模型是这样的：

model <- lm(Sepal.Length ~ poly(Petal.Length,2)*Species, data=iris)