# Charger les données
inverts <- read.csv("inverts.csv")
str(inverts)

mlm1 <- lm(PLD ~ temp*feeding.type, data=inverts)
anova(mlm1) # toutes les variables sont significatives