# Regression within categorical variable in R

Examining the factors that affect dung beetle abundance, including type of dung (categorical; 4 species) and the mass of dung (continuous):

glm.nb(abundance ~ dung + mass)


But how can I test for differences in abundance due to the effect of mass within given dung type?

Thank you in advance!